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AGENDA 

Western Region Summer Meeting of Common 
Denver Hilton Hotel, Denver, Colorado 
July 6, 7, 8, 1966 



Wednesday, July 6 

General 1130 Session 
9:00 a.m. 



9:^5 a.m. 
11:15 a.m. 
1:30 p.m. 

3:30 p.m. 
kzOO p.m. 



Welcome 
IBM Announcements 

1130 Monitor System 

1130 Programming Techniques 

Program for Optical Systems 
Design 

1130 Project Control Systems 

1130 General Discussion Users 
Panel 



Parallel U30 Sessions 
11:15 a.m. 

1:30 p.m. 



Graphic Report Generator 

1130 STRESS - Structural 
Engineering System Solver 



Evening Session 
7:00 p.m. 

7:30 p.m. 



Hew Users Session 
Sound off Session 



General Session 
9:00 a.m. 



10:30 a.m. 



Thursday, July 7 



1130 Scientific Subroutine 
Package 

Solution of the Induction Well 
Logging Response Equation, 
G* Copland, Halliburton 
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General Session - cont. 
11:00 a.m. 

11:50 a.m. 

1:30 p.m. 
2:00 p.m. 

2:20 p.m. 
3:30 p.m. 

Parallel Session 

Education 
10:30 a.m. 

11:00 a.m. 
11:30 a.m. 

1:30 p.m. 



The Direct and Cyclic JacoM 
Methods with Fadeeva Correction 
Algorithm for a Real Matrix 
H. Fettis - J. Caslin, Wright 
Patterson AFB 

Luncheon Speaker - Ed Schwarz, 
Product Administrator - 1130 
System 

Autospot Post Processor Development 
D. Oliver, ACF Industries 

Design of ME-LTV Evaporators used 

in Sea Water Conversion 

D. Kays, Steams - Roger Corp. 

Development of a Public Works 
Engineering D.P. Center 
J. Hunter, L.A. County Dept. 
of Engineering 

FLBSPS AMD STOVE 

K. Jones , Colorado Dept. of Highways 



University of Mississippi 

Test Scoring Program 

R. Ross, U. of Mississippi 

Student Scheduling on the 1620 
G. Crumley, The Citadel 

Decoding and Use of the Special 1230 

Marking Code on 1620 

Without Column Binary 

C. Stallings, Jr., U. of Hew Mexico 

Whitman College 

Registration and Grading System 
G. Purcell, Whitman 0o.ilege 
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Thursday, July 7 - continued 



Parallel Session - cont. 
Education 
2:00 p.m. 

2:30 p.m* 

3:30 p. aw 
l*-:00 p.m. 



Whitman College Budget Forecasting 
G. Purcell, Whitman College 

University of Mississippi 
Floating Point Subroutines 
R. Ross, U. of Mississippi 

Austin College SPS - Improved 
D. Musser, Austin College 

Student D.P. System at 
Christian Brothers College 
Employing an IBM 1620 
J. Wegener, Christian Brothers 
College 



Structural 



l:lf5 p.m. 
2:15 P*m. 

3:30 p.m. 



Matrix Structural Analysis 

E. Cook, Wichita State University 

Design of Heavy Multi-Story 
Rigid Frame Industrial 
Frame Structures 
V. Arndt, Stearns Roger Corp. 

General Discussion 

Period (Structural and Civil) 



General Session 
9:00 a.m. 



9:30 a.m. 



10:30 a.m. 



Friday^ July 8 



Library , Magazines and a Computer 
G. Ahlborn, U. of Calif omia, 
Berkeley 

Accurate Solution of Systems of 

Linear Equations 

D. Musser, Austin College 

How a Programmer Should Write a 

Program for an Engineer 

L. Mahoney, R. W. Beck Associates 
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At the outset, I would like to say that this paper to be presented 
does not concern itself so much' with the answer solution to the induction 
logging response equation as ' it does with the organization of the problem 
itself. The problem is not new and, no doubt, has been set up and solved 
on numerous occasions, not only by ourselves but by others in the same 
field of study. About fourteen years ago, our Company rented a 6ok from 
IBM and some preliminary work was done. Like airplanes, the greater the 
model number, the better, faster and higher the computer, so soon we were 
obliged to try a 65O. The problem seemed to be more of a card handling 
project than a computer solution to equations. From the results of many 
months of computing, some studies were made and a logging tool was designed 
and introduced in the field. This tool, in modified form, is essentially 
the type used today. 

In the last several years the induction logging tool has been more 
widely used and accepted. The Engineering group now had available a 1620, 
20K machine with a 1622 card read punch. This computer was used to optim- 
ize the design of our logging tools. Many various multi-coil systems were 
investigated in this optimization program. 

Halliburton Company is primarily a service company to the petroleum 
industry. We perform various services the customer may require. These in- 
clude cementing, fracturing, chemical services as well as electrical log- 
ging of the formations. These may be either DC -potential contact logs, 
radioactive logs or induction type logs using AC currents. This paper is a 
discussion of the computer programming required to design induction logging 
tools. The next series of slides show the generalized form of the equation 
involved. 

For integration of the curve, let us examine the left position of 
the curve showing this stepwise incrementation. As you can see, the error 
is always positive, so the solution is to just reduce the increment and im- 
prove accuracy. Well, it helped, but not much! So, I tried something dif- 
ferent and I'm sure it's not original, but if you calculate every other step 
and multiply by two, as shown in the right side of the graph, it works nicely. 
Since we knew the answer and the accuracy desired, some idea was known as to 
how far out to infinity to go before the change in the incremented value be- 
came insignificant. This turned out to be from five to seven minutes per 
value timewise. 

Upon examination of the Response Equation, one finds the main time 
consuming portion of the problem will. be in the integration of the R values 
and Z values. Since we were still operating with the 1620 BC 20K, it was 
necessary to compute these tables and store them in card form and use manual 
look-up. First of all, we needed to know what steps to do the integration 
and how far we needed to carry out the calculations since we were going to 
infinity. I don't know any other way to do this but to start running and 
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and see how well it works. I would like to recognize the mathematicians 
present who have already figured out thst one can integrate the equation 
directly without resorting to long drawn-out incrementation techniques. 
This did give ma check and is quite direct, however, it didn't work - for 
me at least - when I tried to integrate the Z equation. So, I still had 
the problem of generating tables and a quick way of checking myself as I 
went along. It is somewhat discouraging to know the answer to the R equa- 
tion is 1.0 and take so long to prove it on a 1620. 

Another problem was how to set the program up so anyone could break 
into the problem and use the computer, since most work was to be done after 
regular hours. 

The program was arranged so the output card carried all the informa- 
tion required by the problem and. could be used to input data for the next value. 
If, say at night, when someone wanted to use the computer, all one had to do 
was stop the computation a:nd clear the punch. The last output card was added 
to the program deck and when one wanted to restart the problem, all that was 
necessary was to reload this deck. This required reloading the program but 
since it was a small deck, it required only a few minutes. In" case of card 
jams at night, the last card in the punch output hopper could be used to re- 
start the calculation with nothing lost. 

I found that the first card which I called the zero card was easily 
handled by a shorter program and therefore a separate program was used to 
calculate the first card which in turn started the second .continued program. 
Obviously, I did not wish to integrate any further than required in either 
the R or Z direction. The answers were periodically plotted to obtain the 
graph shown. You will note, at approximately R = 12, the curve has begun to 
flatte'n. 

Most of you are probably aware of the curve fit program in the 1620 
Users Group Library. This program was used to check the plot and it was found 
a curve fit could have had values of from k to 12 and 10,000. I assumed, for 
all practical purposes, the answer was 1.0 at 10,000. This worked out rather 1 
nicely and a hyperbolic curve was determined that gave us quite good correla- 
tion in the area, especially from 8 to 12. Some values were calculated and by 
slight adjustment of the resulting constants in the hyperbolic equation, I 
could match the curve to the accuracy required by the tables. These constants 
were also used in the problem solution, as I will discuss a little later. The 
R tables were calculated stepwise integration from 0. to 12. and by a hyperbolic 
equation from 12. to 60. A set of cards was assembled covering this range and 
listed on a i+07 for tabular look-up. 

The Z table calculations were handled in a similar manner as the R 
tables. As it turned out, it was easier to calculate the zero card separately. 
By using this card as a starting data for a program, I started calculating Z 
tables in a likewise fashion. When the computer got to Z = 1.0 it went into 
the error routine. After checking the problem, program and equation, I found 
that the equation vas discontinuous at Z = 1.0. Since the computer was unable 
to calculate this value, I simply told it what the value of R = 0.0 and Z « 
1.0 was, and continued on. Since I was unable to integrate the equation in 
respect to Z directly and after plotting some of the values as the calculation 
progressed, I decided' that the answer to the problem was .50 at infinity. 
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I decided to make a 3-dimensional graph of the vertical response equation us- 
ing swab sticks cut to lengths and a piece of 1/2" black plastic sheet for a 
base. The equation is quickly solved when you are solving for points directly, 
I wrote a program that calculated these points and found that about 20 minutes 
running time was all that was required to give us the model shown. I includ- 
ed the 1/2" dimension for thickness of the board and the answer so that the 
sticks could be cut directly. An interesting point here was a short program 
was written for our automatic drill press in the main production facility to 
drill the base. In this program, it turns out it is faster on this particu- 
lar machine to drill every fourth hole - back up to the beginning, increment 
a quarter unit and drill another set of holes k units apart. 

By now I had decided our accuracy adequate but I desired to reduce 
the time required for calculations . By examining the three dimensional form 
we obtained from the plot, I was able to estimate to some degree a better 
technique to speed up the incrementation. If you will note the bulk of the 
graph weight lies in the volume from R = U.O and Z = 2.0. One might be able 
to speed up the calculations by increasing the size of the incrementation 
after exceeding these limits. This was done and we saved about 25% in com- 
puter time over the R table calculations in this manner. I am sorry we had 
not made an R table three-dimensional graph to begin with. As in the R tables, 
we used the curve fit program and obtained different values for the constants 
in the hyperbolic equation. These were adjusted slightly and values were 
calculated out to Z = 60.0. 

So now I have a deck of cards for R values from 0.0 to 60.0 and Z 
values from 0.0 to 60.0. I also listed these so they could be used in tabu- 
lar look-up. 

• We decided to check our previous work done on induction logging tools. 
I wrote a program which used the typewriter to tell me which cards were needed 
to be looked-up from data storage. This amounted to values for GR and 12 
values for GB. We would load these cards into the card read punch and read 
them into the program. These cards were checked and interpolated to give a 
little closer reference values . These solutions check nicely with our prev- 
ious work and actually showed some inaccuracies in the previous work. It 
turned out that it was rather tedious searching for cards , loading them in 
the card read punch, processing and inserting them "hack into the card boxes. 
The computer could calculate the problem faster than we could handle the 
card look-up system. 

By this time our 20K computer had become overloaded from two sources. 
The time required for compiling and loading programs from cards and also our 
problems were simply requiring more than 20K storage Fortranwise to solve. 
Some of our larger programs were written in SPS and even these exceeded our 
core capacity. Rather than get additional 20K on our existing computer, we 
decided to get a 1620, U0K with disk drive. Along with this equipment came 
the Clock, Monitor and Fortran IID. 

The first problem facing us in our induction log problem was to put 
our R and Z tables on disk, so we could retrieve it under program control. 
As you know, the Fortran IID is a formated input program that will not accept 
the free form output that we had from our Forcom routine. We therefore were 
faced with taking our original free form output cards and converting them to 
a format suitable for our computer. The AFIT program would allow us to use 
free form input and E FORMAT output. By putting a record mark in 20,000, we 
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usc;-d this program to convert the decks, consisting of cards from 0. to 60. for 
both R and Z values in steps of .01. This gave us a total of 12,000 cards 
which would require 60 cylinders to store in the form used by the Record, Fetch 
statements of IID. In storing data on a disk you need to store the data as a 
block in some section and move it down into working storage when you run the 
problem. This work would have required 120 cylinders which is in excess of the 
100 cylinders available on a disk. Looking at the method, we obtained the 
values of from 12. to 60. by hyperbolic curve fit. I decided to calculate the 
higher values directly instead of table look-up. When I ran the multi-coil prob- 
lem I found values as high as ^5 not uncommon but we rarely went over 60. 

The DFINE statement is the heart of the working storage assignment. You 
will recall that in the Manual, it says you can reserve up to a maximum of 99 
work cylinders. I decided to use 32 work cylinders for working area. I had de- 
fined the disk sector address of the work cylinders as sector zero. When I tried 
to load this control card, the computer gave an error 13 message, "insufficient 
available storage for specified work cylinder 1 '. You will remember the DIM table 
starts in Cylinder 2h. I decided the work cylinder area could not be greater 
than 2k. If it was greater, then it would require it to jump over the DIM and 
equivalence table. Whether this is true of not, I did not determine. I decided 
to use 20 working cylinders for data storage and this worked nicely. A short 
program was written using the Record statement which assigns the data card in 
sequence in the working cylinder area and assigned a record number to this data. 
By defining the work length 10 or less, it put a card per sector per record num- 
ber. Also, we were told by IBM that record number 1 starts in sector 219. I 
have yet to find this statement written "in the Manual as such. 219 from 3999 
leaves us 3780 sectors for data storage and this is the number of cards we used 
to load into data storage. If you note in the 1311 Manual, the time required to 
execute a Write Disk statement is two seconds. Since each card is handled indi- 
vidually, the 3780 cards took a little over two hours to load two boxes of data. 
Dollarwise, this takes $70.00 to load two boxes of cards onto disk. Kind of ex- 
pensive when you think of it this way. 

After loading the data cards into working cylinder areas on the disk, 
I wished to move them to permanent storage area on the disk so the program 
could recall them when running the problem. This is required because the 
Fortran compiler compiles in the first working cylinders and destroys the 
first few cylinders of tables. I also found that a couple of other working 
cylinders were used during this compilation.' I will get to this in a minute. 
You have two choices to move data from working cylinders up to permanent stor- 
age under Monitor control. One is by DL0AD and the other is by DC0PY. I used 
the DL0AD first to transfer this block of data up into permanent storage. I 
chose to move 3780 sectors by differences in sector address. The routine 
moved 780 sectors so I assumed that a maximum of 999 sectors could be moved by 
DL0AD control card. In using the DC0PY control card, you can successfully move 
3780 sectors and file protect. However, when you type out the Availability list 
you will note that it shows this area to be open and usable "even though it has 
been file protected and has data stored therein. I was careful to avoid this 
area, so as not to destroy the tables and have to reload. 

I had to rewrite the multi-coil program on the induction logging re- 
sponse equation in Fortran IID, I decided to make this a generalized program 
which would encompass U, 5> 6 and 8-coil designs. Since I did not want to get 
involved with the complexities of subscript notation, I defined all of my vari- 
ables directly, that is, RlTl f R1T2, etc. After compiling this program, it was 
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found the program required in excess of 50,000 core positions. As shown by the 
block diagram, the program was split up by making the 8-coil portion of the GB 
section a subprogram on call as this would be used the least. The data R. and Z 
tables were split up in two smaller groups of from to 18.99 for R tables and 
to l8.80 for Z tables. If the problem required values look-up greater then 
either of those stored in the R and Z tables, it used the hyperbolic equation 
constants previously derived to calculate these values. This master program 
was stored on Disk under the name MULTIC. 

The procedure to run the problem is to do a *DC0PY and move the data 
from permanent storage to working cylinders. Next, do a XEQS MULTIC which calls 
the program off of disk and moves it into core for execution. The input data 
cards are loaded as required by the program. , A 1620 Users Group program is 
available called GETREC which will do this DCOPY under Fortran Program control. 
At the time I wrote the MULTIC program, the GETREC program was not available. 

Several things you cannot do are these: You cannot do a *DC0PY after 
XEQS program because the DUP routine destroys the program.. You cannot do a 
FORX because in compilation, this evidently destroys an area of Cylinder 3 and 
Cylinder 10. I discovered this after I couldn't use the Z tables which were 
stored above Cylinder 9. If you watch the index wheel in the 1311 when com- 
piling, you will notice it goes to Cylinder 10 at least once. Sometimes the 
input data would be out of range and we wished to restart the program by a 
branch command on the typewriter. The program communication area starts in 
position 02218. The address beginning at 02222 is the restart address of the 
program. A branch using this address* or ^9022& will start the program if the 
check stop light is not on. 

Information available in this area starting with 02218 is Floating 
Point Word Length, Fixed Point Word Length, Starting Address of Main Line Pro- 
gram, Starting Address of Common Area, Number of Words in Logical Area, Number 
of Logical Records, Word Length and Record Length. 

The running time of this program for a 6-coil design was about k minutes. 
In generation of R tables, it took approximately 250 hours computer time and 
about 200 hours for the Z tables. A good average for the values that had to be 
looked-up would have required about 100 hours of integration time per value. For 
an 8-coil design problem you would have had 16 R values and 32 Z values to look- 
up. This would have been k8 values if integrated directly at 100 hours per 
value, which gives us about 1*800 hours for value determination. Since it only 
takes about k minutes to run this problem by using table look-up, we effectively 
have a time savings of about 72,000 to one. If I had had a 7090 or equivalent 
type computer, I might have gone ahead and done the proglem the long way, not 
aware of the time savings available by table look-up technique. We not only 
have a solution which is quite useable to make configuration studies, but we also 
have a set of tables which can be used without the computer by an engineer to 
solve this problem. 

Our' computer time runs in the order of $*i0.00/hr. and at approximately 
U800 hours solution the long way is some $200,000. One time through (using 
the Disk Storage) costs about $3.00. 

We have run many hundreds of solutions to optimize this design, so you 
see. the justification of this approach. 
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TH2 DIRECT AND CYCLIC JACOBI KETIIODS WITH FADE EVA * 3 
CORRECTION ALGORITHM FOR A REAL SYMMETRIC MATRIX 



One way of formulating the eigenvalue-eigenvector problem for a square- 
matrix A is to find a matrix T such that 

T- 1 AT = A (1) 



where A ::. r- •••ither a diagonal or "Jordan-Normal tf matrix. If A is diagonal * 
the dia^o,:. ... elements are the eigenvalues and the columns of T the eigen- 
vectors of A. If A is real-symmetric, then T is orthogonal which implies 

?~ 1 = T. In this case the existence of the diagonal form is assured, even 
for multiple eigenvalues • 

The Jacobi method finds the matrix T "by an iterative process. The 
diagonal form is realized by performing successive rotations on 2x2 sub- 
nat rices of A in such a way as to eventually anihalate all the off-diagonal 
elements (to within a specified tolerance). Two schemes for doing this are 
explained later. To understand the basic transformation the corresponding 
problem for a 2x2 matrix is considered. Let 

A - [ M2 ], r. ( *' ). (2) 

ai2 &22 X2 



Consider now the quadratic form associated with A, namely 

Q(x lf x 2 ) = xAx = a n x x 2 .+ 2a 12 x^ x 2 + a 22 x £ 2 (3) 



The curves Q(x^ , x 2 ) = const represent conic sections in the (x^ x 2 ) plane. 

From analytic geometry, we know that there exists a rotation of axes such 
that, in the new system, the x x x 2 term is absent. Hence, in the rotated 

system the associated matrix must be in the form 



* For arbitrary matrices, the diagonal form. is not always possible. 



The transformation which rotates the axes through an angle is, evidently, 

given by 



x * x 



/-coso smo\ 
" -sin9 cosO^ 



and the new quadratic form by 

Q(x x , x 2 ) = xAx - x(0A0)x 

Kence the matrix 

A = 0A0 

must be diagonal. The expression for 9 is well known: 

2 a l2 



tan 2G = 



ill - c 2 2 



The above principle constitutes the basis of the classical Jacobi 
method and its modifications. In all versions, the given matrix A is 
transformed to a new matrix A} defined by 

A 1 = TAT 

The matrix T- is made up of a 2x2 submatrix of the form. 0, with one ! s on the 
remaining diagonal position - and zero ! s elsewhere.^ the elements of C 
being so chosen that some pair of corresponding off-diagonal elements of A. 
will vanish. The matrix A 1 is again transformed in such a way as to create 

zeros in another pair of off-diagonal positions. In so doing, the elements 



originally zeroed by the first transformation will again attain n&#-zero 
values. However, it can be proven that, by repeated applications of such 
transformations , all of the off-diagonal elements can be reduced to 
arbitrarily small numbers . When this point has been reached , the eigen- 
values can be read from the diagonal elements, and the eigenvectors found 
by combining the successive transformations matrices into a single one. 

Two versions of the above mentioned have been programmed for the 1620. 
The first, based on the classical approach suggested by Jacobi, always 
anihiiates the off-diagonal element of greatest magnitude. The second 
version, sometimes called the cyclic method "anihalates" the off-diagonal 
elements in some systematic order, for example: 



skipping any which are less than a specified tolerance, and making as many 
passes as are necessary to reduce all to the required smallness. 

In general, no appreciable difference in time between the two versions 
could be detected, although in some cases the cyclic method was somewhat 
faster. 

Several refinements are included in the present programs. Equation 
(9) was replaced by the following two 




a 2 3 
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where 



ot 



(a ii- a jj )/2 



(13) 
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end (i , j ) the indeces of the particular elements to "be anihalated. These 

wore only usee as . long . as p was large enough that no loss of accuracy 
resulted from the subtraction under the radical. Thereafter cos 6 and sinO 
were computed from the following series 



sir .O 1 A. [i . I ( ^- ) 2 + |i ( 1- )3 ... 



(no 



Finally a device due to Fadeeva( [l j > p 1*8*0 was added which improves the 
final approximations to the eigenvalues and eigenvectors to within a higher 
tolerance than that attained by the method itself. The correction takes 
the form 



X. * a. 4 



with a similar correction to the components of the computed eigenvectors. 
The correction requires modification if there are multiple eigenvalues. 
Also , the algorithm is only valid if the off-diagonal elements are small 
compared to the lowest eigenvalue. When dealing with matrices such as the 
Kil'bert matrix in which there is wide separation in the greatest least 
eig ■ lue, a knowledge of the magnitude of the. lowest one is helpful. 

The method has been tested on segments of the Hilbert matrix. A check 
on the accuracy of the results, can be made by comparing the product of the 
eigenvalues with the value of the determinant , for which the exact expression 
is known [2]. A second test matrix due to W. Frank [3] was tried. In this 
examp~ e the exact eigenvalues were known which allowed a direct comparison 
to be made between the uncorrected and corrected values. In each case the 
results proved that the corrections were of the order of magnitude of the 
original tolerance squared. 
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EXHIBIT MAIN 



C MAIN PROG TO TEST VARIAbLi PRECISION JCLiI oUu 

DIMENSION A (2A, 12 ) ,B( 12,12) , EE ( 12) ,0< 12; 
COMMON A, 2, 22" 
1 ACCEPT 1C2,2R 

ACCEPT 1 1 » M 
i\ i l = N-.-l ■ 

NN-N+N 

I r ( SENSE SWITCH 1-) 200*203 
20 DO 201 I = 1 , N 
DO 201 U = 1 , I 
READ 102 >A ( I >J) 

201 A ( J j I ) = A ( I » J ) 
GO TO 202 

203 CALL LI ( D E T » N ) 
CALL M A T X ( N) 

202 CALL JCQI(A,N,ER) 

205 CALL IMPEV(N»EE»Q) 
S2=1.0 

51=1.0 

DO 206 K=1,N 
S1=S1*EE(K) 

206 S2=S2*A( K,K) 
S3=S1-DET 
S4=S2-DET 
PUNCH 120 

PUNCH 121,DET,S2»S1 
PUNCH 122 , ' S4 ,S3 
CALL VEC(A»B»N) 
PUNCH 112 

PUNCH 115, ( A{ I , I ) ,Q( I > ,EE( I ) , 1=1. N) 

' PUNCH 109 

DO 26 J=1,N 

DO 25 I=N1»NN 

K=I-N1+1 
2 5 PUNCH 110,A( I »J) ,B(K,U) 

26 PUNCH 106 

I F ( SENSE SWITCH2)6»10 

6 PUNCH 114 

114 FORMAT ( 15X25HINPUT MATRIX D I AGONAL I ZED/ ) 
DO 7 1 = 1, N 
DO 8 J=1,N 

8 PUNCH 113,A(I,J) 

7 PUNCH 106 
10 PAUSE 

I F ( SENSE SWITCH 3)9,1 

9 DQ 40 I = N1 , NN 
DO 40 J=1,N 
X-I-Nl+l 

IF( I-J-N)41» 42,41 
42 A(K,K)=EE(K) 

GO TO 40 
41 A( I » J ) =B ( K » J ) 

40 CONTINUE 

GO TO 205 /3 



EXHIBIT MAIN (CONTINUED) 



106 



FORMAT ( E 26.19 ) 
FORMAT (/)■ 



112 FORMAT (4X17HAPPROX-EIGENVALU£:'14X10riCORRECTIOM5X14riNcW-wI6cNVALi! 

10 9 FORMAT C / 1 5X7HVECTORS/ ) 

10 2 FORMAT ( E 14. 7) 

101 FORMAT (1.2) 

104 FORMAT ( / 15X17HEXAC7 EIGENVALUES/) 

110 rORMAT(2E30.19) 

115 FORMATS £26. 19 » 1XE26 .19 ♦ 1XE26» 19 ) 

120 FORMAT { 3X1 7HEXACT DETERMI WANT 10X12HBEFGRE CORR.15Xiir.AFTER.CORR, 

121 F0RMAT(E26.19*2(1XE26.19) ) 

122 FORMAT ( 26X * 2 ( IXE26 . 19 ) ) 



END 
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EXHIBIT A 

DIRECT JACOBI METHOD 

LARGEST OFF DIAGONAL ELEMENT METHOD 

A=AiI>J) 9 A REAL SYMMETRIC MATRIX DIMENSIONED (2N>N) 
C V:r\ER£ i\= ORDER Or MATRIX* £R=TOLL« 

C AFTER CONVERGENCE THE DIAGONAL OF THE UPPER NXN MATRIX 

C, WILL CONTAIN THE APPROXIMATE EIGENVALUES OF All 9 J) 

THE COLUMN OF THE LOWER HALF WILL CONTAIN THE APPROXIMATE EIGENVECTORS* 
SUBROUTINE JCBI ( A »N • ER ) 
DIMENSION A(24t 12) 
N 1 - N -v" 1 
N 2 - N — I 
NN-N+N 
27 DO 7 I-NI*NN 



■7 



J=l >N 



I F ( I —J—N ) 6 * 5 » 6 
5 A ( I > J ) = 1 • ' 

GO TO 7 
o A ( i ? J ! - « 

7 CONTINUE 
2 5^A£2»i) 

M=2 

) — "! 
L- — X 

I F £ N-2 ) 26.26.25 
2 5 DO 9 1=3 »N 
11=1-1 
DO 9 J=1,U 

IFC A; I »J)*A( I.J)-3*3)9»9.8 

8 S = A £ I » J ) 
M=I 

L=J , 

9 CONTINUE 

26 IF(B*B-ER*ER) 24*24*21 

21 AL=('A(L»L)-A(M»M) J/2.0 
I F £ AL ) 28x22.28 

2S X=B/(2.0*AL) 
Y = X*X 

I Ft'Y-.OOOl ) 23 .23 >22 

23 O1.0-.5*Y+1«375*Y*Y-4.3125*Y*Y*Y 

S = X* ( 5*Y+3.875*Y*Y-11«6875*Y*Y*Y) 

GO TO 12 

22 R2=AL*AL+B*B 
R=SQRTF ( R2 ) 
C2=.5*( 1.0+AL/R) 
S2=l*0-C2 
C=SQRTF(C2) 
IF£B) 10. 11.11 

10 S=-SQRTF£S2) 
GO TO 12 

11 S=SQRTFCS2) 

12 DO 13 1=1. NN 
D=C*A( I »L)+S*A( I »M) 
E=-S*A< I »L)+C*A< I »M) 
A( I »L> =D 

13 A(I.M)=E 

14 DO 15 J=1,M 

D=*C*A( L, J)+S*A(M»J) ■ /D 

E=-S*A£L,J)+C*A(M,J) 

A(L.J)=D 

15 A(M»J)=E 

40 I F( N-2 ) 24.24,20 

24 RI~URN 



CYCLIC JACOB I METHOD 

A=A< I»J) » A REAL SYMMETRIC MATRIX DIMENSIONED (2f<».N) 
VJHER E N= ORDER OF MATRIX? E.K = TOLL» 

AFTER CONVERGENCE TME DIAGONAL OF THE UPPER NXN MATRIX 

WILL CONTAIN THE APPROXIMATE EIGENVALUES OF A(I»Ji 

THE COLUMN OF THE LOWER HALF WILL CONTAIN THE APPROXIMA 

SUBROUTINE JCBI(A»N,ER) 
DIMENSION A (24, 12) 
Nl=N+l 
N2=N-1 
N N = N + N 
N.T=(.N#N2)/2 
DO 7 I«N1»NN 
DO 7 J=1»N 
IF( I-J-N) 6,5*6 
A( I » J ) = 1.0 
GO TO 7 
• i I » J ) = • 
CONTINUE 
K=G 

DO 17 L=1»N2 
L1=L+1 

DO 17 M=L1,N 
B = A ( L » M ) 

IF(B*B-ER*ER) 16,16,21 . 
AL=(A(L»L)-A(M»M) )/2.Q 
IF(AL) 28*. 26*28 
C2=.5 
GO TO 25 
28 X=B/(2.0*AL5 ' 
'Y = X*X 

IF( Y-.0001 }23 »23,22 

C=1.0-.5*Y+1.375*Y*Y-4.3125*Y*Y*Y . 

S=X*( l.-1.5*Y+3.8 75*Y*Y-ll .6875*Y*Y*Y) 

GO TO 12 
22 R2=AL*AL+B*B 
• R= SORT F < R2 ) 

C2=.5*( 1.0+lAL/R) } 
25 S2=l.-C2 

C=SQRTF(C2) 

IF(B) 10*11*11 

10 S=-SQRTF(S2) 
GO TO 12 

11 S=SQRTF(S2) 

12 DO 13 1=1, NN 
D=C*A( I »L)+S*A( I ,M) 
E = -S*A( I tL')+C*A( I »M) 
A ( I » L ) =D 

13 A ( I ,.M } = E 

14 DO 15 J=1,N 
D=C*A( L» J )+S*A(M, J) 
E=-S*A(L»J)+C*A(M,J) 
A(L,J)=D 
A ( M , J ) = E 
IF(N-2) 24,24,17 

16 K=K+1 // 

17 CONTINUE 
IF(K-NT) 19,24,24 

24 RETURN 
END . 



EXHIbIT C 



C F ADiiEVA » S EIGENVALUE CORRECTION ALGORITHM 

SUBROUTINE IMPEV ( N » EV ,RL$UD ) 
DIMENSION A(24»12) >EV( 12) -»RESUD( 12) 
COMMON A ■ 
DO 6 1=1 >N 
W = 0.0 

DO 5 J=1>N 
IF ( I -J) 7.5.7 
7 Y= A ( I ♦ I 5 -A ( J * J ) 

X=A( I ♦ J ) 
Z=X*X 

IF(Y)3»4>3 

4 PRINT 100 
GO TO 8 

3 W-W+Z/Y 

5 CONTINUE 
RESUD ( I ) =W 

6 E V ( I ) =A ( I » I )+W 
3 RETURN 

100 FORMAT { 56HMATRI X HAS MULTIPLE EI GEN VALUES -CORRECT I ON METHOD FAI l.3 • 
1/) 

END 



n 



EXHIBIT D 



FADEEVA » S EIGENVECTOR CORRECTION ALGORITHM 

SUBROUTINE VEC(A»B»N) 
DIMENSION A<24, 12) ,b< 12,12) ,EV<12) 
COMMON A*B»BV 
NN=N+N 
' N1=N+1 

DO 6 K = 1,N 
DO 4 M = 1»N 
IF(K-M)2»3»2 
EV(.v'!)=1.0 
GO TO 4 

EV(M)a-A(M»K)/(A(K»K)-A(M*M) ) 

CONTINUE 
LL=0 . ' 
DO 6 I=N1,NN 
LL=LL+1 
S=0.0 

DO 5 J=1»N 

S=S+A( I » J ) *EV ( J } 

B(LL,K)=S 

RETURN 

END 



If 



EXHIBIT E 

COMPUTES THE DETERMINANT OF A HILBERT MATRIX (EXACT) 

SUBROUTINE LI(DET.N) 
G=1.0 ' 
FAC=1.0 

C = N 

CN=C**N 

Ci = C*-*2 

N1 = N-1- 

DO 2 K=1»N1 

FK=.< 

FAC=FAC*FK. 
J = N-!< 

G = G* ( ( FAC**2 ) / ( Ci-FK**2 ) **J ) 

CONTINUE 
DET=G/CN 
RETURN 
END 



11 



iXnlLil'T. F 

C COMPUTES ELEMENTS OF HILBERT MATRIX 

SUbROUT I NE MATX ( N ) 
I MENS I ON A { 24- > 12 ) 
COMMON A 
DO 2 i = 1 s N 
X - i 

^0 2 J = I * I 
Y=J 

A( I *J) = 1#/(X + Y-1«0). 
2 A(J»I)=A(I*J) 

RETURN 
END 



TRANSMISSION LINE SAGS AND TESNIONS 



Perhaps the first person faced with the problem of tension and corresponding 
sags in a line stretched between two points may have been named ADAM* It is 
conceivable that Eve, having gotten Adam to sample the apple, caused an awareness 
of the need for fig leaves, and a consequent need for a fig line on which to hang 
the newly washed apparel . A grape vine being in sight, spanning some distance 
between two trees, may have offered the best solution Adam could master under 
the circumstances* Most likely, being a male, he would rather see the fig leaves 
on the wash line than on Eve - but be that as it may - a line strung between two 
rigid supports has both sag and tension and because of these characteristics, Adam 
may have had some trouble* Perhaps one day Eve hung one too many wet fig leaves 
on the grape vine and lo and behold several of the freshly laundred fig leaves 
began to drag in the dirt, making Eve very unhappy indeed* Ole Adam, seeing the 
predicament, figured he would fix that. So he quietly slipped around to one end 
of the grape vine and hitched it up a notch. Whoosh! The vine parted and Eve f s 
wash really in the dirt - you can imagine what Adam got for that! 

Ever since this fictitious event, we have been worrying about sags and correspond- 
ing tensions and vice versa* Through the intervening years, mathematicians and 
physicists have formulated the laws governing lines strung between supports* This 
turned out to be what is called a catenary and it is fairly simple to calculate 
the tension for a certain sag, or, if necessary, the sag for a given tension* 
Since temperature caused the line to expand as it got hotter and contract when 
it got colder, the sag would not stay put and neither would the tension* About 
this time some one found out that the line stretches some when you pull on it and 
if you pull too hard, it won f t return to its original length. And, of course, 
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if it is pulled harder still, it breaks, and you start over. Mathematicians and 
physicists are very smart people* They found a way to get the temperature and 
stretch into a formula for determining the sag and tension when the temperature 
and the stretch characteristics of the line are given. 

The line crews being people, and not caring to string a line in freezing weather 
with 1/2-inch of ice and with an 8-mile per hour wind blowing on it just to 
provide adequate ground clearance for a certain hi -voltage under these conditions, 
asked them (the physicists and mathematicians) to go back and figure some more. 
They preferred to string this wire on a nice sunny, 70° day at some tension that 
would guarantee that when that same wire was subjected to zero-degree weather, 
and had a certain amount of ice formed on it, the tension would not exceed a given 
value. Also, that the sag could be depended upon to still clear the ground by 
the desired amount. 

This caused the physicists and mathematicians no end of trouble. About the only 
way they could do this was to start with the specified tension at the loaded 
condition, determine the permanent stretch, the elastic stretch, the long time 
creep and then the corresponding sag. Working backward from the loaded condition 
to the condition of a nice sunny day to determine the tension to use, became a 
cut-and-try problem until the conditions of stringing on the warm sunny day would 
produce the desired result on the worst anticipated day. The cut-and-try method 
involved a lot of laborious calculations and good guessing to minimize the 
number of trial calculations . 

Many methods have been devised to simplify this chore. One of the methods was 
devised by Mr. Francis J. Hubert, of the Division of Water and Power, City of 

71 
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Los Angeles* The method of computation used in this computer program is based 
on his technical paper, "Simplified Sag-Tension Equations for Power Lines* 11 
The equations developed by Mr. Hubert are empirical equations that give results 
within 1% accuracy and are easily adapted to computer programming. See EXHIBIT 1. 

The equations used are based on the diagram shown here. From this diagram, it 
can be seen that L is the span length, the right-hand support is elevation -high, 
the support on the left elevation-low. Tension will be that in the conductor at 
the high support. The solution gives this tension, the vertical sag below the 
low support , and the distance from the low support to the lowest point of the sag. 
Of course, for supports of equal elevation, the low point is mid-point of the span. 

It was desirable in developing this program to have a direct and simple input 
format, a complete and readable output, covering the normal range of expected 
temperatures, and to provide automatic conformance to the National Electrical 
Safety Code Rules. Rule 1 of the safety code stipulates that the maximum tension 
shall not be more than 60% of ultimate strength. Rule 2 . the tension at 60°F, 
without external load, shall not exceed 25% of ultimate tension in the final 
unloaded condition. Rule 3 . the tension at 60°F, without external load shall 
not exceed 35% of ultimate tension in the initial unloaded condition. 

Most utilities are somewhat more conservative and, to provide more safety from 
Aeolian vibration, apply the last two rules, not at 60°F, but at the lowest 
temperature expected (usually zero degrees). 

The Program Flow Diagram, Exhibit 2, shows the method used to insure conformance 
to the safety rules previously mentioned. TMAX is an input quantity based on 
Rule 1. If TMAX is entered with a value larger than Rl percent of ultimate, the 

2 3 
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program automatically reduces it to this value* A calculation is then made to 
find the corresponding tension for the unloaded final condition at the same 
temperature* This value is compared with R2 percent of ultimate* If the calcu- 
lated tension is less than B2% of . ultimate, the program proceeds* If the calcu- 
lated tension exceeds R2 percent of ultimate, TMAX is reduced and recalculation 
takes place until Rule 2 is satisfied* 

Rule 3 percent of ultimate is then used to calculate the stringing tension at 
the same low temperature* This tension is then compared with the previously 
used TMAX. If the calculated value of tension T4 is less than that produced 
by TMAX, the program automatically makes TMAX = Rule 3 percent of ultimate with 
the assurance that if TMAX passed Rule 1 and Rule 2, and R3 times ultimate was 
the largest value that would pass Rule 3, the new value of TMAX would consequently 
pass all three rules* If T4 were greater than that produced by TMAX, then Rule 3 
would automatically be passed* With the rules having been cleared with an 
acceptable TMAX, the calculations proceed through the specified number of tempera- 
ture increments to produce the output results by following the Program Flow Dia- 
gram, Exhibit 3* 

Data input is accomplished by 5 input cards — 3 for identification, 1 for conductor 
data, and 1 for span data* A form was developed for punching data as illustrated 
in Exhibit k* 

This exhibit is a picture of 5 input cards with the data arrangement shown. Note 
that the first 3 cards are identification. Any wording punched will be reproduced 
on the output. Cards 4 and 5 are conductor and span data cards. They include: 
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Conductor data 


Span data 


Type 


No. 


Size 


Length 


Strand 


ELevation-Hi 


Diameter 


KL e va t i on-Lo 


Area of Conductor 


Wind pressure 


Area of steel 


Thickness of Ice 


Area of aluminum 


xni uicti bcinpcrdbUi 


Weight 


1 einperci uurt; _lih^ i tziucsii u 


Ultimate tension 


Final temperature 


Maximum tension allowed 


Rule 1 % 


Creep constant 


Rule 2 % 


Largest Span No. 


Rule 3 £ 




Loading 




Loading Code 



The program is stored on disk pack and when called into the computer, excecutes 
the program with the above data to produce the output as follows: 

Alternate No. 1 EXHIBIT 5 

- only final sags and tensions calculated 

Alternate No. 2 EXHIBIT 6 

- without including long term creep 
Alternate No. 3 EXHIBIT 7 

- Enter exact values of modulus of elasticity and coefficient 
of expansion - otherwise it is calculated internally. 

Switch 4 is not used in the program. 
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Alternate No* 4 EXHIBIT 8 

With all switches off: 

- Full output, final and stringing tensions and sags, 
with 10-year creep, using calculated values of modulus 
of elasticity and coefficient of expansion. 

Calculation of each span with full output takes about 3^ minutes with the 1443 
on-line printer. 

Accuracy has been proven by comparison with other calculations. One notable 
comparison was on a proposed long span of 10, 000 f . 

The results compare very closely with the results obtained from the Copperweld 
Company for the same span and conditions. EXHIBIT NO. 9* 



TMAX - 


RLxTu 


Final, 


Loaded 


i 


STC 


TA 


Pinal, 


Unload" 




ed 



TAS R2xTu 



TA>R2xTu 



T » R3 
Initial 


x Tu 
., Unload. 




STC 


Final, Loaded 




T4<TMAX 



Reduce TMAX 



TMAX - T4 



Use last TMAX and STC to compute 
sag and tension at different 
temperatures for final and 
stringing conditions 



Fig» 3 Abbreviated Plow Diagram for 
Sag-Tension Computation 



(MAXIMUM TENSION) 



Eq. I 



TRIAL 
SAG 



LL 



v/w 2 + w d 



8T-LI55L Z (W z + W*)-*Wh 



Eq.2, Iteration 



NON-LEVEL SPAN 



I * -f hs l" 

LL„ 



8T-4Wh-8Vw d Z +1-33 (W^ 



At 



ijf Eq. 3, Iteration 



LEVEL SPAN 



100 + 



Z, ^ $ Eq.! 



WW 



iJkr- 100 SAT 



V/W) 2 + '(W d L c ) 2 

AE 




% cr = ^5625 + 4 (%S) Z -75-_Z_ H 

800 |_ %S + 1450 _ 



* cr, = $ cr + Z$£or 



Eq. 10 



#s 2 

LEVEL SPAN 



1 



800 



Eq.3 



NON-LEVEL SPAN 



is. 2 

37.5 ^ 



too + 



h i s 



TTTV 



Eq.2 



3.76 (| + 

LW - 



WL 



2 r 



83. 



,5hS 



+ S VW . 2 + 



L fi L v LL C J 



FIGURE 2 SIMPLIFIED FLOW DIAGRAM OF STC 



WL ' 

1.33 (—) + Wh 



EXHIBIT 3 



IBM DATA FORMAT SHEET- SAG AND TENSION PROGRAM 



SHEET 



/ 



.OF 



/ 



10, 



15, 



20, 



25, 



30i 



35, 



40, 



45, 



50, 



55, 



60, 



65, 



70, 



75, 



80 



:ard 1 
:ard 2 
:ard 3 



JUL 



, ,ST,EAR,^-,R,0,GEP, ,C,ORP,., 



I I I I I 1 I I I 1 M i I I 1 I I I I 



-L 1.1 



.S.Afi .AMD, ,T,E i N t S | 1 ,0,N, .PftO.GftAft 



till 



1 1 1 ' 1 1 1 1 1 1 { i 



i I 1 I 



I 1 I 



1 1 I I 



CONDUCTOR DATA (1-CARD) 



x 

TYPE 



SIZE 



STRAND 



DIAMETER 



AREAC 



AREAS 



ARE AL 



WEIGHT 



ULT. TENSION 



MAX. TENSION 



CRPK 



SPN 



CARD 4 



CARD 5 



4 



10. 



15) 



20 



25. 



30 



35i 



40 



45 



50* 



55i 



60| 



651 



70, 



30 



1 1 1* 1 -^"° 1 1 A*e 



lure i /aA*& 



NO. 



LENGTH 



EL- HI 



EL-LO 



W. PRESS. 



ICE THK. 



I. TEMP. 



INC. 



F. TEMP, 



RULE \ 



RULE2 



RULE3 



LOAD 



o 



X 

u 



3 



CONDITIONS AND SPAN DATA (1-CARD) 



in 
a 
cr 
< 



* 

+ 

T 

X 



LOADING CODE HEAVY - t. MEDIUM '= 2, LIGHT =3 

CREEP CONSTANT STEEL AND COPPER - .5, ACSR = 1.0. ALL ALUMINUM 
AREAS CANNOT BE ZERO - IF NONE. ENTER ,0001 

FIXED POINT (WHOLE NUMuEkS} 
ALPHABETIC CHARACTERS 

ALL CHARACTERS TO BE RT. HAND JUSTIFIED 



~ 1.5 



EXHIBIT 4 



FORM El 3. 1 



o 



o 



STi-ARNS-KGM.'R COKr. SAG AND TENS inn PROG RAH E]3 

jou no. ui-i?j)2 custumer common noting proj. -xi-nrm r > 

SObJlGl SANR!!: C A I. G l J !.. A T 1 1"; N - MEDIUM SPAfi RY J.k. SMITH OATf: '(/'!/(•■ 



AM LENGTH EL-HI EL-LG COND SIZE STRAND OF I GOT AREA D I A . LOAD TEN;, 
1 1200.0 13!?. 125.0 AW 248NGM 19/9 .5658 .1954 .572 HEAVY 1314 



I b;'i PfckA 1 URE 




FINAL 




DEGREES 


TENSION 


SAG 


LOW PT. 


0.0 


13171.5 


21. 528 


542. 9 


0.0 


8570.6 


7.423 


473.9 


10.0 


8330.9 


7.752 


477.4 


20.0 


8095.8 


8.094 


480.9 


30.0 


7865.7 


8.449 


484.3 


40.0 


7640.7 


8.818 . 


487.6 


50.0 


7421.2 


9.200 


490.8 


60.0 


7207.4 


9.596 


494.0 


70.0 


6999.5 


10.005 


497.0 


80.0 


6797.6 


10.427 


500.0 


90.0 


6601.9 


10.862 


502.9 


100.0 


6412.6 


11.308 


505.7 


110.0 


6229.7 


11.766 


508.4 


120.0 


6053.4 


12.23 5 


511.0 



TENS ION 



STRINGING 

SAG LOW PT 



c 



c 



3/ 



EXHIBIT 5 



T /' i ; C — "• r. c rz p 


CORP. 






SAG AND TENSION PROG R A 


;.; f 1 3 


JOB MO. OL-128 


2 CUSTOMER 


COMMON 


MEETING 


PRO J. 


t~ 
L 


XH1BIT 6 




SUB 1FCT SAf'i D LE 


CALCULATION - MEDIUM SPAN 


BY J.R. SMITH 


DATE. 


7 / V / 6 6^ 




EL-LU C0MD SIZE 


STRAND 


HEIGHT AREA 


DIA. LOAD 


T Hi-! ST f 


1- 1200.0 135.0 


12 5.0 AW 


2 48MC M 


19/9 


.5658 .1954 




572 HEAVY 


13145 


1 L lire r\ A I U r\ C 




FT NAL 






s 


T RINGING 




U c b is c tl o 




O H vj 


LOW PT. 


TENSION 




SAG LOW PT. 





1 3145*7 


21 . 579 


543.0 










n o 


8 570-6 


7.423 


473.9 


10118.5 




5.695 


451.1 


inn 


W _J V # J 


7. 752 


477.4 


9855.6 




5.948 


4 54.9 


? - 


8095.8 


8.094 


480.9 


9595 .9 




6.213 


458.3 


^ n n 

J? U • u 


7865.7 


8 • 449 


484.3 


9339.6 




6.489 


46 2.5 


&n n 


7640.7 


8.818 


407.6 


9087.0 




6,778 


466.3 




7421.2 


9.200 


490. 8. 


8838.3 




7.079 


469.9 


a n n 


7207.4 


9.596 


494.0 


8593.7 




7.393 


473.5 




6999.5 


10.005 


497.0 


83 53.5 




7.720 


477.1 


8 0.0 


6797.6 


10.427 


500.0 


8118.0 




8.061 


480.6 


90.0 


6601.9 


10.862 


502.9 


7887.3 




8,415 


484.0 


100.0 


6412.6 ■ 


11.308 


505.7 


7661.9 




8.782 


487.3 


110.0 


6229.7 


1 1 . 766 


508,4 


7441.9 




9.163 


49.®.. 5. 


120.0 


6053.4 


12.2 35 


511.0 


7227.5 




9.556 


49 3.7 



UPLIFT AT LOWER SUPPORT OCCURS WHEN BOTH SAG AND LOW PT ARE ZERO 



o 



EXHIBIT 6 



S "i E A P !\ ,:, -KG GE R CORP. SAG AND TENSION PROGRAM El 3 

.'''Oi : OL 12><2 COSTGMER COMMON MEETING PROJ. EXHIBIT 7 

SUb/LCO SAMPLE CALCULATION - MEDIUM SPAN . BY J.R. SMITH DATE- 7/7/,' 

... ; : L L-H I EL- L COMD SIZE STRAND WEIGHT AREA 01 A. LOAD TEN', 
TN ! ( . 1.35.0 125.0 AW 248MCM 19/9 .5658 . 1934 .572 HEAVY 13i 



i EN PE 



FINAL 





T C M C T n M 

1 b N i> i U N 




i nil r\ t" 

LOW PT« 


T ENS ION 








CA 5 O 

H c. 9 O 




C * 


8570.6 


7.42 3 


473.9 


10127.9 


1 ,0 


8319.6 


7.768 


477.6 


9852.4 


2 • 


8073.7 


8.127 


481.2 


9580.4 


J ® 


7833.4 


8.501 


484. 8 


9312.2 


r . 1 * N.-' 


7598.8 


8.889 


488.2 


9048.1 


5(5. C 


7.370.2 


9.293 


491.6 


8788.3 
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1. INTRODUCTION 



1.1 Numerical Control (N/C) is the operation of a device under control of numerical data 
(usually recorded in punched tape). It is the process of controlling machinery by auto- 
matic means rather than by manual means. An example would be the controlling of a 
drill press with the data recorded in a punched tape. 

1 .2 AUTOSPOT (AUTOmatic System for Positioning Tools) III is a general purpose computer 
program, developed by IBM for use with the 1620 Monitor I System, designed to aid in 
preparing data primarily for point-to-point machine tools. The program also provides 
limited contouring capabilities. 

1.3 This report presents a general approach to organized postprocessor development for the 
AUTOSPOT III program. The approach is centered around a structural outline with an 
emphasison modularization and standardization. This outline is demonstrated in the form 
of a postprocessor input module (AUTOIT) thatallows for the expansion of the AUTOSPOT 
III input language and a postprocessor program for the 3 BHT-B Burgmaster Turret Drilling 
Machine. 

2. AUTOSPOT III POSTPROCESSOR OUTLINE 

2. 1 Postprocessor Modularization 

The functionsof a postprocessor are divided into five major modules or sections. Similarly, 
each major section can be divided into various sub-sections. Modularization can best be 
attained using the FORTRAN ll-D language; therefore, whenever possible, this language 
will be used in coding the various sections of the postprocessor. 

2. 1 . 1 Control Section. The Control Section initiates execution of all five major sec- 
tions and maintains control of the flow of data in the postprocessor. This section 
should be very near the same for all postprocessors. 

2. 1.2 Input Section. The Input Section reads input for the postprocessor and transfers 
it into a buffer area in a prescribed format. This section can be used to expand 
the AUTOSPOT language through the use of the POSTPR/ statement. (See 
Section 3.) 

2. 1.3 Geometry Section. The Geometry Section computes all of the axes movements 
and converts this data to the machine tool coordinate system for the Output Sec- 
tion. Any adjustments to the coordinate data will also be provided in this sec- 
tion. An example would be the adjustments required for coordinates that are 
outside the machine tool limits. 

2. 1.4 Auxiliary Section. The Auxiliary Section executes action specified by auxiliary 
instructions and develops auxiliary command blocks for the Output Section. An 
example would be a value specified for scaling a plot or an auxiliary command 
generated to set a particular machining mode, The tables used by this section 
should be standardized as much as possible. 
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2. 1.5 Output Section . The Output Section provides for editing of all output informa- 
tion. The outputdata will be printed and plotted asitis received from the vari- 
ous other sections. The information to be punched in the control tape will be 
stored on disk until processing of all of the AUTOSPOT III output is completed. 
If no unrecoverable errors have been detected, the information will then be 
punched in the control tape. 

2 f 2 Standard Common 

In order to provide interchange of data between the sections, certain internal data (varia- 
bles, parameters, switches, etc.) will be part of common storage. Common will be divided 
into areas for each major section of the postprocessor. The organization of the common 
area along with the standard sections will simplify communications between AUTOSPOT 
III postprocessor programmers. 

2.3 System Description 

2.3. 1 The Control Section initalizes the postprocessor then calls the Input Section to 
provide information in the input buffer area. Control will identify the record 
type and determine the section required to initiate processing of the information. 

2.3.2 The Geometry Section is called to process coordinate information. A discrete 
block of associated data (machine tool coordinate, spindle specification, auxil- 
iary function, etc.) is generated for each coordinate record input. Once a 
block of data has been generated, its relationship to neighboring blocks of data 
is determined. (An example would be if there is a minimum or maximum departure 
required from the previous set of coordinates.) If the relationship is acceptable, 
the data is printed, plotted, and stored on disk to be punched in tape following 
processing of the entire AUTOSPOT III output. If the relationship is unaccept- 
able, suitable correction is obtained if possible, and the proper diagnostics 
printed. If suitable correction is impossible processing will continue but the 
control tape will be eliminated. 

2.3.3 The Auxiliary Section will process the non-motion instructions in the input buffer <, 
Certain auxiliary instructions will simply set program mode and parameters, others 
will execute special logic, and still others will develop auxiliary command 
blocks. Thus the Auxi liary Section wi 1 1 be a set of cal Is to sub-sections designed 
to react for a specific instruction. 

2.3.4 The Output Section will be called by the Geometry and Auxiliary Sections to 
plot and print the output information and to store on disk the information that is 
to be punched in the control tape. When processing of the AUTOSPOT Ell out- 
put is complete, Control will call Output to punch the information stored on 
disk in a control tape, if no unrecoverable errors have occurred. 

4 
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3. AUTOIT (AUTOspot Input Translator) 



3. 1 General 

AUTOJT is an input subroutine for AUTOSPOT II! Postprocessors, The program reads 
AUTOSPOT IH output and translates it into a prescribed form. Also, the subroutine pro- 
videsa method for increasing the AUTOSPOT III language through the use of the POSTPR/ 
statement. The program is written in SPS ll-D and can be used by both SPS !I-Dand 
FORTRAN ll-D programs. 

3.2 POSTPR/ Statement 

AUTOIT will accept a maximum of 79 characters following the slash in the POSTPR/ 
statement. A combination of postprocessor words, fixed point numbers, floating point 
numbers, and special characters can follow the slash. The postprocessor words must be 
defined in the table provided in AUTOIT. The fixed point numbers can be a maximum of 
4 digits and the floating point numbers can be a maximum of 8 digits. The special char- 
acters "(,*$) can be used to separate the above three type of fields. 

3.3 Subroutine Description 

The output from the AUTOSPOT III program consists of variable length records wri tten on 
disk. The contents of these records are processed as follows: 

(a) Fixed Point Numbers. The fixed point numbers are from two digits to six digits in 
length. The fixed point numbers in each record are converted to four digit fixed 
point numbers and stored sequentially in an array provided by the users. If the 
number is larger than four digits, it is separated into two fixed point numbers, the 
four lower order digits stored as one number and the remaining digits as a separate 
number. 

(b) Floating Point Numbers. The floating point numbers have either eight digit or thir- 
teen digit mantissas. All of the floating point numbers are converted to eight digit 
mantissa floating point numbersand stored sequentially in an array provided by the 
user. 

(c) Alphameric Data. Alphameric data is separated into four digit fixed point fields 
and stored sequentially in the same array used by the fixed point numbers. 

(d) Postprocessor Records. The floating point and fixed point numbers are stored as 
described above. Each postprocessor word is decoded according to the table pro- 
vided in the subroutine. The number code corresponding to each postprocessor word 
is stored sequentially in an array provided by the user. If the postprocessor word is 
PARTNO, thealphameric information following the slash in the POSTPR/ statement 
is stored in an array provided by the user. 
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3.4 Linkage 

3.4.1 FORTRAN ll-D 



CALL AUTOIT (IDSKRD, IFIX(l), FLOATf!}, IPP(l), IPTNO(l)) 

(a) !DSKRD = Initialize AUTOIT to read first AUTOSPOT record. 

= 1 Continue processing the next AUTOSPOT record. 

(b) ! FIX (1) The first field of the fixed point array. Minimum array 

size is 10. 

(c) FLOAT (1) The first field of the floating point array. MinEmum array 

size is 10. 

(d) IPP (1) The first field of the postprocessor codes array. Minimum 

array size is 10. 

(e) IPTNO (1) The first fieldof the array provided for the PARTNO state- 

ment. The minimum array size is 38. 



3.4.2 SPS li-D 



(a) The following statement must be entered in the user's program. 

(1) CALL LOAD, AUTOIT, @@SUB 

This statement loads AUTOIT into memory from disk and must precede 
the calling sequence. 

(2) @@SUBDSSSIZE 

@@SUB is the location where AUTOIT is to be loaded and must be 
defined in an even numbered location. SIZE is the length of the 
subroutine. It is variable because of the additions that may be made 
to the postprocessor word table. 

(b) The following sequence of instructions are used to "call" AUTOIT: 
BTM @@SUB, 11 

DSAIDSKRD, IFIX1, FLOAT1, IPP1, IPTNOl 

where the above fields are defined as follows; 

IDSKRD DC 4, 

IFIX1 DSB 10, 4 

I FLO ATI DSB 10, 6 

IPP1 DSB 4,10 

IPTNOl DSB 4,38 
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Each entry to the postprocessor word table requires six alphabetic and four numeric char- 
acters. Additional postprocessor words can be added to AUTOIT. These additions must 
be entered in the postprocessor word table following the remark: 

* ENTRIES TO THE POSTPROCESSOR WORD TABLE SHOULD BEGINE HERE 

A postprocessor word of six characters such as MACHIN with a number code such as 0002 
would be entered as follows: 

DAC 6, MACHIN 
DC 4, 0002 

A postprocessor word of less than six characters such as RAIL with a number code such as 
0010 would be entered as follows: 

DAC 2, XX 
DAC 4, RAIL 
DC 4, 0010 

4, AUTOSPOT III POSTPROCESSOR EXAMPLE 

The following is an example of an AUTOSPOT III Postprocessor that was written using the AUTO- 
SPOT Postprocessor Organization Outline described previously. The postprocessor is written for 
the 3 BHT-B Burgmaster Turret Drill. 

4. 1 Postprocessor Description 

4.1.1 AUTOSPOT Part Program Requirements 

(a) Specification Statements 

All of the specification statements can be used except the following: 

(1) BKC 

(2) SAFE 

(b) Major, Minor, and Geometry Statements 

All of the point-to-point major, minor, and geometry statements can be 
used. The following POSTPR/ statements can be used: 

(1) POSTPR/PARTNO, NNNNNN 

This statement specifies the part number 

(2) POSTPR/MACHIN, BURG 

This statement specifies the Burgmaster Postprocessor to be used. 
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(3) POSTPR/PPLOT, START, XX. X, YY.Y, SCALE, SS. S 

This statement initiates plotting at the point (XX. X y YY. Y) at a 
scale of SS. S. If this statement is omitted from the part program, no 
plot will be drawn. 

(4) POSTPR/PRESET 

This statement must be entered in the part program if the computed 
Z and W values ore to be output. If this statement is omitted from 
the part program, a Z of 2 and W of 1 will be output for each indi- 
vidual pattern generated. 

(5) POSTPR/RAIL, XX. X 

This statement specifies the rail height XX. X. If this statement is 
omitted from the part program, a rail height of 39 inches is assumed. 

4. 1.2 Spindle Selection. The spindle (1-8) will be selected by the SP(N) statement. 
The following statement would cause spindle 6 to be selected. 

TAP,20,SP(6)... o .... 

4. 1.3 Mode Selection. The machining operations tap and bore generate the mode E4. 
All other machining operations generate the mode E3. 

4. 1 . 4 Diagnostics 

Error Number Cause 

1 Thiserror occurs when the machine tool axis limits 

(X00000-600000, Y-00000-480000) are exceeded. 
The point in question is plotted and listed but is 
eliminated from the control tape. 

2 This error occurs when the minimum movement 

along the X or Y axis is not satisfied. That is, 
movement along both the X and Yaxes must be 
zero or greater than or equal to . 100. If this 
error occurs, the point in question will be stored 
and then processed when it is possible to do so 
without error. If it is not processed before the 
next spindle change, a "dummy" (X> Y) move- 
ment of .200 will be inserted followed by the 
point in question. 

3 Illegal input code. 

4 Illegal postprocessor word* 
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Illegal spindle number. 

Storage table for minimum movement points is 
exceeded. Maximum is 10. 



7 Illegal machining operation. 

The error messages are output on the printer and processing continues. 

4.2 Input Requirements. The input to the program is the output from the AUTOSPOT III pro- 
gram^ The input is read from disk using the AUTO IT subroutine. 

4. 3 Output Requirements 

4.3.1 Listing. See Figure 1. 

4.3.2 PLOT. See Figure 2. 



4,4 Control Tape 
Output Format 

1 . XnnnnnYnnnnnDn 4= 

2. XnnnnnYnnnnn =+= 

3. ZnnnnnWnnnnnEn 4= 

4. ZnnnnnWnnnnn #= 

5. ZW=N 

6. M#= 



Remarks 

This format is used when an X, Y coordinate and 
d new turret number are to be output. 

This format is used to output X, Y coordinates 
without a turret number. 

This format is used when a Z and W value is to be 

outputanda changein the mode of operation has 
occurred. 

This format is to be used when a Z and W value 
is to be output without a change in the mode of 
operation. 

This format is used to output a Z and W value that 
is identical to the previous Z and W values output 
and there is no change in the mode of operation. 

Stop command. 



Ten channel five punches will precede each output block containing a turret change (for- 
mat 1). Format 3 will be used any time a mode command is output. If a part number is 
input, it is imaged in the control tape preceding and following output data. 
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Figure 1. SAMPLE LISTING 
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Figure 2. SAMPLE PLOT 
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I. INTRODUCTION 



The increasing need for potable and ultrapure water, created by ur- 
banization, industrialization and population growth has caused local 
and regional shortages to arise. The extensive planning and forecast- 
ing carried on by several national governmental agencies and inter- 
national organizations has indicated that shortages would become limiting 
to at least economic growth of extensive geographic areas in the period 
beginning not later than 1980. One source of water, which matches the 
population concentrations in the coastal areas of many nations, is the 
sea. Desalting of this water for human and industrial purposes is not 
entirely new, but what is new is the attempt to do this on a large and 
relatively economical scale. 

The development of storage for and transportation of large quantities 
of imported water requires decades of planning and construction, as 
well as large capital investments. Therefore the U. S. Government 
has undertaken a program of desalting technology development and de- 
monstration to circumvent some of the costs and the time lag. This 
program aims to (1) decrease thfc cost of desalting water; (2) develop 
the experience and technology for many sizes of plants and processes; 
and (3) provide new water sources close to the points of need in a rela- 
tively short period of time. The program is under the direction of the 
Office of Saline Water, U. S. Department of the Interior. 

Stearns-Roger Corporation was named Management and Operations Con- 
tractor in April 1961, for the first production size Sea Water Desalting 
Demonstration Plant built by the Office of Saline Water, at Freeport, 
Texas. For the past five years, the purposes of plant operations have 
been, the continuous improvements in performance, and the development 
of reliable design data and hardware. About two and one half years ago, 
the operations reached the stage where optimization became important, 
and enough data was available to consider computer design programs. 
Stearns-Roger initiated advanced hardware design efforts which have re- 
sulted in a presently active, expansion program at the Freeport Plant. 
Several important process problems were solved which resulted in 
superior performance of the ME- LTV process, creating a high demand 
for proposals of completely designed plants to compete with the pre- 
viously popular Multi-Stage Flash process. This situation has resulted 
in the Computer Programs presented in this paper. 

II. THEORY AND PROCESS 

The Freeport Plant is an initial design of a Multiple Effect- Falling Film 
Sea Water Evaporator. It is a several- fold expansion of Multiple Effect 
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chemical solution evaporators. The process has a successful industrial 
history extending over more than fifty years. The plant represents an 
extension of twelve effects and has achieved a major improvement in 
thermal economy. It also, has successfully demonstrated scale free 
operation using the highly effective falling film evaporation in long ver- 
tical tubes. The process name is derived from these characteristics, 
since ME -LTV is the representation commonly used for Multiple Effect- 
Long Tube Vertical sea water evaporators. 

To fully identify the most useful sea water distillation form of the pro- 
cess, it is necessary to specify forward series feed, with regenerative 
feed heating. Many modes of operation have been investigated at the 
Freeport Plant, resulting in much simplification of the process and 
equipment. The advanced design equipment incorporating the simplifi- 
cations, is now being installed at the Freeport Plant and is represented 
by Figure 1, Process Flow Diagram -- ME Falling Film Plant. This 
diagram is also self-explanatory of the forward, series fed, regenera- 
tively heated feed, multiple effect evaporator. The diagram is also the 
process basis for the computer programs which follow. Symbols appear- 
ing on the diagram are listed in the appendix. 

Prior to the establishment of the preliminary computer program, it was 
necessary to make detailed process calculations on an effect by effect 
basis, starting with Effect No. 1. These calculations were necessarily 
based on a number of assumed flows and conditions. It required approxi- 
mately 24 man hours to complete a set using a desk calculator. If the 
results of the assumptions were not acceptable, it was necessary to ad- 
just a flow rate, temperature, etc. , and recalculate the entire set, just 
to bring the extraction ratio, duty, or some other result into an acceptable 
range. Arithmetic transposition errors so frequent in repetitive calcu- 
lations have the effect of invalidating all successive calculations in a set, 
since the performance of Effect n + 1 is totally dependent on the perfor- 
mance calculated for Effect n. This method of design calculations was a 
severe limitation on the consideration of alternates for any set of proposal 
conditions, and thus precluded most optimization other than the judgment 
of an experienced engineer. Geometry calculations had to follow a final 
set of acceptable process calculations. 

In the ME-*LTV system, the most important considerations which are un- 
der the designer's control are, to divide the total temperature spread be- 
tween heating steam and heat sink in some manner, and to fix the extrac- 
tion ratio. These two basic considerations must result in a plant which 
operates within boundary values fixed by the site, such as the production 
rate, sea or brackish water supply, available land area, and available 
energy sources. A more concise statement of the independent design 
variables is given in Table 1. 
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Most of these process variables can be changed arbitrarily by the de- 
signer by varying such things as the temperature profiles, the number 
of effects, the approach of stream temperatures in the heat exchangers, 
etc. , within the boundary limits previously mentioned. The costs of 
fuel and money are generally fixed by the user and are related to the 
site. There are a greater number of dependent variables which change 
each time one of the independent variables is changed. This illustrates 
the complexity of the design problem whose ultimate goal is the minimum 
cost of water, for most situations. Other desired end results are pos- 
sible for less common situations, and are listed in Table 1. 

Specifically, it was desired that a computer program be developed which 
would take a given set of terminal conditions, capital and fuel costs and 
determine the optimum plant characteristics in accordance with Table 2. 

The lack of adequate design correlations and reliable cost data still ham- 
per all phases of the effort to develop a complete program. These two 
problems are receiving much attention at present, with particular em- 
phasis on cost records for the module now being installed at Freeport. 

III. THE PRELIMINARY PROGRAM 

The development of simplified equations, from thermodynamic consider- 
ations, was followed by a check of their accuracy. This check was made 
by comparing the predicted performance to actual performance of the 
Freeport Plant, These equations and assumptions are contained in a 
previous paper by the author (Reference 1). They form the basis of the 
preliminary program depicted by Figure 2. This program requires con- 
siderable input by the designer, but does eliminate the problem of arith- 
metic and transposition errors. It runs on the 1620 IBM computer in 
slightly less than two minutes. 

To prepare the initial set of input data requires approximately two hours 
of an experienced engineer's time. The required input information is 
listed in Table 3, Additional input, where variations to be investigated 
are set forth, can be established in less than a half hour. Manipulations, 
based on results of the initial sets of input data, are made rapidly by re- 
punching only a few cards in a deck. Thus, it is possible consider eight 
to twelve alternatives in a period of one day; if the computer is-available. 

The output of this program is limited to the most .critical portions of the 
process and geometry portions of the desired end results. It has only 
one process loop, which adjusts the preheating temperature profile to 
increase or decrease the thermal performance around a target value, 
using the input number of effects. Even so, it provides the engineer 
with the necessary end results to make a logical judgment optimization 
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1. The Optimum Process Parameters 

a. N, the number of effects and operating temperature 
profiles. 

b. The duties, temperature differences and heat transfer 
coefficients of all heat exchangers. 

c. The velocity, quantity, quality and location of all pro- 
cess stream flows, for the required production rate. 

d. The fuel, electricity and chemical requirements for the 
required production rate. 

e. The Net thermal efficiency; extraction ratio, feed water, 
heating steam, product, and waste brine quantities for 
the optimum plant process. 

2. The Geometry of the Optimum Plant 

a. Heat exchanger dimensions, 

b. Vessel dimensions. 

c. Pump requirements; flow, head, driver horsepower. 

d. Vacuum equipment requirements; volume, pressures, 
cooling water and energy requirements. 

3. Cost Data for the Optimum Plant 

a. Capital cost of the complete plant to battery limits. 

b. Actual product water cost, with the total broken down 
into operating, maintenance and fixed costs. 

4. Sensitivity of the Selected Optimum Design 

a. To variations in fixed charge rates. 

b. To variations in energy costs. 



TABLE 2. DESIRED OUTPUT FROM A COMPUTER PROGRAM 
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READ TITLE CARD 
READ HFCP TABLE 

READ CONTROL CARD 




PASS 3 



PASS 1 

* 2 Ol^LY 



PASS 2 



PASS 2. ADJUST TSFI £ N'2 
PASS 3. ADJUST TSFI @ N 



READ N * I EFFECT CARDS 



PRINT INPUT DATA FOR 
IDENTIFICATION 
AND VERIFICA- 
TION 





CALCULATE FINAL PRODUCT 
TEMPERATURE AFTER PRODUCT 
COOLER 



CALCULATE ACTUAL 
VAPOR. SWF AND BRINE 
TEMPERATURES 



n 



FRI"Ni MESSAGE 
KX200 APPROACH 
IS UNRE ASON- 
ABLE 




CONVERT UNIT SURFACE AREAS 
TO TOTAL StrR FACE ARE AS 
FOR EACH EXCHANGER. CONDENSER, 
AND EVAPORATOR 



PRINT MESSAGE 
HX200 APPROACH 
IS REASON- 
ABLE 




SUM PREHEATER AREAS. INCLUDE 
HX200 

SUM EVAPORATOR AREAS. INCLUDE 
HXRJ 



CALCULATE UNIT DUTIES FOR 
EVAPORATORS AND FINAL 
CONDENSER 



CALCULATE TOTAL SEA WATER 
FEED RATE TO EVAPORATOR 



SELECT VAPOR ENTHALPIES 
CORRESPOND fNG TO TVO, RTGP fl) 



CALCULATE MEAN TEMPERATURE 
DIFFERENCES FOR HEAT 
TRANSFER, -- PREHEATERS, 
PRODUCT COOLER, FINAL CON- 
DENSER AND EVAPORATORS 



CALCULATE PERFORMANCE 
FOR EACH EFFECT (I) , 
WP. WPS, er. WCP. WVP 



CALCULATE UNIT HEAT TRANSFER 
SURFACE REQUIRED FOR EACH 
EXCHANGER. CONDENSER. AND 
EVAPORATOR 



CALCULATE OVERALL 
PROCESS PERFORMANCE 
ERN 



(5 




RETURN TO 
START 















LEGEND 








LEGEND 




N 


Number of Effects for this pass 


TVI 


Temperature of Vapor Into 


er 


Individual Effect Economy Ratio 


HXRJ 


Heat Rejection Condenser (Final 


PR 


Production Rate 




Effect I 


ER 


Assumed or Economy Ratio 




Condenser) 


RST 


Raw Sea Water Temperature 


TVO 


Temperature of Vapor Out of 


ERN 


Calculated Net Economy Ratio 


I 


The Effect under consideration 


SWF 


Sea Water Feed Rate 




Effect I 


HFGP 


Latent Heat of Vaporization. 




where I * 2. 3, 4 ... N 


TCPO 


Temperature of Condensed Pro- 


WCP 


Weight of Condensed Product 




Water 


Q 


The number of Effects under 




duct Out of Product Cooler 


WP 


Product (Unit) Effect I 


HSR 


Heating Steam Rare 




consideration a * 1, 2. 3 


TSFI 


Temperature of Sea Water Feed 


WPS 


Sum of Unit Products 


HX 


Heat Exchanger 




. . . N 




to Preheater in Effect I 


WVP 


Weight of Va^or Product 



FIGURE 2. BLOCK DIAGRAM, PRELIMINARY PROGRAM 



OVER- ALL CONTROL DATA 


DATA FOR EACH INDIVIDUAL EFFECT (N + 1) 


Number of Effects 


Boiling Point Rise 


Feed to Steam Ratio 


Vapor to Vapor Temperature Difference 


Target Economy Ratio 


Sea Water Preheating Temperature Rise 


Steam Inlet Temperature 


Vent Loss 


Heat Rejection Condensing Temperature 


Heat Loss to Environment 


Specific Heat of Brine 


Boiler Feed Water Return Quantity 


Sea Water Feed Inlet Temperature 


Over- all HT Coefficient for Preheater 


Vent Loss From Final Condenser 


Over- all HT Coefficient for Evaporator 


Raw Sea Water Temperature 


Temperature Rise up to the tube due to pres- 


Temperature Rise of Sea Water in Final 


sure change 


Condenser 




Heating Steam Rate 




or 




Production Rate 





TABLE 3. PRELIMINARY PROGRAM INPUT DATA REQUIREMENTS 

© o © 



as far as total required heat exchanger surface required for given pro- 
duction rate and thermal economy is concerned. Since this surface area 
accounts for approximately 50 percent of the total cost of a plant, very 
rough cost optimization can be manually accomplished. The cost of fuel 
and the applicable fixed charge rate broadly categorize the usable thermal 
economy ratios, as a matter of experience. For the reasons stated above, 
this preliminary program has been very useful and serves quite well as 
the starting point for heat exchanger and vessel design. 

IV. THE INTERMEDIATE PROGRAM 

This program is in the development stage. Additional equations have been, 
written, and modifications to the existing equations have been developed. 
This program becomes possible as a result of the acquisition of sufficient 
operating data to develop temperature related curves for heat transfer 
coefficients and specific heats. These are empirical relationships limited 
to normal sea water feed and an extraction ratio of approximately 0. 67 
(2/3 pound of water extracted from each pound of feed). Minor variations 
of these two conditions can be tolerated without much error. The Block 
Diagram of Figure 3 illustrates the flow pattern for this program. 

The program has five built in optimization loops. Three of these operate 
on process conditions to achieve (1) the optimum preheating profile, (2) 
the correct extraction ratio, and (3) the desired thermal economy. Two 
loops operate on geometry factors to provide (1) the correct number of 
passes and (2) tube length. The requirements for input data, given in 
Table 4 are considerably changed. The time consuming process variable 
inputs have been significantly reduced. The specification of broadly ap- 
plicable geometry considerations as input have been increased. 

The output of this program will largely fulfill the requirements listed in 
Table 2, for items 1, process and 2, geometry, with some exceptions. 
The principal process outputs still missing are (1) an optimized number 
of effects, and (2) electrical and chemical requirements. The geometry 
section will include all the desired information for the input design except, 
the vessel dimensions, pump heads and horsepowers, and vacuum equip- 
ment requirements. The items provided will be a large percentage of 
those required for final design and costing of a proposed plant. The en- 
tire program, with its reduced input data, will provide a much greater 
insight to the engineer for manual optimization, with little, if any, in- 
crease in the time required. This program is expected to be operational 
by September, 1966. 

Computations with the intermediate program will provide a great deal of 
insight into the design of the costing step, and eventually the programmed 
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READ TITLE CARD 
READ HFGP TABLE 

READ CONTROL CARD 



2 



CALCULATE OVER-ALL PROCESS 

PERFORMANCE, TCPI 
ERN, Over-all Economy Ratio. 
XR, Over- all Extraction Ratio 



ESTIMATE REASONABLE SWF 
TEMPERATURE RISE IN PRODUCT 
COOLER AND CALCULATE TSFI (N) 




CALCULATE TCPO 
TEMPERATURE OF PRODUCT 
OUT OF PLANT 



(PASS 2A AND SUBSEQUENT 
\ IF R EQUIRED 



INCREASE OR DECREASE 
TSFI(N) 1° FROM LAST PASS 



ALL PASSES 
PRIOR TO 
ENTERING 
LOOP B 



PASS 1A 

""adjust^ 
"temperature profileT" 

(Restore to Pass 1 Before 
Adiustment) -- ^^ , * , * , 

(PASS 2B 
1 YES (PASS 3B 



PASS 2A ADJUST TSFI @ N/2 
PASS 3 A ADJUST TSFI @ N 



READ N + 1 INPUT CARDS - PASS 1 
CALCULATE NEW PREHEAT 
TEMPERATURE PROFILE 
PASSES 2B, 3B 



PRINT INPUT DATA PASS 1A 
PRINT ADJUSTED INPUT DATA 
FOR ALL SUBSEQUENT 
PASSES FOR VERIFI- 
CATION 



CALCULATE ACTUAL VAPOR, 
SEA WATER FEED, AND 
BRINE TEMPERATURES 



ALL PASSES 
PRIOR TO IB 




PASS IB 



ADJUST SFR DOWN TO INCREASE 
FINAL CONCENTRATION 
ADJUST SFR UP TO DECREASE 
FINAL CONCENTRATION 



SELECT VAPOR ENTHALPIES 
CORRESPONDING TO LATEST 
TVO(I)»b [HFGPm] 



CALCULATE CP(I)'s 
CORRESPONDING TO TBO<I)'s 



CALCULATE PERFORMANCE FOR 
EACH EFFECT (I) OF (N) 



PASS 2C AND 




PASS 3C 




LOOP C 








PRINT MESSAGE 






ERROR MAGNITUDE 






AND ADJUST 






OR EXIT 





CALCULATE UE (I) FOR N EFFECTS 
CALCULATE UH FOR 300 HEATERS 
AND HX200 
CALCULATE UE N + 1 FOR HXRJ 



CALCULATE UNIT DUTIES FOR 

HX200 AND HX's 300+ n 
CALCULATE UNIT DUTIES FOR 
EVAPORATORS AND FINAL 



CALCULATE MEAN TEMPERATURE 
DIFFERENCE FOR HEAT TRANSFER - 
PREHEATERS, PRODUCT COOLER, 
FINAL CONDENSER, AND 
EVAPORATORS 



CALCULATE UNIT HEAT TRANSFER 
SURFACE REQUIRED FOR EACH 
EXCHANGER, CONDENSER 
AND EVAPORATOR 




CALCULATE HSR 



IS HEATING 
STEAM RATE (HSR) OR 
OR PRODUCTION RATE (PR, 
^GIVEN? 

JlFHSR 



CALCULATE PRODUCTION RATE 



CALCULATE TOTAL SURFACE 
AREAS FOR EACH EXCHANGER, 
CONDENSER AND EVAPORATOR 



DOES NOT COMPARE 



PRINT MESSAGE 




CALCULATE NO. TUBES REQUIRED 
TO NEXT HIGHER INTEGER 
FOR EACH EXCHANGER AND 
EVAPORATOR, BUT NOT CONDENSER 



EXIT 
INVESTIGATION 
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FOR EACH PREHEATER 
AND HX200 

ADJUST 
2 PASS TUBE LENGTH? 
(For Preheaters) 



REDUCE TUBE LENGTH ONE 
FOOT FROM PREVIOUS PASS 
EACH TIME AROUND LOOP D 
AND CALCULATE NEW TUBE 
NO. TOTAL 



CALCULATE FVR FOR PREHEATER 
ON 2 PASS FLOW FVR (I), 




FVR (I) » FVR ±0.5 



CALCULATE 2 PASS TUBE SHEET 
DIAMETER 
HTSD(I) 




REDUCE TUBE LENGTH ONE 
HALF FOOT FROM PREVIOUS PASS 
EACH TIME AROUND LOOP E 
AND ADJUST TUBE NUMBERS 



CALCULATE FVR FOR 4 PASS 
FLOW PREHEATER FVR (I) 



FVR (I) is< 
FVR 




(I) = FVR * 0. 5 



PRINT MESSAGE 



CALCULATE 4 PASS TUBE 
SHEET DIAMETER 
HTSD (I) 



Q. 



SUM PREHEATER AREAS INCLUDE 
HX200 

SUM EVAPORATOR AREAS 
INCLUDE FINAL CONDENSER 



CALCULATE EVAPORATOR 
TOTAL SEA WATER 
FEED RATE 




RETURN TO 
START 



NOTE 1 Values to be calculated for each 



Effect 


are: 


WP 


- Product (Unit) Effect I 


BPWT 


- Blowdown Concentration 




WT% 


WPS 


- Sum of Unit Products 




Through This -Effect 


er 


- Individual Effect Econ- 




omy Ratio 


WCP 


- Weight of Condensed 




Product 


WVP 


- Weight Vapor into Effect I 


WBO 


- Weight Brine out of Ef- 




fect I 




LEGEND 



BCF 

BWTP (I) 

er 

ER 

ERN 

FVR 

HFGP 

HSR 

HX 

HXRJ 

I 

MHTL 



N 

PR 

RST 

SWF 

TBO 

TCPI 



TCPO 

TSFI 

TVI 

TVO 

WBO 

WCP 

WP 

WPS 

WVP 

XR 



Blowdown Concentration Factor 
"Percentage of Solids, by Weight, in 

Blowdown from Effect I. 
Individual Effect Economy Ratio 
Assumed or Economy Ratio 
Calculated Net Economy Ratio 
Feed Velocity Range 
Latent Heat of Vaporization, Water 
Heating Steam Rate 
Heat Exchanger 

Heat Rejection Condenser (Final Condenser) 

The Effect under consideration where I = 
2, 3, 4 ... N 

Maximum Preheater Tube Length 

The number of Effects under consideration 
n = 1, 2, 3 . . . N 

Number of Effects for this pass. 

Production Rate 

Raw Sea Water Temperature 

Sea Water Feed Rate 

Temperature of Brine Leaving an Effect 

Temperature of Condensed Product En- 
tering Product Cooler 

Temperature of Condensed Product Out 
of Product Cooler 

Temperature of Sea Water Feed to Pre- 
heater, in Effect I 

Temperature of Vapor Into Effect I 

Temperature of Vapor Out of Effect I 

Weight of Brine Out 

Weight of Condensed Product 

Product (Unit) Effect I 

Sum of Unit Products 

Weight of Vapor Product 

Extraction Ratio 



CALCULATE EVAPORATOR 
TUBE SHEET DIAMETERS 
BASED ON TUBE NO. AND 
DIAMETER 
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FIGURE 3. BLOCK DIAGRAM, INTERMEDIATE PROGRAM 

-1 



© o o 



OVER -ALL CONTROL DATA 



Number of Effects 

Initial Steam to Feed Ratio 

Target Economy Ratio 

Steam Inlet Temperature 

Heat Rejection Condensing Temperature 

Temperature Rise of Sea Water in Final 

Condenser 
Target Preheater Brine Velocity Range 
Inlet Sea Water Concentration 
Vent Loss From Final Condenser 
Raw Sea Water Temperature 
Heating Steam Rate 
or 

Production Rate 
Target Extraction Ratio 
Maximum Preheater Tube Length 
Evaporator Tube Length 



DATA FOR INDIVIDUAL EFFECTS 
Boiling Point Rise 

Vapor to Vapor Temperature Difference 
Sea Water Preheating Temperature Rise 
Vent Loss 

Heat Loss to Environment 

Boiler Feed Water Return Quantity 

Temperature Rise up the tube due to 

pressure change 
Evaporator tube diameter and gauge 



TABLE 4. INTERMEDIATE PROGRAM INPUT DATA REQUIREMENTS 



optimization of the over- all plant, which is presented in the next section. 
V. AN ULTIMATE PROGRAM 

The ultimate program required for design of the optimum plant may be 
represented in a generalized form by Figure 4. The subordinate de- 
sired end results listed previously in Table 1, are not all provided, so 
some manipulation of the input, based on preliminary output will still be 
required for unusual situations. However, the principal desired result 
of minimum water cost will certainly be provided within the accuracy of 
the input data. The subordinate desired results of minimum surface area, 
optimum thermal economy, minimum capital cost, and optimum plant 
life, will be provided by the same design as minimum water cost in a very 
high percent of situations encountered. Sensitiveness of the optimum de- 
sign to fuel cost changes and interest or fixed charge rates can be investi- 
gated very easily by changing only two input data points. 

The availability of reliable process correlations cost data, and the de- 
velopment of realistic optimum search techniques are the two most im- 
portant obstacles to the early completion of a useful ultimate program. 
The small computer user is presently faced with a need to store large 
blocks of empirical data for sea water processing, in the present state 
of the art, which makes it necessary to consider off line storage or other 
difficult and/or expensive means of expanding computer capacity. Cor- 
relations now developed for heat transfer coefficients, such as those in 
use by Houston Research Institute (Reference 2) require complex sub- 
routines and are thus better suited to large computers. 

The author, through the analysis of operational data from the Freeport 
Plant, does not consider any of the present correlations of heat transfer 
or thermodynamic data to be sufficiently reliable to warrant program- 
ing equations instead of storing tabular experience data. Thus, the 
capacity of the IBM 1620 presents another restraint to further develop- 
ment towards the ultimate program. 

Cost data, now being acquired in volume, must be carefully related to 
process and geometry. The geometry section output must be indexed or 
coded for capital costing. Progress is being made in these areas. It 
is expected that the intermediate program will be extended to the cost 
section by the end of Calendar Year 1966. Optimization search techniques 
are also being developed to guide the manual manipulation of input data 
arid thus provide a degree of optimization which is consistent with the 
state of the art. 
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READ CARDS 



PROCESS 



Given Target performance, optimize all inde 
pendant variables within requirements of site. 
Calculate all process flows and conditions. 
Select tube size and calculate Heat Transfer 
Coefficient. 




MAJOR ADDITIONS REQUIRED 
AND REMARKS - 

Dynamic programming required. 
Establish and adjust number of 
effects. 

Establish and adjust temperature 
profiles and preheating ap- 
proach temperature, 
Establish tube sizes (diameters) 
Calculate Heat Transfer Coeffi- 
cient. 

Optimize for minimum surface 
yielding desired performance 
ratio (Gradient Search). 



Dynamic programming required. 
Establish vapor flow spaces and 
areas. 

Establish pipeline sizes and va- 
por flow areas. 

Establish pressure drops and 
size pumps. 

Coding system for cost program, 

Test for geometry acceptance 
and method of instructing pro- 
cess changes. 



ALL NEW 

No dynamic or search proce- 
dures. 



ALL NEW 

Gradient search procedure to 
instruct for geometry alternates 
and then to instruct process 
change pattern in msjor loop. 
Comparative technique must also 
be established. 



FIGURE 4. GENERALIZED BLOCK DIAGRAM, 
ULTIMATE DESIGN PROGRAM 
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Symbol List 

n - An effect of evaporation train containing N effects, where 

n = 1, 2, 3, . . . N. 

W - Mass rate of flow per unit time. 

W f p - Mass rate of steam flow per unit time into an effect. 

Wp - Mass rate of steam flow per unit time out of an effect. 

W f b - Mass rate of brine or sea water flow per unit time into an 
effect. 

Wb - Mass rate of brine or sea water flow per unit time out of 

an effect. 

W ! f v - Mass rate of steam flow per unit time out of a condensate 
flash chamber. 

W*c - The mass rate of condensed water flow into a flash chambe 
FT - Flash chamber for hot condensate. 

HXRJ - Heat Rejection Condenser 

HX30n - The exchanger preheating sea water feed in effect n. 
HX2Q0 - The product cooler-feed preheater. 
W cc - Cold condenser cooling water. 
W wc - Warm condenser cooling water. 



NOTE: All other symbols used in this paper are defined in the legend 
appearing in Figure 3. 
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DO NOT ENTER COMMAS. 

Enter one only, HSR or PR. in 

THOUSANDS. 

An N+1 card must be entered as 

LAST CARD WITH AN EFFECT NO. 
OF N+1 . 

Enter UH-200 on N+1 card. 

Enter UE-HRJ on N+1 card. 

Any decimal point can be over- 
ridden. 

Zero values may be left blank. 



EXPERIENCES WITH FLBSPS AND STOVE 



or They Won't Approve a Faster Computer, What Now? 
By K. W. Jones, Department of Highways, State of Colorado 



The primary purpose of this paper is to pass on three table 
handling techniques used at our installation to speed up 
existing programs as follows: 

1. Optimum instead of Logical Ordering 

2. Random instead of Sequential Storing and Searching. 

3. Pocket instead of Whole Array Sorting 

These techniques were utilized in two programs - FLBSPS - 
"Faster, Less Bulky Symbolic Programming System" (1 and 2) 
written in SPS and - STOVE - "Sampling, Testing and Overall 
Variance Etc." (3) written in FORTRAN. 



Chart I displays the ordering of the operation codes in the 
op code table in SPS II. This table is used to check for 
valid op codes and to branch to the appropriate handling 
routine for each accepted code. Also displayed on Chart I 
are the results of a program written to accumulate the fre- 
quency of occurence of each code. The frequencies shown are 
the results of twenty two source decks which were analyzed 
by this program. It is apparent that the searching time can 
be reduced by a rearrangement of this table noting particu- 
larly TF f TFM, B, DC and DORG and the blocks of unused 
operations* 

Chart II displays the optimized arrangement in FLBSPS and a 
tested speed gain for the first minute on a symbol table read 
for' a patch assembly. The speed gain shown is based only on 
the rearrangement of this table. 



TECHNIQUE ONE 



Optimum Instead of Logical Ordering 




TECHNIQUE TWO 



Random Instead of Sequential Storing and Searching 



Chart III is an attempt to demonstrate how the two above 
methods utilize a table area. The Sequential Method starts 
at the first slot in the table and continues down through 
the table examining each slot until an open slot or dupli- 
cate is found when storing or until the sought for label is 
found when searching. This is demonstrated for the 11th 
label on the left side of the chart by following the arrows. 
Hie Random Address Method uses the actual label being stored 
manipulated in some way to compute the address of some speci- 
fic slot somewhere in the table where the label is stored or, 
if searching, found directly. A problem arises when two 
different labels happen to generate the address of the same 
slot and in this case we revert to a sequential search from 
that slot forward as shown by the arrows on the right side 
of the chart for the 11th entry. 



Chart IV details the program steps used in SPS II to store 
the label table (Pass I) and to search for operand labels 
after the table has been completed (Pass II) using the 
sequential method* 



Chart V displays the methods and results of a program written 
to compute a random address from a label using different 
methods from which we wanted to select the best. This pro- 
gram was run with a 931 label symbol table as the source of 
labels. The method shown in the leftmost column (Addition) 
wherein the 12 numeric digit labels are broken into four 3 
digit parts which are added together and adjusted to fit the 
table area is obviously superior in speed and at least equal 
in the minimum number of duplicate addresses computed and 
therefore was selected. 



Chart VI details the program steps used in FLBSPS to store 
the label table (Pass I) and to search for operand labels 
after the table has been completed (Pass II) utilizing the 
random Address method selected. 



Chart VII displays actual test results of the difference in 
Pass I card reading speeds, again using the symbol table for 
a patch assembly, and shows a maximum speed gain of approxi- 
mately 11 times faster with 700 labels when using the random 
addrsss method of storing instead of the sequential method* 




TECHNIQUE THREE 



Pocket Instead of Whole Array Sorting 



The left column of Chart VIII details statements originally 
written to perform a core sort which, for each argument, 
sequences up through previously sorted arguments until a 
greater value is found and then finds the end of the sorted 
array finally sliding all greater values ahead one slot and 
storing this argument. Core storage requirements and test 
times for 200 numbers are shown and were used as a basis for 
evaluating the following improvements. The right column 
details the second method tried and although no speed improve 
ment is shown the memor requirement is far less than the othe 
methods tried. 



Chart IX shows the methods and results of two additional 
methods tested where a double length array for sorting is 
broken into pockets and each argument is sorted into its 
proper sequence in its appropriate pocket and finally the 
used portions of the pockets are put together to form the 
sorted array. Using 10 pockets we find that storage require- 
ments have more than doubled but we have achieved a time 
requirement of only 0.12 of our original. 



Chart X details the statements used in implementing the 10 
pocket sort and reviews the four basic steps necessary. 



These three techniques are not put forth as the best ways to 
handle similar problems but rather as significant improvement 
aver the first thoughts a programmer might have. 
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SF 
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LDHED 


CM 
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TFM 


85140 




30 


... 






86110 


SFIFIN 


IF 


86120 




TO 


S6130 




TF 


86140 




AM 


86150 




A 


86160 




c 


86170 




BNE 


96180 




C 


86190 


* 




86200 




BE 


87010 




TF 


87020 




AM 


87030 




B. 



LABLE STORING (Pass I) 

ROUTINE TO LOAD LABELS INTO SYMBOL TABLE 
CLERER+11»INPUT+10 

***** COUNT NUMBER OF CHARACTERS ***** 
LAp-H 

***** CHECK EACH FOR VALIDITY ***** # » # 
LABCT.R»6»10 

***** ADO HEADING DIGIT IF IS ONE ***** ♦ 
*+23»SYMTBL START AT FIRST OF TABLE. 

SEIFIN» , IS THIS OPEN SLOT — « — 

***** STORE LABEL THEN EXIT ***** « ## 



*+23»IT+23 
♦♦35 

*+71»*-l 
♦♦59»*10 
*+47»*-l 

LABCTR»SEIFlN+47 
*+36 
INPUT 3 
MULTIPLY DEFINED LABEL 
ER10 



NO» MOVE ADDR OF LABEL* TABLE ) 
LENGTH AND PEEL IT OFF* 
BUILD ADDRESS FOR CHECK IF 
LABEL IS SAME LENGTH* 

IS LABEL SAME LENGTH.' 



YES. IS THIS 



IT+23»*-13 
IT+23*6»10 
IT+12 




NO» THEN INCREMENT TO NEXT 
LABEL AND GO CHECK IT FOR 
A DUPLICATION. 



0.005 second I** I 3 be/ 
0.4 second 100* I tie I 
0.8 second ZOO* / 3 ble 
Z.I seconds 500- hb/e 



# J / Instructions' Executep 
for Evexr Fkeviou*Lr 
StorEP Lable Until Empty 
Slot z$ Founo. 

4k Z Aopitional Instructions? 
Executed for Each Label 
FaS*£p of Sams Iength. 



N 



LABLE 5EARCWING*(passII) 



35020 


* 






35030 


* 




THE SYMBOL TABLE IS 


35040 


» 






35050 


LBADD 


TFM 


-•+23.SYMTBL 


35060 




BD 


LBADDS 


35070 


*c 




UNDEFINED SYMBOL IN 


35030 




8TM 


EVALER»50000 


35090 




DC 


1.-.* 


35100 


LBADDS 


TF 


*+23*LBADD+23 


3 5110 




TD 


* + 3 5 


35120 




TF 


LA8C0M+11.*-! 


35130 




TFM 


♦♦47» »10 


35140 




A 


*+35»«-l 


35150 




A 


LABC0M+l»*+23 


35160 




CM 


C0LL-17 


35170 




BNE 


*+36' 


35180 


LABCOM 


C 


COLL-2 


35190 




BE 


LABOK 


352 OG 




*w 


» » . « •» _ * ~% t\ 


36010 




TP 


LBADD+23 »LABC0M+1 1 






o 





START AT 1ST OF TABLE. 
HAS WHOLE TABLE BEEN LOOKED AT. 
OPERAND 
YES. THEN TYPE ERROR 5. 

NO* MOVE ADDR OF LABEL (TABLE) ^ 

LENGTH AND PEEL IT OFF, 
MOVE ADDR AGAIN FOR 3 DOWN. 
STORE LENGTH THEN 
DOUBLE IT. 

FINISH BUILDING ADDR LABEL ( TABLE ) . 
IS LENGTH OF LABEL LOOKING FOR 

SAME AS THIS LABEL ( TABLE ) . 

YES. IS THIS THE LABLE WE 

ARE LOOKING FOR. 

KC? INCREMENT TO ADDR OF NEXT 

LABEL LENGTH* MOVE It AND GO 
CHECK NEXT LABEL (TA fj # 



0.005 second /-search 
OS second /00* search 
0.9 second 200- search 
Z.3 seconds 500- search 



j t2 Instruction Executep 
for Each Lable Storep 
Aheap of the LablE 3ein& 
Searched For. m± 



CHART 3C 





ETHODS &. KESULTS 



K ~ . ■ — 

Addition 


Quotient 


Combined (&cr,AJJ) 


ReMAXN-DEfc 


xxxxxxxxxxxx 

A ~B C D 


xxxxxxxxxxxx 

OVR 


xxxxxxxxxxxx 

LBL 


xxxxxxxxxxxx 

00 


[(a+b<+d)*l]+s 


[(ft) - l ] + s 


[(LBL+K)*lUS, 






MlLLE SECONDS PI 

If) £ 


r Label Calc. 




Average Numser 

L33. 


>f Duplicate At 

7£*> Many 
to CouA/r 


'DRES5ES (93 J L 


\belS - 1000 Suits) 


L= Lenglli of fable slot K " Constant to prevent bunching 
S - SJarf stldress of table NS =- Number, of slot* in table 



*MDOM ADDRESS ffl 



CHART 3ZT 




ft*? 00 
87010 • 
87020 • 

87030 LDLBL C 

♦ • • 

87110 zspec tf 

• • • 

88020 LOHfAO BNF 



LABEL 5T0RING (PASS n 



ITS 



* « * 

88120 
88130 
88140 
88130 GRP1 
88160 



ROUTINE TO LOAD LABELS INTO SYMBOL TABLE 

CLERER+1 1 t INPUT+10* • 
»••»• COUNT NUMBER OF CHARACTERS 
LABtlNPUTMO** 

CHECK EACH FOR VALIDITY ♦»♦#• 
NOHEDt 1SHE0* • 
*»•»» ADD HEAOING DIGIT IF IS ONE •••«• 
TR ARGL0C-4tGRPl-4*» GO TO SUBROUTINE TO FIND 

B GETLBLt »* TABLE SLOT AND CHECK FOR 

DORG ♦-4»»t DUPLICATE LABEL. 

DS A LAB* STOLBL * ER 10 • ER9 
DC !•- 



.008 second I s * febe I 
. 010 second avenge label 



33050 * 

35060 • 
35070 * 
35080 LBADD 
33090 
33100 

35110 MVLBLB 
35120 
33130 
35140 
,3130 GRP2 
160 



LABEL 5EARCHING (PASSU") 



THE SYMBOL TABLE IS SEARCHED FOR EQUIVALENCE 



TF LB.CLERER+11 • t 

TFM MVLBLB+6*LB-12** 

A MVLBLB+6*C0LL-1? 
TF •C0LL-2 

TR ARGLOC-4*GRP2-4** 

B GETLBL* »• 

DORG »-4*** 

DSA LB • NOTDEF tGOODLB *NOTDEF 

DC lt- 



BUILD 
DIGIT 



A 12 
LABEL. 



GO SEE IF 
LABEL HAS 
BEEN DEFINED. 



.01 Z second j < d /tiel 
.015 second sverjge label 



SUBROUTINE COMMOM TO BOTH ABOVE 



Q 



89030 
89040 
69030 
89060 
69070 
69080 
89090 
69100 
89110 
69120 
89125 
89130 
89140 
89150 
89160 
89170 
89180 
89190 
89200 
90010 
90020 
90030 
90040 
90050 
90060 
90070 
90060 
90090 
90100 
90110 
90120 
90110 
90140 
150 
160 
"0170 
90180 
90190 
90200 



* 
* 
* 
• 

ARGLOC 
FREXJT 
DPEXtT 
TILT 

GtTLiL 



SUBROUTINE TO SEARCH LABEL TABLE FOR SAME 
LABEL OR UNUSED SLOT OR NEITHER. 



ISFREE 

606K 

NXTSLT 
LOOPS* 



DS 

OS 

DS 

DS 

DS 

TF 

CF 

SF 

SF 

A 

A 

A 

SF 

MM 

AM 

TF 

AM 

TDM 

BNR 

B 

DORG 

AM 

C 

BE 

NOP 

CM 

BL 

TFM 

TOM 

AM 

B 

DORG 
CM 
BWC 
8 



WK*ARGLOC*ll* 






MAKE PLUS* 




PARTITION 


WK-5*»* 


INTO 4 


WK*WK-9** 


SECTIONS THEN 


WK*WK-6»» 


ADD THEM UP 


WK*WK-3»* 


TO COMPUTE 


WK-2*** 


SYMBOL TABLE 


WK»17*10* 


ADDRESS. 


99*20000 




L00PBK+11*99*» 


SET LOOKED 


LOOPBK+lltll** 


THRU WHOLE 


GOBK*l*l*t 


TABLE CONTROL. 


*+20*99*il* 


IS THIS LOCATION FREE. J 


FREXIT*»6» 






#-3 

99#U»* 






99 ♦ ARGLOC #611 # 
OPEXIT06 , 


LABEL IDENTICAL. I 




LOOPBK»M 


NO* M ! WRAPPING AROUND* 


99» 9973»» 


, NO* SHOULD I START 


NXTSLT 


WRAPPING AROUND* 


99*1999** » 


1 YES* 


G0BK+1«9 


99*6** 


NO* THEN JUST GO 


tSFREEMi 


ON TO NEXT SLOT* 


•-3 




99 » M 


« f BACK AT 


NXTSLT 


STARTING POINT* 


TILT»#6* 


YES* 



><S>A 







/7* /.3 m ± 22 Average 
Number of Instructions 

Executed No Matter How 
Man/ Label* Were 
Previously Storeo. 

±23* 1.3* ± 30 Average 
Number of Instructions 
Executed Mo Matter 
Where in the. Table we 
Label Being Looked For 
-Was Stored. 
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ABM 50ST MIKE S SORT 



29 



DIMENSION 
FRSORT(l) 



STOI 400*4)* 
99999999* 



NN » NBRSET(I) 
DO 35 J * It NN 



FRSORT<401> 



Z SORT THE WHOLE ARRAY 
Z THAT WAS READ IN* 



30 
31 
32 
33 
34 



00 30 K • Iff 400 Z LOOK FOR GREATER OR 

IF I FRSORT < K ) - STO(J*IJ) 30O1.31 Z EQUAL* 

CONTINUE 

^ 

00 32 L * K*400 Z NOW FIND TOP OF 

IF I FRSORT ID - 99999999*) 32*33*94 Z TABLE* 

CONTINUE 

FRSORT (L+D - FRSORT ( L ) Z NEXT SLIDE WHOLE 



IF (L - K) 94*35*34 
L * L - 1 
GO TO 33 

FRSORT (L ) = STO(J*I> 



Z TABLE UP ONE SLOT 
Z TO OPEN SPACE FOR 
Z THIS ENTRY. 

Z NOW STORE THIS AND 
GO DO NEXT IF ANY. 




S/tic Tzi>/e 
up One Notch 
Men Jterc 
Me* Argument 




OOO0022222 
00QQ03 33 3 3 
0000044444 



0000055555 
0000077777 
0000088888 



0000099999 
9999999999 



Looking for 
► Proper S/or 



Finding BnJ 
f of TaB/e 



c 

C 

c 
c 
c 



558 

559 



DIMENSION STOf400»4) 
MM « NBRSETt I ) -I 
LIMIT » NBRSETin 
DO 559 MIKE - 1* MM 
LIMIT • LIMIT - 1 



2 SOffT ALL *E*D IN* 

2 GO CLEAR TO TOP 1ST 

TIME WITH LARGEST* 
2 SORT WHOLE ARRAY. 



Z GO 1 LESS FAR THIS 
TIME. 

DO 559 MIKE1 * 1* LIMIT Z IS THIS VALUE GREATER 
IF (ST0<MIKE1*II - ST0(MIICE1*1*I) 559*559.558 

THAN THE NEXT ONE* 
HOLD » ST0IMIKE1*I) Z YES* JHEN INVERT 

ST0(MIKE1*I) » STO(MIKEl*l*n Z THESE 

ST0CMIKE1+1*I> » HOLD Z TWO SLOTS* 

CONTINUE 



/- Move largest value to top of if 
array exchang/ny inverted 
pairs on the way. 

Z.~ Move next largest value to 
1 zlot below the top of the 
array exchanging inverted 

pairs on the way . 

ETC.- down to the sma/test 
value in the first slot . 



Requires ± 970 
Core Locations 
Plus Extra 

rt &AND 



16 MINUTES 

Sort Time 
for 200 

ITEMS 



Requires ±610 
Core Locations 

(U$e$ Same Bmp) 



17 Minute? 
Sort Time 
for ZOO 
Items 
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4 Pockets 



1 



I<r~ Initialize 
Limits o-f Pockets 

from Range of Data. 
n FOR EACH ARGUMENT: 

A. 

XEo~ Determine which 
Pockef to Put -Anj. In. 

Mo" Sort Arg. inlo 
Proper Place In 
Pockef. 

AFTER ALL DONE : 



§ BL Slide Pocket* 
-(g ). j , TogeMier info'l.anJ 
2 for Sorted Array. 



RANGE « 00OOO TO 999^9 
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10 Pockets 

t 

The same four 
steps detailed on 
the left are used 
here except a 
"Normal Distribution 
should be used 
for setting the 
limits of the pockets 
under step I. 

This Method 
Required Oml f 

'0.12 OF THE 

Original 
Sort Times 





4 l/Z MINUTES 


Require* ± 37DO 


2. Minutes 


Cofce Location^ 




Cof?e Loo*tions> 


^ORT Tine 


Fiog Extka Pou^i.e 


R>* zoo 


Flu£ Ext*a Double 


Fo* Zoo 


Length *dokt Bawp 




Length £ort 3ahP 





VQSORT DETAILED 
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DIMENSION STOUOO.M* PSI800I 



I FS SAME AS FRSOUT* 



900 



*0l 
504 

sot 



910 



512 
5U 

515 



516 

61? 
717 



517 
518 

519 
520 



5*5 



o 



10 
5* 



THF DATA IS SONTCD INTO TEN POCKETS THEN THESE POCKETS 
ARE PUT TOGETHER INTO THE FINAL. SORTED ARRAY* 



D! MFNS ION IFAR(IO) 
DO 500 IANYTM ■ 1 t 10 
IFARMANYTM* - IANYTM • 



Z FARTHEST WAY GONE IN EACH POCKET • 



60 - 79 



Z INITIALIZE TO FIRST SLOT EACH ONE 



7Z> 



DIMENSION B0TH1I Z CALCULATE MINIMUM VALUE FOR EACH POCKET, 
*NG • SAVMAX(I) - SAVMINU) 



ROT« 
«OTI 
WOTC 
BOT( 
ROT< 
BOTI 
KOT( 
BOTI 
»OT( 
POT t 10) 
B0T<11> 



1) 
2) 
3) 
41 
51 
6) 
7» 
8) 
9t 



SAVMINMI 
SAVMINi I I 
SAVMINf II 
SAVMINU! 
SAVMINU* 
SAVMI N( I I 
SAVMINU* 
SAVMI N( 1 1 
SAVMINU* 
SAVMI N( I * 
SAVMINU) 



IXENO • NBRSETU I 
DO 520 IX It IXENO 
ARQ * STOUX*!) 



•24»RNG 
•35*RNG 
•40«RNG 
• 4MRNG 
*50*RNG 
•55«RNG 
•60»RNG 
•65«RN6 
•?6«RNG 

RNG 4 I (SAVMlN 1 ! **RNG J • 

Z SORT THE WHOLE 
2 ARRAY READ IN* 



Z LIMITS REVISED AFTER TEST DUMP* 



.1) 



IPOC "4 2 DETERMINE WHICH POCKET THIS GOES IN* 

IF ( ARG - BOTIIPOO) 301*504*304 2 (START IN MIDDLE) 

IPOC • 1POC - 1 

IF I ARG - BOTUPOCU 501*310*510 
IPOC « IPOC ♦ 1 

IF i ARG - BOTUPOC)* 309.304*304 

IPOC • IPOC -1 

LOOK THRU THIS POCKET AND PLACE *N PROPER POSITION IN ORDER MAKING 
ALLOWENCE FOR OVERFLOW INTO NEXT POCKET IF THIS ONE IS FULL* 



2 SET BEGINNING OF LOOK* 



SET LAST SLOT PREVIOUSLY USED 
AND INCREMENT FOR NEXT TIME* 



I BEG • IPOC * 80 - 79 
tENDPC * IPOC • 60 
1FR2G0 • IFAR(IPOC) 
IFARUPOC* » IFARtlPOCI ♦ 1 
DO 514 IPO * TBEGt IENDPC 

If UPO - IFR260) 512*319*94 
1FIARG - FSUPOH 516*316*514 

CONTINUE 

IF (IPOC - 10) 515*96*96 
IPOC « IPOC ♦ 1 
GO TO 510 

NEXT SLOT VALUE IS GREATER THAN ARGUMENT SO SLIDE REST 
OF POCKET FORWARD LEAVING A SPACE FOR THIS ARGUMENT* 



2 AM I AT PREVIOUS LIMIT* 
2 NO* IS THIS ARG GREATER THAN SLOT. 
2 NO* GO UP TO NEXT SLOT IF ANY LEFT* 
Z NONE LEFT THIS POCKET SO GO TO 
I NEXT ONE UNLESS IT IS THE LAST ONE* 
— 



IF UFR2GO - 801) 617*96*94 

IF (IFR2GO - UPOC»80*1*J §lt#?lt»n7 

IPOC ■ IPOC ♦ 1 

IFR2GO • IFARUPOC) 

IFARUPOC) * IFARUPOC) ♦ 1 

GO TO 516 



2 FIRST IS STORAGE EXCEEDED. 
2 NO* THED DID POCKET OVERFLOW* 
2*YES* THEN GO UP TO END OF 
2 NEXT POCKET BEFORE BEGINNING 
2 SLIDE* 



FSUFR2GO) • FSUFR260-1) 

IF UFR2G0-1 - IPO) 94*519*518 

IFR2GO ■ IFR2G0 - 1 

GO TO 517 

FSUPO) • ARG 

CONTINUE 



2 SLIDE PORTION OF POCKET UP 

2 MAKING ROOM FOR THIS ARGUMENT. 



Z STO THIS ARGUMENT 
Z AND GO GET NEXT ONE. 



PUT ALL POCKETS TOGETHER IN THE LOW END OP THE *FS* ARRAY 



!T0 • IFARU) 
DO 539 IPMV - 2* 10 
I FROM • IPMV •60-79 

IF UTO - I FROM) 534*531*551 
!T0 ■ 1FARUPMV) 
00 TO 539 

IF UFARUPMV) •» IMGM) 949359*339 
FSUTO) ■ FSUPROH* 
!TO • 1TO *1 

ifrom • trno* ♦ t * 

f P UFARUPMV) ■*> fP*0*t 94*539t3l6 
fp HtMWMO+11 - I FROM) 94*999*599 
CdNTtNUf 



Z MOVE POCKETS 2 THRU 10 



Z 010 THIS POCKET FILL UP OR OVERFLOW 
Z YES* SKIP TO NEXT ONE THEN. 

Z NO* WAS THIS POCKET EVEN USED* 
Z YES* SLIDE THIS POCKET DOWN. 



Z ft WE PAST LAST SLOT USED* 

Z NO* ft WE PAST LAST SLOT TABLE* 



i 



1. Initialize 

Set location uled irW far in 
each pocket to * None ? 



Compute and set Minimum 
Value 4o be placed in each 
pocket. 

2. Which Pocket? 

Pick up nexf argument. 

Determine which pocket this 
argument shoulJ be Sort* J into 
s far ting at center tor |ea*t 
tries * 

3. Sort Into Pocket 

Determine where in pock*4/* 
fhi* argument belong* 

(Sirr?)/3r io the, \ 
"5-trortg Arm * Met ho J J 

51 iJ* yrester valued in 
pocket/ s forward one slot 
Uaving B hole 4or argument. 

Store 4hi^ argument gni 30 
get next one. 

4.. Slide Pockets 

TbGETHBMNTO 

One Sorted 
Army 



77 
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STUDENT SCHEDULING ON THE 1620 



by 

George L. Crumley 
The Citadel 
The Military College of South Carolina 
Charleston, South Carolina 

Presented at 
The Western Region Meeting of COMMON 
Denver, Colorado 
July 6, 7, 8, 1966 

At the beginning of the summer vacation last year, we began 
thinking seriously about a computer registration program for The 
Citadel. This v/as undertaken strictly as an academic problem, to see 
whether or not we could do it. There was no pressure, or even approval 
at the beginning, from the college administration. 

We began by looking at a couple of programs used by other 
colleges- After considering traditions and certain established ways 
of doing things at The Citadel, we decided that it would be easier 
to write our own program than to try to adapt some existing program 
to our needs . We tried to arrange our disc files and course numbering 
systems to agree with existing systems and procedures, so that the 
faculty and administration would not have to accept a complete 
revision of all their habits. We decided also, near the beginning, 
that conditions at The Citadel were such that we did not need to place 
any particular emphasis on speed of the program. Dealing with 
approximately 2000 cadets, we were not under any great pressure from 
the point of view of time limitations. Therefore, whenever there was 
any question of doing a job thoroughly or doing it quickly, we elected 
to do it thoroughly. 

We have a Model 1 1620 computer with 20K memory, one 1311 
disc drive, and a 1622 Card Reader- Punch. Me also have a 407 accounting 
machine which we use for printing card output. 

We were able to devote an entire disc pack to the scheduling 
program, and began by loading Monitor I with the SPS assembler and the 
Disk Utility Program. 

The scheduling system consists of a number of programs stored 
on the disc, all operating under the control of a supervisor called 
EXEC. It and the Monitor input/output routines remain in core at all 
times. The input/ output routines are necessary because we have used 
the GET and PUT statements for disc operations. All programs in the 
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scheduling system arc called by use of EXEC control cards. Input can 
be stacked, except in a few special cases , where the last card indicator 
is used. 

The master schedule is stored in two forms in different 
places on the disc. The full schedule for every section of every course 
is stored as a two-sector record. This record is used primarily for 
printing the cadet schedules and the master Schedule. A short form of 
the section record- -giving only the section identification, its 
meeting times in coded form, and the number of seats available and 
already assigned— is stored as a one-sector record. This short record 
is used in the scheduling program itself. Additionally, we have 
provided a course directory so that the record of any course can be 
located without a complete search. 

Each cadet has a two-sector record showing his identification 
number, name, his academic major, his class, and his advisor's code 
number. This is followed by a list of the sections in which he is 
currently enrolled. There is room for him to be registered for 12 
courses, which is certainly more than the Registrar would ever allow 
him to take in any one semester. Each course that he is taking is 
represented by a seven-digit entry. Finally, there is an alphameric 
date which shows the last time that this schedule was changed and a 
two-digit number which starts at zero with the cadet's original regis- 
tration and is increased by one every time a change is made in his 
registration. 

The directory that we have provided for the cadet records 
serves two purposes. First, when we want to find the record of a 
particular cadet, we can locate it without searching all 2000 records. 
Second, because we know at all times who is registered, we cannot 
register the same cadet twice. If we attempt to do this, the act of 
making a new entry in the directory gives us error messages. This 
directory is kept in order by cadet number so that we can produce 
various outputs alphabetically by going through the entire directory. 

There is a short faculty directory so that instructor or 
advisor codes can be translated into alphameric form for printing on 
cadet schedules, the master schedule, or other output. 

To begin each semester's registration, department heads are 
asked to send the master schedule for their departments to the proper 
administrative office. There the schedules are completed and coordinated, 
and are finally punched on cards. By loading the master schedule in the 
computer, the cadet directories from the last semester are cleared, and 
both the long and short forms of the section records, as well as the 
course directory, are stored on the disc. 
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By using DUMPSK [Dump Schedules] , we arc able to retrieve the 
full form of the schedule which was loaded, as well as the name of the 
instructor of each section, who was identified by code number. These 
cards are used to cut mimeograph stencils on the 407 so that the full 
master schedule can be .reproduced. Using SHSKD tshort Schedule], we 
can obtain cards for cutting the stencils for the List of Courses, 
k'e have found that both the long and short forins of the Master Schedule 
have their uses. The average student in enrolling is interested merely 
in getting the code numbers for the courses for which he wants to 
register— the short version of the schedule is best for his purposes. 

We can pre-pUnch Course Request Cards for all students who are 
expected to register the following semester, and these cards, together 
with the master schedules 5 are sent to each cadet's raculty advisor. 

Each cadet, in conference with his advisor, fills out a Course 
Request Card, using a five-digit number to identify each course that he 
wants to take. When these cards are received by the Computer Center, 
they are key punched, and first run through a program called COUMTM.* 
This program makes a t^lly of the number of requests fbr each tourse 
and gives us a typewriter message every time a cadet requests a rton- 
exi stent course. It also punches a blank Course Request Card for those 
Cadets who ask for nonexistent courses. The outpiit from C0UHTT1 gives 
us the number of seats available for each course, the number of requests 
for each course, and then the number of seats still available, or the 
number of seats short. He examine the original cards of those cadets 
who asked for nonexistent courses and find that some of these were 
caused by key-punching errors. Usually we find that the error was 
made by the cadet. In some cases, we can assume what he probably meant, 
and we make out a new Course Request Card for him. In other cases, we 
have no idea what he meant, so we merely delete his request in the 
corrected Course Request Card. Generally we run COUNT?] again after 
making corrections. The information concerning the number of seats over 
or short is sent to department heads, who may wish to make changes in 
the master schedule as a result. 

The next step is to sort the Course Request Cards by class, 
putting seniors first, juniors next, and so on, and then run them under 
SKEDCA [Schedule Cadets], our primary scheduling program. SKEDCA sets up 
an array of available sections from the courses requested. First, it 
ranks the courses requested so that the course with the fewest number of 
sections open receives the highest priority and the course with the 
greatest number of sections open receives the lowest priority. After this 
ranking, the sections within each course are ranked so that the section 
with the greatest number of seats available will be examined first, and 
that with the fewest number of seats available will be tried last. After 
setting up the array, the computer begins to examine every possible 
combination of sections to find a schedule without conflicts. As soon 
as such a combination is found, the cadet record is placed on the 
disc, showing the courses and sections for which he is registered, 
and the section records are adjusted accordingly, the number of seats 
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assigned being increased by one, and the number of seats available 
being decreased by one. Placing the revised section records on the 
disc completes this cadet's registration. His schedule in printed 
form can be obtained at any later time. 

If all possible combinations of the sections available are 
tried without a nonconflicting schedule resulting, the computer then 
types a conflict message. Because of the way our program is written , 
we feel confident that when our computer declares a conflict to exist, no 
one on the faculty can produce a schedule by hand that will woH<. 
We kflow of programs used in some colleges that will try 300-500 or some 
other fixed number of combinations, and then declare a conflict. In a 
large college with a limited time for scheduling, we can see reasons 
for following such a procedure— but not in our case. SKEDCA can 
schedule 2000 cadets in about eight hours , or one cadet approximately 
every fifteen seconds. 

When the initial enrollments have been completed under SKEDCA, 
we can produce a summary of the enrollment statistics and we can Use 
OUTPUT to produce cadet schedules. These enrollment statistics may give 
department heads reason to change their schedules either by adding or 
deleting certain sections. We have a program called ADDSEC [Add a 
Section] which can be used to change the master schedule. 

When the first schedules are received by the cadets, inevitably 
there will be, for various reasons, required changes. We have two 
programs available to make changes. The first one, called DROPAD, uses 
a modified form of SKEDCA to completely revise the cadet's schedule He 
is first dropped from every section in which he ifc enrolled, and then a 
list of courses requested is built up, first using the courses just 
dropped with the first Add request. We know that there is room for him 
in all courses except the new one, because he was in them a few micro- 
seconds ago. The computer goes through the same process it did under 
SKEDCA in an attempt to find a schedule with the first additional course. 
If it is successful in this attempt the course just added becomes a 
permanent part of his registration, and the second course to be added is 
tried. If no schedule is found with the last course added, then that 
course is deleted, and the cadet is given a schedule without it. This 
procedure is followed for each course to be added. 

We have found that the procedure used by SKEDCA usually gives 
enrollments which are well balanced among the several sections of a 
given course. In those few cases where the distribution of cadets 
Is noticeably uneven, it seems that the fault lies with the master 
schedule. The sections have not been offered at the times when they 
are needed. 

We have also found that there is a tendency for cadets to be 
quite optimistic in any pre- registration. Very few of them predict 
that they will fail any course. Therefore when the grades come out at 
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the end of the semester, there is a general tendency for most changes 
to be from higher to lower numbered courses. This one-way shifting 
keeps the lower numbered courses well balanced, but sometimes leaves 
imbalances in those courses from which there was a relatively large 
exodus. 

We can rebalance such sections by running blank Drop-Add 
Cards for the cadets in those courses. 

This DROPAD procedure can be used until and on Registration . 
Day. After Registration Day, it is not desirable, of course , to change 
the cadet's entire schedule unless it is absolutely necessary. We 
use the second program, called MAHUDA [Manual Drop-Add], for changes 
in schedule after classes have begun. 

After Registration Day, we turn control of section enrollments 
over to the departments concerned. When a cadet wants to drop or add 
a course, he and his advisor, having determined what section he wants, 
fill out a Section Change Card, which must be signed by the head of the 
department concerned. By signing this card, the department head gives 
us the authority to register" the -student in the section specified. 
MAMUDA does this, without changing any other part of the cadet's 
schedule, and regardless of whether or not the section is still open. 
P.ANUDA does not even check to determine whether or not a Schedule con- 
flict is created by such a change, for we assume that the cadet and his 
advisor have done this. 

We do not check in this instance because the primary concern 
here is that we want to give complete control to the department heads 
concerned. We have found rare cases where, by special arrangements, 
cadets were to be scheduled for courses that appear to conflict. 

We have an artificial department, number 75, which offers non- 
credit courses designed to help give athletes light schedules in the 
afternoon. We set up schedules for three sections of one of these 
courses, the first of which has the greatest number of seats available, 
and whose time schedule takes up most of the late afternoon; the second 
has a smaller number of seats available and takes up somewhat less 
afternoon time and the third has even fewer seats available, but does 
not take up any time at all. By trying to fit these "classes" into the 
cadet's schedule, SKEDCA will schedule his real classes in the morning 
hours if it is possible. In the end, however, he may be scheduled for 
the third section, which keeps us from getting a 'conflict. In this 
manner, we attempt to save the athlete's afternoons for practice but 
we will not sacrifice an enrollment in any academic class in order to 
do so. 

The possibilities for misuse of this device are obvious to 
anyone who is familiar with the resourcefulness of college students. 
The course number for this course is not published in the master 
schedule and is not printed on the schedules of the cadets who are 
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enrolled in it. The Director of Athletics sends a list of eligible 
cadets to the Registrar, who then fills out Drop-Add Cards for the cadets 
involved. To provide greater security, the number of this course is 
changed from one semester to the next. In the event of a "leak", it 
could even be changed during an enrollment, if necessary. 

At the present time this is the only use we make of our 
Department 75. Other colleges could use it to schedule a lunch hour, 
or to free certain hours for outside work. 

At the end of Registration Day, we punch Course Cards for every 
course (except those in Department 75) in which every cadet is enrolled. 
These cards are sorted by the Sorting code, and then are used to print 
section rosters for the instructors* 

After Registration Day there will be a certain number of 
changes by Section Change Cards. When the Registrar has stopped such 
changes, we punch a new set of cards, called Grade Cards, which arc 
printed differently but are punched exactly the same as Course Cards. 
These cards are sorted by sections and are used to provide final section 
rosters for the instructors. The same cards are kept for use as Grade 
Cards at mid-semester s and copies of them are made for use as Grade Cards 
at the end of the semester. 

Since the enrollment of every cadet at The Citadel is on the 
computer disc throughout the semester, we are able to produce rosters 
of the entire Corps of Cadets for the Registrar's Office, rosters of 
cadets majoring in their departments for all department heads, and almost 
any kind of roster for anyone who has a legitimate use for it, at almost 
any time. These rosters are automatically kept up to date as the 
result of searching cadet records for desired information. The cadet 
records are kept up to date by the use of Drop-Add or Delete Cards. 

Changes of major are completed by sending a Change of Major 
Request Card to the computer. The program, MAJORX, looks up the former 
major and advisor from the cadet record and then substitutes the new 
information on the cadet record. Then it produces cards for printing 
notices to be sent to everyone concerned. 

Near the end of a semester, when we are ready to begin pre- 
registration for the following semester, we can obtain pre-punched 
Course Request Cards. 

Since this set of programs has been written as an academic 
exercise and not as a regular part of our duties, we have been under 
no particular obligation to provide documentation. We have drawn 
flow charts for only a few programs where they seemed unavoidable. 
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We have written a scheduling manual which has been distributed 
to the faculty at The Citadel and to those attending the COilslOM meeting 
in Denver. The manual was designed to acquaint The Citadel faculty with 
the fundamental steps in planning and setting up courses and sections 
for the coming semester 5 as well as to provide general and specific 
guidance concerning details of registration and related procedures* 
It also contains operating procedures for Computer Center personnel. 
This manual is not printed as a part of the Proceedings, but may be 
obtained by writing to 



George L. Crumley 
Box 79, The Citadel 
Charleston, S. C. 29409 




WHITMAN COLLEGE REGISTRATION AND GRADING SYSTEM 



Whitman College is an independent four year liberal arts insti- 
tution with an enrollment of approximately 1100 students. As part . 
of its commitment to high quality education, Whitman College established 
in the spring of 1964 a computer center for faculty, student, and 
administrative use* The : original computer installation consisted of 
an IBM 1620 computer with a 1622 card reader-punch, sorter, keypunch, 
and a 407 accounting machine. The computer installation has recently 
acquired a 1311 disk drive, a 1443 on-line printer, a gangpunch de- 
vice on the original keypunch and an interpreting keypunch* 

The first version of the registration and grading system used 
the original computer configuration. A Class Card and a Student 
Master File were designed and two programs were written to do most 
of the desired work* The registration information program printed ' 
student information such as his schedule of classes, a ddresses, data 
on parents, etc* for the use of other administrative offices. The 
grade processing program printed grade reports, probation lists, honors 
lists, rank in class, and fraternity-sorority grades. Other opera- 
tions such as class lists were done using only the 407. The results 
of the initial system were encouraging — a large amount of work was 
being performed satisfactorily on the computer. There were some 
limitations in the system because the original computer equipment 
was inadequate* Without the on-line printer it was necessary to 
punch cards and carry the cards to the 407 for printing. This be- 
came a problem when the number of print lines in the registration 
information report reached 25,000 lines and required the punching 
of 50,000 cards. Both the registration information program and the 
grade processing program required the merging of the Student Master 
File with the Class Cards* This is a cumbersome operation without 
a collator. 

When the Computer Center acquired additional equipment, an^. 
improved registration and grading system was made possible. In planning 
for the improved system it was necessary to decide whether or not 
the Student Master File should be permanently stored on a disk pack 
and if so, whether the File should be stored on a monitor pack or be 
independent of the monitor. Programs could be written which use the 
disk drive and which are independent of the monitor disk input-output 
routines, but this would have complicated the programming for the 
disk file. Therefore, it was decided to use the monitor disk input- 
output routines, thus requiring storage of the Student Master I? ile 
on a monitor disk pack* 

The Student Master File contains a large amount of alphabetic 
information* This combined with the two-digit alphabetic representa- 
tion of the 1620 means that the disk records will be of large size* 
Presently the disk records of the Student Master File are 16 sectors 
or 1600 digits of disk memory* Since Whitman College has an enroll- 
ment of approximately 1100 students, the disk form of the Student 
Master File requires 1,760,000 digits of disk memory* In addition 
there must be memory available for expansion of the disk record* 




Since there are only 2,000,000 digits of memory on a disk pack 
and the monitor program requires about 440,000 digits of memory, 
approximately 1,560,000 digits of disk memory are available for data 
storage. Therefore, more than one disk pack is required to store the 
disk form of the Student Master File* Any data file using more than 
one disk pack on a single drive system cannot be handled readily as a 
permanent disk stored data file because the programming of maintenance 
programs for new records would be too complex. 

To get around the difficulties of a permanent disk stored file, 
it was decided to make all corrections, additions, and deletions to a 
card form of the Student Master File and to load the card form onto 
a disk pack prior to processing registration information or grades. 
This eliminates the need for a collator to merge Student Itfaster File 
and Class Cards for processing. In practice card handling has almost 
been eliminated. Also selective runs on various sets of class cards 
can be made with unused disk master records being bypassed. 
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STUDENT REGISTRATION 




Class 
Header 
Cards 



c 



Start 



3 



Student Header 
Program 



Student Header 
Cards are merged 
manually in front 
of class cards for 
classes in which 
a student enrolls 








Gangpunch student 
information from 
student header card 
into class cards 



Remove Student Header 
Cards from merged deck 




Class List 
Program 



c 



End 




p, Merged 
Student 
Header and 
Class Cards 



Merged deck 
with student 
information 
in class cards 



Class Cards 
with student 
information 



Sort Class Cards by 
student sequence 
number and then 
by class 


4- - 


i 






— 



Class cards 
sorted by 
class 




PREPARATION OF CLASS SCHEDULE AND CLASS CARDS 
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SYSTEM SUMMARY 



I. 



Data Files 



Glass Master File 
Class Header File 
Student Master File 



d # Student Reader File 

II« Reports 

a* Class Lists 

b # Registration Information 

c w Grade Reports 

c a Probation Lists 

e. Honors Lists 

■Z+ Fraternity and Sorority Grade Lists 

£«, Rank in Class 

III* Boo tenants 

a. Class Schedule 

b* College Directory 

IV* Programs 

a. Class Schedule Program 

b„ Cl a ss Heading Program 

c # Class Card Program 

d # Student Header Program 

e. Class List Program 

fi # Disk Load Program 

g # Registration Information Program 

fc. College Directory Program 

i. Grade Processing Program 




CLASS MASTER -FILE 



There are master cards for each class and each laboratory section. 
The comments card is used for printing of the class schedule and 
applies to both classes and laboratory sections* 



Lecture Master Cards 



OcHTVl v>Oi CLuLll 




1-4 


Department Code 


5-6 


1st Semester Course Number 


7 


1st Semester Footnote Symbol 


8-9 


2nd Semester Course Number 





2nd Semester Footnote Symbol 


11-29 


Course Description 


30-32 


Hours of Credit 


33-36 


Course Number and Section 


37-47 


^ime 


43-62 


Instructor 


63-71 




76 


Special Case Code 


77 


Time C Q de 


78-79 


Department Number 


ao 


first card for each class 



succeeding lecture cards for 
each class 



Laboratory Master Cards 



Card Column Contents 

I- 4 t Department Code 

5-6 1st Semester Course Number 

7 1st Semester Footnote Symbol 

8-9 2nd Semester Course Number 

10 2nd Semester Footnote Symbol 

II- 29 Course Description 
33-36 "Lab" and Section 
37-47 T ine 

43-62 Instructor 

63-71 Room 

76 Special Case Code 

77 T ime Code 
80 "o" 
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CLASS MASTER PILS (Cont.) 



o 

Comment Cards 



Card Column Cont aits 

1-4 Department Code 

5-6 1st Semester Course Number 

7 1st Semestc:; ^otnote Symbol 

8-9 2nd Senester Course Number 

10 . 2nd Semester Footnote Symbol 

12-29 Commentary 

78-79 Department Number 

SO Card Code 



CLASS HEADER FILE 



Card Column 


Contents 


1-11 


x ime 


12-20 


&ocm 


21.-35 


Instructor 


38-41 


Department 


42-43 


Course Number 


44-45 


Section 


47-48 


Department Code 


49-67 


Course Name 


68-69 


Hours (scaled (1,1)) 


77-79 


Number of Places in class 


80 


"o u card code 


STUDENT HEADER FILE 


Card Column 


Contents 


1-9 


Social Security 


10-13 


Sequence Number 


14-35. 


Name 


36 


Sex 


37 


Classification 


80 


f, o H card code 



STUDENT HASTER FILE 
KAME CARD NO. 1 (Gray) 



1-9 Social Security Number 

10-13 Sequence Number 

14-lo .J.u 

36 Se:c 

37-42 Birthdate 

43-57 Birthplace 

53-63 High School (CSKB Code) 

64-67 Date Graduated from High School 

68-78 Blank 

79 Classification at Entrance 

(7 indicates freshman with credit) 

80 "1" punch 



RELATIVES CARD NO. 2 (no Stripe) 

Not in all student packets 
Essentially free form 
Abbreviation followed by name 
Different people separated by commas 



80 "2" punch 

TEST SCORE CARD NO. 3 (green) 

1-9 Social Security Number 

10-13 Sequence Number 

14-17 High School GPA 

18-21 High School GPA (Solids) 

22-30 Rank in High School Class 

(Slash in Column 26 must be there 
for statistical studies) 
31-34 Grade prediction 

35-38 Date of CSEB Scores 

39-41 Verbal Score 

42-43 Percentile for Verbal Score 

44-46 Mathematical Score 

47*48 Percentile for Mathematical Score 

49-52 Date of Achievement Exams 

53-73 Data for three tests 

2 digit Identifier 

3 digit Score 

2 digit percentile 

74-7$ 1 - . ..Matriculation Date 

80 "3" puncli 
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STUDENT RASTER FILE (Con'O 



HOME ADDRESS AT ENTRANCE CARD NO. 4 

1-9 ^ceial Security Number 
10-15 ■ Sequence Number 

14-34 KvrAcr and Street* 

35-47 Ci-y^ 

43-61 Stav.,, Zip (^ountry)* 

62- 67 Continenu*' 

63- 69 Blank 

70-79 Church Preference 
30 punch 

~IIc;::.a addr:..;s at entrance. ' Lc.?t 
adjusted on cards w 

HOME ADDRESS CARD NO. 5 

1-9 Social Security Number 

10-13 Sequence Number 

14-34 Number and -Street -21 "eft 

35-47 City -13 Adjusted 

48-61 State & Zip (Country) -14 on Cards 

62- 67 Continent , -5 

63- 69 State Code C*or geographic distributi 
70-72 Blank 

73-79 Home Phone (No area code) 

80 "5" punch 

LOCAL ADDRESS CARD NO. 6 (Blue) 

1-67 Same Forraat as Kane Address Card 

68-72 . Blank 

73-79 Local Telephone Number 

80 u 6" punah 

PARENTS NAKE CARD NO. 7 (ttose) 

1-9 ^ccial Security Number 

10-13 Sequence Number 

14-30 T itle . % 

31-52 Name ( irst, Initial, Last) 

53-76 • Blank " 

77 M or F (Mother or Father) if last 

name of parent is different from 
last name of student 

73 G if Guardian (only if guardian is 

not father) 

79 "1", "2", or "3" depending on order 
of name and address cards** 

80 "7" punch 
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STUDENT MASTER ICLIa (Cent.) 



•">1V^2N?S ADDRESS CARD NO. 3 (Rosa) 

1-67 Sjirae Format as Heme address card 

60-78 BXar.lc 

?A "I", w 2*% or t7 ^ M depending or order 

of n-^a and addrec cards .* 

1 hcrc are possibly 3 sets of .v7 Si #8 cards* 
They are distinguished by 1 9 2, 3 in column 
79 of the 7 and S card, 

PARENTS INFORKATION CARD IvO . 9 (Salmon) 



1-9 Social Security Number 

10-13 Sequence Number 

14 Father living (y or n) 

15-30 Father's Occupation 

31-46 Name of ^ira 

47 Mother Living (y or n) 

48-63 Mother'* s Occupation 

64-79 i^aine of Firm 

80 "9" punch 

PERSONAL DATA CARD NO. 10 (Brown) 

1-9 Social Security Number 

10-13 Sequence Number ••* 

14 Marital Status (m or S) 

15-29 Name of Spouse 

30 Car Operation (Y or N) 

31-32 . Fraternity 

33 Veteran Status 

34-37 Degree of Self-support 

38-46 Re-entry dates (3 digits each) 

47-52 Counselor Code (2 3-digit codes) 

53-67 Counselor (s) 

68-71 Major (2-digit code) 
Two majors possible 

72-79 Other Colleges attended (CEEB codes) 

80 V unch 



** This card appears only if parents t 
address is different from home address* 



^7 



1-9 

10-13 
14 



15-21 



22-25 

26-29 

30-35 
36-37 
33-40 
41-42 

43 
44 

45-48 
49-51 
52-54 
55-57 
58*60 
61-63 
64-66 
67-75 



76-79 
80 



STUDENT MASTER FILE (Cent.) 



GARD-KO. 11 (Purple) 

^ocial Security Number 

Sequence Number 

Grccv itiovi on Schedule 

1 - yes 3 - Late 

2 - early blank- Did not graduate 
■ Ratk in Glass 

(Punch rank in class at er of each 
semester f but by-pass incomplete 
and deferred grades) 

Cumulative GPA to date 

(By-pass incomplete and deferred grades)* 

Cumulative Grade Points 

(By-pass incomplete and deferred grades) 

Graduate Record Examination Date 

Test Identifier 

Score 

Percentile 

Classification 

Number of Semesters in Residence 
Hours to Date (Including Credit) 
Non^graded Hours to Date 
Cumulative Hours of F 
Previous Semester Graded Hours 
Previous Semester Grade Points 
Transfer Hours 
Hours of F in last semester 
Scholastic Gases 

1 - Dean f s List 

2 - Undergraduate Honors 

3 - Scholastic Probation 

4 - Dropped 
Blank 

punch 
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SOCIAL SECURITY NO. 


SEQ. NO. 


STUDENT NAME 


i 




CEPT. 


NO. 


>EC 


3 




COURSE NAME 


HRS 






points 


i 


YtAB 
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WHITMAN COLLEGE 

TO THE STUDENT: 

FILL IN THIS PORTION OF THE CARD. 



UAST MAIM 



rmtr nam* 



tMOOLt MA9#* 



CLASSIFICATION : FR SO JR SR GR SP. 



DO NOT FOLD, TEAR, OR DAMAGE THIS CARD. 



** 47 414190 31 $2333435313790111011121)141351 17 St» 10 71 7271 74 7$ 71 71 71 TIM 

Qa 



TO THE INSTRUCTOR: 

TO REPORT FINAL GRADE, USE 
SHARP PENCIL TO PUNCH DESIRED 
GRADE. IF AN ERROR IS MADE, PUNCH 
THE CORRECT GRADE AND CIRCLE 
CORRECTION. 



□ b 

□ c 

□ o 
Of 
Oct 
Qwf 

. Owp 

INSTRUCTOR'* INITIALS ^ IP 

IF 



r 



SOaAi. SECURITY MO. SEQ. HO. 



STUDENT NAME 



iU\ DEPT.] NoJs EC.T 



, ^ w j t . . . COURSE NAME ■ • 

1 2 1 4 5 t 7 t I 10 11 12 11 14 g H T7 WH2>2122n2425M 27 2l2f30313233M 35 3t3PM39404142434 45 41 47 41415051 S233S4 MW37H59I0H 12 03 149316 17 10 49 70 71 72 
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WHITMAN COLLEGE 

TO THE STUDENT: 

FILL IN THIS PORTION OF THE CARD 



3 



LAST NAME 



Ml DO IK MAMI 



CLASSIFICATION : FR„ SO JR _SR GR !_SP. 



DESIGNATE LECTURE SECTION, IF ANY . 



(8MJ74A70 



73 74 75 78 77 71 71 10 



>- Q 

o 2 

CO I— 
< U 



f1 



o 



SEX VETERAN 



STUDSMT'S NAME, home address * TELEPHONE 



FRATERNITY 



COURSE 



NO. 



AUTOMOBILE 



STUDENT'S LOCAL ADDRESS * TELEPHONE 



NAME OF SPOUSE 



PARENT S OH GUARDIAN'S NAME * AODRESS 



CR. 



FATHER LIVING 



CLASSIFICATION 



ROOM 



HIGH SCHOOL 



INSTRUCTOR 



MOTHER LIVING 



DATE GRADUATED 



OTHER COLLEGES ATTENDED & DATES 



BIRTH DATE 



CHURCH PREFERENCE 



BIRTHPLACE 



SELF-SUPPORT. 



RELATIVES. INCLUDING SPOUSE. IF ANY. WHO ATTENDED 
WHITMAN AND THEIR RELATIONSHIP TO YOU. 



CAMPUS ACTIVITIES A HONORS IN COLLEGE 



FATHER'S OCCUPATION 



NAME OF FIRM 



MOTHER'S OCCUPATION 



NAME OF FIRM 



e 



e 
• 



©. 



JONES JOHN 





HIST 


27 


ENGLISH HISTORY 


3.0 


B 


9.0 


WHITMAN COLLEGE 


HIST 


55 


HIST OF FAR EAST 


3.0 


B 


9.0 


WALLA WALLA, WASH. 


HIST 


71 


HISTORY OF RUSSIA 


3.0 


A 


12.0 


SEMESTER GRADE REPORT 


ART 


11 


ART APPRECIATION 


3.0 


B 


9.0 




GEOL 


21 


HISTORICAL GEOLOGY 


4.0 


A 


16.0 



MR AND MRS ROBERT JONES 
4336 PENSHURST CT 
SACRAMENTO 
CALIF 95825 



SEM GRADED HOURS 16.0 

SEM GRADE PTS 55.0 

SEM GPA 3.437 

TOT WHITMAN GR HRS 094.0 

CUM GRADE PTS 286.0 

CUM GPA 3.042 

TOTAL HRS EARNED 128.0 



DOE JOE 





SOC 


64 


POP - HUMAN ECOLOGY 


3.0 


B 


9.0 


WHITMAN COLLEGE 


HIST 


51A 


AMER REVOL - CONST 


3.0 


C 


6.0 


WALLA WALLA, WASH. 


CHEM 


91 


PHYSICAL CHEMISTRY 


4.0 


B 


12.0 


SEMESTER GRADE REPORT 


BIOL 


51 


VERT EMBRYOLOGY 


5.0 


B 


15.0 




CHEM 


93 


SEMINAR 


1.0 


C 


2.0 



MR AND MRS JAMES DOE 
7004 SE 20TH 
MERCER ISL 
WASH98046 



SEM GRADED HOURS 16.0 

SEM GRADE PTS 44.0 
SEM GPA 2.750 

TOT WHITMAN GR HRS 122.0 

CUM GRADE PTS 320.0 
CUM GPA • 2.622 

TOTAL HRS EARNED 126.0 



CRS SMI SM2 NAME - PREREQUISITE CR SEC 

POLS 84 SEE HIST 84 

POLS 87 88 SPECIAL PROBLEMS 1-3 
PRE CONSENT 

POLS 89 90 THESIS— SEMINAR 1-3 
PRE CONSENT 

POLS 98 HONORS THESIS 3 

PRE CANDIDATE FOR 
MAJOR HONORS 



TIME 



7 8 T 



INSTRUCTOR 



FLUNO 



FLUNO 



FLUNO 



ROOM 



R315 



R103 



P315 



PS YC 11 PRINCIPLES 



PSYC 11 PRINCIPLES 



PSYC 12 EXP ANAL BEHAVIOR 

PRE 1 1 



PSYCHOLOGY 

3 A 
B 

3 A 
B 



10 MW 11 F 
9 MTTH 

8 M 9 WF 

9 MTTH 



9 MTTH 

LAB W 1 2 MW 

LAB X 3 4 MW 

LAB Y 1 2 TTH 

LAB Z 3 4 TTH 



EACKER 
FOGARTY 

EACKER 
FOGARTY 

MEYER 
MEYER 
MEYER 
MEYER 
MEYER 



S120' 
M31 1 

S120 
M31 1 

B408 



PSYC 37 SEE ECON 37 

PSYC 43 SOCIAL PSYCHOLOGY 

PRE 1 1 OR SOC 3* 
SAME AS SOC 43 

PSYC 53 ABNORMAL PSYCH 

PRE il - JR 

PSYC 54 EDUCATIONAL PSYCH 

PRE 11 - JR. FOR 
TEACHER CANDIDATES 



A 
B 

A 
8 



11 MTTH 



1 MWF 
10 TTHF 

1 MWF 
10 TTHF 



CHERTOK 



FOGARTY 
FOGARTY 

FOGARTY 
FOGARTY 



B408 



B208 
B308 

B208 
B308 



PSYC 63 CHILD PSYCHOLOGY 

PRE 1 1 



11 MTTH 



MEYER 



B208 



PSYC 67 TESTS - MEASUREMNTS 

PRE 37 



1 2 TTH 



EACKER 



B302 



PSYC 71 



PERSONALITY THEOR 
PRE CONSENT 



3 MWF 



MEYER 



B302 



PSYC 80 

PSYC 81 

PSYC 

PSYC 83 

PSYC 96 
PSYC 



THEOR OF LEARNING 
PRE CONSENT 

SEMINAR MOTIVATION 
PRE CONSENT 

.82 SEMINAR 

PRE CONSENT 

SEM HIST-SYST PSYCH 

96 THESIS 

98 HONORS THESIS 

PRE CANDIDATE FOR 
FOR MAJOR HONORS 



2 MWF 



MEYER 



2/30 3 TTH EACKER 



1 MWF 
1 MWF 

♦ * 



EACKER 

EACKER 

FOGARTY 

FOGARTY 



B302 

B302 

B302 

B302 
R306 
R306 



15 /*?2-~ 



COOPER SUSAN JEAN FR ANDERSON HAUL SEC C 

COPELANO CHARLES M JR 1005 ISAACS 

COPEL AND THOMAS A FR LYMAN HOUSE SEC E 

CORNUE MARTHA JANE SR PRENTISS HALL SEC B 

COTTLE JAMES SIDNEY SR 715 ESTRELLA 

COURTN AGE JOHN CREEGAN FR JEWETT HALL 403 

COWEN MARILYN JEAN SR PRENTISS HALL SEC A 

CRAIG PATRICK EDGAR FR LYMAN HOUSE SEC E 

CRAMER SYDNEY ANN SO PRENTISS HALL SEC E 

CRAVEN DAVID ROGER JR 949 ISAACS 

CRAVEN HENRY TRUXTON SR 949 ISAACS 

CRAWFORD KAREN HALEEN PR SAME AS HOME ADDRESS 

CROSBY JOHN LLOYD FR JEWETT HALL 402 

CROWE LL ROBERT IKUA ' FR JEWETT HALL 210 

CRUTCHER DONNA JEANE SO SAME AS HOME ADDRESS 

CUDNEY CHARLES HENRY FR JEWETT HALL 406 

CULHAM WILLIAM M M JR 809 VALENCIA 

CULLEY JAMES ROBERT JR 949 ISAACS 

CUMMINS GALE CATHERINE M SP SAME AS HOME ADDRESS 

CUNNINGHAM PAUL ALAN FR JEWETT HALL 325 

CUNNINGHAM TIMOTHY J SO 1005 ISAACS 

CURFMAN GEORGE HARDIN SO 949 ISAACS 

CURRENT NANCY FR PRENTISS HALL SEC B 

CUSACK CHARLES ROOS JR SO 715 ESTRELLA 

CZINGER JOHN LOUIS III FR LYMAN HOUSE SEC B 

DAI NES HOLUS LOUISE SR PRENTISS HALL SEC G 

DALE FORREST PAUL SO 31 WHITMAN 

DALLA IRE MARSHA ANN FR ANDER SON HALL SEC B 

DALY KATHLEEN DIANA FR ANDERSON HALL SEC F 

DANCE DOUGLAS ELORED FR JEWETT HALL 428 

DANIELSON SUE JEAN SO PRENTISS HALL SEC G 

DANISHEK STEPHAN D SO 716 BOYER 

DAUGHERTY EILEEN P SO PRENTISS HALL SEC C 

DAUGHERTY PATRICIA A JR PRENTISS HALL SEC F 

DAUGHTERS EUGENE A JR M JR 19 E BIRCH 

DAVIDSON CHRISTOPHER SR 429 LINCOLN 

DAVIES HEATHER MARIE FR ANDERSON HALL SEC F 

DAVIS BARBARA JOAN SP ANDERSON HALL SEC C 

DAVIS BRUCE VANARSDALE FR JEWETT HALL 416 

DAVIS JOHN M JR SR 715 ESTRELLA 

DAVIS PAUL FREDERIC FR JEWETT HALL 312 

DAWSON ELIZABETH ANN JR PRENTISS HALL SEC B 

DAWSON MARY SR PRENTISS HALL SEC G 

DAY MICHAEL JUDO SO 716 BOYER 

DAYTON KENNETH ALLAN JR 925 ISAACS 

DEAN WILLIAM WENDELL SO 715 ESTRELLA 

DEBUTTS BOBETTE FR ANDERSON HALL SEC E 

DECHERT HEDY SUSAN SR PRENTISS HALL SEC G 

DECOLAN RALPH RODNEY SO LYMAN HOUSE SEC A 

DEGRASSE MICHAEL E FR JEWETT HALL 319 

DELZELL DAVID ATHOL FR JEWETT HALL 202 

DEMEULES JAMES HEAD JR 1005 ISAACS 

DEMMY STUART DOUGLAS FR JEWETT HALL 430 

DENNEN RODGER S TAYLOR FR JEWETT HALL 109 

DENNIS WILLIAM G JR 716 BOYER 

DERMOND DONNA CAROLYN SO PRENTISS HALL SEC G 

DESHLER DONALD DAVID SO 715 ESTRELLA 

OETWUER ANNEKE JAN FR ANDERSON HALL SEC E 

DIBBLE WAYNE ROBERT SO LYMAN HALL SEC B 

DICKEL CHARLES TIMOTHY SO 925 ISAACS. 

DION JANE SHAW SO PRENTISS HALL SEC C 

DITMARS FRANK R JR FR JEWETT HALL 401 

DOCHEZ MARC WELCH FR JEWETT HALL 217 

DORAN GARRETT JOHN FR JEWETT HALL 106 

DOUGLAS DAVID MICHAEL SR 401 E WHITMAN 

DREHER ANNA LESLIE JR PRENTISS HALL SEC G 

DRUMMOND HOWARD B JR 1005 ISAACS 

DUNCAN DIANA LYNN SO PRENTISS HALL SEC B 

DUPREE JOHN ROBERT FR 1005 ISAACS 

DUSENBERY THOMAS J FR JEWETT HALL 313 

DWIGGINS MICHAEL BELL SO 949 ISAACS 

DWIGGINS STEVEN FRANK SO 949 ISAACS 

DYE DAVID GARY SO 715 ESTRELLA 

EDDY JOHN WHITTEMORE FR JEWETT HALL 409 

EAGLESON JEANNE ANNE SO PRENTISS HALL SEC C 

EASTON ANN LOUISE FR ANOERSON HALL SEC E 

EBY BEN WILDER SP SITTNER HALL COLL PL 

EDWARDS GEORGE MICHAEL SR 16 S PARK 4 

EICHELBERGER DONALD N FR JEWETT HALL 217 

ELDERK IN RICHARD H JR 949 ISSACS 

ELLIS WILLIAM L SO 1005 ISAACS 

ELLS WILLIAM ARTHUR SO 715 ESTRELLA 

ELMORE RICHARD F SR 135 S PARK 

EMEL JOHN ERIC SR 633 MILITARY RD 

EMPEY SUSAN LEE SO PRENTISS HALt SEC B 

ENGLAND LAUREL AWN FR ANDERSON HALL SEC E 

ERICKSON KENT LYNN SO LYMAN HOUSE SEC C 

ERNST CHARLES JOHN SO 949 ISAACS 

ESSER JILL SR PRENTISS HALL SEC A 

ESTRIN ROBERT NORMAN SO 1005 ISAACS 

EVANS CAROLYN JANE SR PRENTISS HALL SEC A 

EVANS MICHAEL JOHN JR LYMAN MOUSE SEC A 

EVANS SUSAN ETHEL SO PRENTISS HALL SEC B 



JA5-9898 157 TROY 
JA5-3780 944 BRYANT 
JA5-9887 5460 BON ITA PL 
JA5-8972 140 KING 
JA5-6542 12300 MELODY LN 
JA5-9839 2213 FORAKER DR 
JA5-3240 1876 OAK KNOLL CT 
JA5-9887 728 E WHITMAN 
JA5-8642 608 OAK 
JA5-5653 3039 DIXON PL 
JA5-5653 KEELER RD 
JA9-2240 2037 CRAWFORD DR 
JA5-9839 9915 S TACOMA WAY 
JA5-9979 BOX 598 
JA5-8763 1920 LOUBECK 
JA5-9839 2009 S 10TH 
JA5-4187 113 SE 10TH 
JA5-5653 BOX 727 
JA5-8667 610 TAYLOR 
JA5-9878 1914 MUKILTEO 
JA5-3780 3353 COTTAGE WAY 95 
JA5-5653 20 ELM 
JA5-8972 9724 MERCERWOOD DR 
JA5-6542 4111 LOS COCHES WY 
JA5-9988 3142 KEMPTON DR 
JA5-8792 3509 WOOD ACRES DR 
JA5-3426 1104 WASHINGTON 
JA5-9896 1401 PERKINS 
JA5-9946 9281 56TH S 
JA5-9788 110 EL MONTE 
JA5-8792 1424 E BROADWAY 
JA5-3822 4714 NE 36TH 
JA5-8232 3026 NE 92ND 
JA5-8592 846 GREENWICH PL 
JA5-1227 1556 30TH WEST 
JA9-5336 637 WARM SPRINGS 
JA5-9946 1818 HOLBROOK 
JA5-9898 16751 17TH NW 
JA5-9839 7662 SE 22ND 
JA5-6542 7662 SE 22ND 
JA5-9962 2722 S LINCOLN 
JA5-8972 RT 1 BOX 42 
JA5-8972 RT 1 BOX 42 
JA5-3822 718 N GEE 
JA5-8443 1318 LAKESIDE S 
JA5-6542 1617 SPRUCE 
JA5-9846 3243 EVERGREEN PT RD 
JA5-8792 110 W OTIS RD 
JA5-9847 3014 F 
JA5-9962 725 YAKIMA 
JA5-9978 65115 VIEWS RD 
JA5-3780 6815 44TH PL NE 
JA5-9783 2904 NE 46TH 
JA5-9949 1139 23 E 
JA5-3882 5167 S CRESTON 
JA5-8792 805 N FOOTE 
JA5-6542 615 W MERCURY 
JA5-9846 10202 SE 28TH 
JA5-9988 BOX 16 
JA5-8443 850 NW POWHATAN TERR 
JA5-8232 818 STUART 
JA5-9839 505 PALOS VERDES DR W 
JA5-9979 3959 LA DONNA 
JA5-9949 912 6TH 

412 BIRCH 
JA5-8792 2305 DALLAS 
JA5-3780 141 NORTH E 
JA5-8972 2154 MCMILLAN 
JA5-3780 306 WEST 
JA5-9962 12122 S E 21ST 
JA5-5653 3 STANHOPE PLACE 
JA5-5653 3 STANHOPE PLACE 
JA5-6542 1261 C SW 
JA5-9839 4515 W RUFFNER 
JA5-8232 625 S PALOUSE 
JA5-9846 7704 OXON HILL 

10312 W COURT 

1135 ALVARADO 
JA5-9979 522 32NO S 
JA5-5653 2537 STATE 
JA5-3780 E9404 MAIN 
JA5-6542 BOX 164 

1706 N EASTMONT 
JA5-8443 7321 172 SW 
JA5-3582 1009 W 21ST 
JA5-9846 2504 BARGE 
JA5-9987 616 S 117TH 
JA5-5653 1055 LEMON 
JA5-3240 125 N PERCIVAL 
JA5-3780 833 S CITRUS 
JA5-3240 501 LAKEV1EW BLVO 
JA5-9847 1034 A 48TH 
JA5-B972 4774 NE 180TH 



JUNEAU 


ALAS 99801 


CHICO 


CALIF 


TUCSON 


ARIZONA 


WALLACE 


IDAHO 


LOS ALTOS 


CALIF 


ANCHORAGE 


ALAS 99503 


LAKE OSWEGO 


ORE 


WALLA WALLA 


WASH 99362 


SANDPOI NT 


IDAHO 83864 


PALO ALTO 


CALIF 


BR I DGEWATER 


CONN 


WALLA WALLA 


WASH 99362 


TACOMA 


WASH 


KAUAI WAIMEA 


HAWAII 96796 


WALLA WALLA 


WASH 


LAS VEGAS 


NEV 89105 


PENDLETON 


ORE 


WEED 


CALIF 


WALLA WALLA 


WASH 99362 


EVERETT 


WASH 98202 


SACRAMENTO 


CALIF 95825 


DENVER 


COLO 80220 


MERCER ISL 


WASH 98040 


SACRAMENTO 


CALIF 95825 


LOS ALAMITOS 


CALIF 


BOISE 


IDAHO 83702 


TOPPENISH 


WASH 


RICHLAND 


WASH 99352 


SEATTLE 


WASH 98118 


CONCORD 


CALIF 94521 


MOUNT VERNON 


WASH 98273 


SEATTLE 


WASH 98105 


SEATTLE 


WASH 98115 


PALO ALTO 


CALIF 


SEATTLE 


WASH 98199 


BOISE 


IDAHO 


EVERETT 


WASH 98201 


SEATTLE 


WASH 98177 


MERCER ISL 


WASH 98040 


MERCER ISLAND 


WASH 


SPOKANE 


WASH 99203 


JOSEPH 


ORE 97846 


JOSEPH 


ORE 97846 


TACOMA 


WASH 


SEATTLE 


WASH 98144 


S PASADENA 


CALIF 


BELLE VUE . 


WASH 98004 


BARR I NGTON 


ILL 60010 


VANCOUVER 


WASH 98663 


WENATCHEE. 


WASH 99801 


TACOMA 


WASH 98407 


SEATTLE 


WASH 98115 


PORTLAND 


ORE 97213 


SEATTLE 


WASH 98102 


SEATTLE 


WASH 98178 


OLYMPIA 


WASH 


BUTTE 


MONT 59701 


BELLEVUE 


WASH 98004 


GRUVr LAND 


MASS 


nnh T I Akin 

PORTLAND 


ORE 97210 


HELENA 


MONT 


nil nc wcnncc 
PALUb VtKUtJ> 


LALI r 9027D 


PALO ALTO 


CALIF 94306 


SANTA MONICA 


CALIF 94103 


RICHLAND 


WASH 99352 


RICHLAND 


WASH 


MISSOULA 


MONT 59801 


EUGENE 


OREGON 


WAITSBURG 


WASH 


BELLEVUE 


WASH 98004 


LONDON W2 


ENGLAND 


LONDON W2 


' ENGLAND 


EPHRATA 


WASH 


SEATTLE 


WASH 98199 


WALLA WALLA 


WASH 99362 


OXON HILL 


MD "20021 


PASCO 


WASH 


WALLA WALLA 


WASH 99362 


SEATTLE 


WASH 98144 


BUTTE 


MONT 59701 


SPOKANE 


WASH 


CUTTEN 


CALIF 


E WENATCHEE 


WASH 


EDMONDS 


WASH 


SPOKANE 


WASH 


YAKIMA 


WASH 98902 


TACOMA 


WASH98444 


MENLO PARK 


CALIF 


OLYMPIA 


WASH 


LOS ANGELES 


CALIF 90036 


SANOPOINT 


IDAHO 83864 


LOS ALAMOS 


NEW MEX 87544 


SEATTLE 


WASH 
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WHITMAN COLLEGE BUDGET FORECASTING PROGRAM 




Each year Whitman College prepares ten year forecasts for the operating 
budget. In the past this has been done by hand computation. These cal- 
culations of income and expense are quite simple, but take considerable 
time* Only a few variations of income and expense could be calculated 
because of the time involved in making calculations ♦ 

This problem of calculation time was solved by xrciting a computer pro- 
gram to do the computation of income and expense* The budget forecasting 
program reads the input data, performs the calculations and prints the 
results in appro:cimat..'.y six minutes # This short time makes it feasible 
to recompute income and expense for any change of input data. Thus input 
data can be varied until a satisfactory solution can be obtained. 

In writing a program of this type there are several planning decisions 
to be made. An example is the consideration of endowment growth. Endow- 
ment growth can be projected as a percentage growth or an estimate can be 
made of the endowment principal for each of the ten years. It was decided 
to use the latter method for the budget forecasting program. In contrast, 
the program calculates scholarship costs based on a single piece of input 
data called the scholarship ratio. The scholarship cost ratio is multiplied 
by tuition income for each of the ten years to yield scholarship costs. 
The choice of method for any item of income or expense depends on the par- 
ticular needs of a college* One can modify any part of this program to meet 
particular needs. 



o 



The following items are input data to the program: 



Income Data 



Fortran Variable 



Number of students 

Number of Resident Students 

Number who board only 

Research expenditures 

Tuition and Fees 

Room Charges 

Board Charges 

Endowment Principal 

Gift Income 

Research overhead ratio 
End owment yield 




There are 
eleven entries 



for each of 
these items f 
one for each 
year of the 
forecast. 



G 




Cost Data 



Fortran Variable 



Faculty Salary Increment Ratio 

Special costs 

Number of Faculty 

First year average faculty salary 

Cost of living factor (non-salary items) 

Administrative Salary Increment 

Ratio of scholarship costs to tuition fees 

Contingency 

Faculty Fringe Benefits Ratio 
Administrative Fringe Benefits Ratio 




11 entries 



F 

H 

HH 
SF 
C0 




11 entries 




-2- 



All costs other than faculty salaries, faculty fringe benefits, scholar- 
ship costs, special costs, and contingency are treated as budget categories . 
Each budget category is split into a salary and a nonsalary part # Further, 
a budget category is considered as independent of enrollment growth or 
directly proportional to enrollment -growth. There is • cvision for adding 
a jump in either the salary or non-salary part of any bucget for any year. 
The Input data for the budget categories are the 03 previously mentioned 
and the salary and non-*salary parts for each budg^u category for the first 
year* 

The nonsalary variables are X(J,X) and the salary variables are Y(J,I) 
■/here J is the budget category and varies from 1 to 12 and I is the year 
and varies from 1 to 11 . The number of budget categories which are inde- 
pendent of enrollment growth is arbitrary and can be changed. Presently, * 
th :?irst five budget categories are independent of enrollment and the last 
Sv^u are directly proportional to enrollment. 

The program calculates for each of 11 years: 



research income 

tuition and fees income 

room income 

board income 

endowment income 

total income 

student faculty ratio 

average faculty salaries 

total faculty salaries 

faculty fringe benefits 

scholarship funds 

budgets for all budget categories 

total expenditures 

surplus 

surplus/student 
expense/student 



All input data and calculated results are printed. There is one 
piece of data per card with a format of F 10.0 • *he order of input can be 
determined by studying the Fortran listings . 




C WHITMAN COLLEGE BUDGET FORECAST PROGRAM 

C VERSION 2 

C 

C FIRST L INK OF THREE LINK PROGRAM 

C 

DIMENSION A(ll) , B ( 1 1 ) , C i 1 1 ) , D ( 1 1 ) , E ( 1 1 ) , F ( 1 1 >,Q(11),K{ 11) 

DIMENSION S(ll) ,T(11),U(11) ,V( 11),W(11) ,X(12,11) ,Y(12,11),Z( 11) 

DIMENSION P(ll) 

DIMENSION AA(ll), AAA(ll) 

DIMENSION FS( 11 ) 

DIMENSION FFB(ll), AFB(ll) 

COMMON A, B,C,D,E,F, Q,R,S,T,U,V, W,X,Y,Z,P 

COMMON L » K 

COMMON H, HH, SF, CO 

COMMON FFB, AFB 

COMMON FS 
1 FORMAT (F10.0 ) 

3 FORMAT { 1H1,23X,8HRESIDENT» IX, 1 2HS TUDENTS WHO) 

4 FORMAT ( 13X, 8HSTUDENTS, 3X, 8HSTUDENTS » 3X, 10HBOARD ONLY » 

1 3X, 8HRESEARCH, 2X, 

2 12HTUITION+FEES, 6X , 4HR00M, 7X, 5HBOARD , 2X, 11HGIFT INCOME, 

3 2X, 9HEND0WMENT) " " 

5 FORMAT ( 5H0YEAR, 13, 6( 7X, FIO.O)) 

- 6 FORMAT ( 24H0RESEARCH OVERHEAD RATIO, F10.2/ 

1 16H ENDOWMENT YIELD, 1 OX , F10.4) 

7 FORMAT ( 1H0,11X, 1 5HRESEARCH INCOME, 3X, 19HTUITI0N+FEES INCOME, 
1 2X, 11HROOM INCOME, 

r ; 2 5X, 1 2HBOARD INCOME, 3X, 16HEND0WMENT INCOME, 2X,12HTOTAL INCOME) 

8 FORMAT! 1H1, 8X, 23HFACULTY SALARY INCREASE, 5X, 14HN0. OF FACULTY 
..- i,5X, 21HSTUDENT FACULTY RATIO) . --' 

9 FORMAT ( 5HOYEAR, 13, 7X , F10.3, 12X, FlO.l, 12X, F10.2) 

10 FORMAT ( / 22 H COST OF LIVING FACTO 
1R, 5X, F10. 3/ 27H ADMINISTRATIVE SACARY INC.,- .... F10.3/ 

218H SCHOLARSHIP RATIO, 9X, F10.3) 

11 FORMAT (5H0YEAR, 13, 9( 2X, FIO.O)) - 

12 FORMAT (5H0YEAR, 13, 5X, F10.3,13X, F10.3) ;. 

13 FORMAT (1H0, 16X, 7HFACULTY, 15X, 14HADM IN ISTRAT ION / 
1 17X, 15HFR INGE BENEFITS, 7X, 15HFRINGE BENEFITS / 

- . 2 17X, 5HRAT 10 ,. 17X, 5HRAT 1 ) - - ,.- .. . _ 

DO 101 I = l\t 

XU,I) =0 " 

101 Y( J,J ) = ™- .. ^ .. rr . : . ... .. . . . ..... rr _, : , 

READ 1,A(I),AA(I ),AAA(I ), B ( I ) , C { I ) ,D ( I ) , E I I ) , F ( I ) , P ( I ) 
100 PRINT 11,1 » All ),AA(I ),AAA<I),B( I), C(I), DU), Ed) , P(I), F(I) 

_ . READ 1 »G , i_ 

:^^:PR-i : f#£^Spi^ 

DO 21 1 = 1,L %~ " """"" ~33-~ - -.3"- . ".- " ' 

Q( I )*&{ 1)*G 



. ■ ■ \ 

R (I )=C(I )*A(I ) f\ 

s ( i ) = u( i )*aa( i ) 

T(I) = E(I)*AAA(I) 
U(I)=F(I)*0 

V (I )=Q(I )+R(I )+S(I )+T{ I )+U( I }+P( I ) 
21 PRINT 5,1, Q(I),R(I),S(I),T(I),U(I),V{I) 
DO 32 1=1, L 
READ 1,D(I) 

32 READ 1,E (I ) 
READ 1,F( 1) 

READ 1, (P( I ) ,I = 1,L) 

READ 1, H,HH,SF - 

READ 1, CO 

READ 1,(FFBU), 1=1, L) 

READ 1 , ( AF B { I ) , 1 = 1, L) 

DO 30 1= 1, L 
30 FS(I) = A ( I ) /P ( I ) - 

PRINT 8 . ■ .. 

DO 33 I = 1,1 

33 PRINT 9, I , D( I ) , P (I ) , FS( I ) 
PRINT 10, H, HH, SF 

PRINT 13 ' W . : : 

DO 600 1=1 ,'L 

600 PRINT 12, I , FFB { I ) , AFB { I ) -. ^ - 

CALL LINK (F0R2) 
END 

■ o 



c 



c 



c 



WHITMAN 
VERSION 



COLLEGE 
2 



BUDGET FORECAST PROGRAM 



SECOND LINK OF THREE LINK PROGRAM 



o 



1 

2 

104 
105 

34 

35 
201 

300 
301 



200 
31 

41 



51 



42 



44 



50 
60 



45 
43 



DIMENSION A( 1 1 ) , B ( 1 1 ) , C ( 1 1 ) , D I 1 . ) ,E( 
DIMENSION S(11),T(11),U(11),V(11),W( 
DIMENSION P(ll) 
DIMENSION FFB(ll), AFB(ll) 
DIMENSION FS(ll) n 

COMMON A,B,C,D,E,F,Q,R,S,T,U, V,W,X,Y 

COMMON L t K 

COMMON H, HH, SF t CO 
COMMON FFB, AFB 

COMMON FS 

FORMAT (F10.0) 

FORMAT ( 12, 12, F10.0, F10.0 ) 

FORMAT ( 1H1 , 9X , 24HAVERAGE FACULTY 

FORMAT ( 5H0YEAR, 13, 10X, F10.0) 

DO 34 J = 1,K 

READ 1,X< J,l) 

DO 35 J=1,K 

READ 1,Y ( J, 1) 

READ 2,M,N,AZ1,AZ2 

IF (M-K) 300,300,200 

IF (N-L ) 301,301,200 

X(M,N) = AZ1 : . ... - ,-. 

Y(M,N)=AZ2 ." ■ 
GO TO 201 

DO 31 1= 1, L _ . ■■: . 

PU ) = FS(I) 
DO 41 1=2,1 

F(I )=F(I-1)*D(I } _ r . . 

PRINT 104 - - 

DO 51 1=1, L . • - I - 

PRINT 105, I, F(H - _-_ -I- ^ 

DO 42 1=1, L 

F(I)=A(I)/ P(I)*F(I ) 

Q { I ) = F ( I )*FFB{ I * . _- - - _. r 3 

DO 49 J«1,K - ~ \ > T^-^.^^-ff -r-- - 

BBB=A(1) -_ 

; DO /.44^^^<S=^i^^^^i^^^^^ 
X{J,I )«X<J,I-1)*H + X(J,I ) 
Y(J,I )=Y(J,I-1)*HH + Y(J,I ) 
Y(J,I )- = -^fUil ) +Y(J,U -*- AFB (I) - 
X{J,1) * *<J,1) +- Y(J»D- -^> : 

■■ o o 43 m^m^Bwmmm;;m^^^g&^ 

X ( J y I ) *X ( J , I ) +Y ( J t I ) ™^ Zr^^Jc—^H-^ 

IF <X(J,I-1)> 50,50,45 
IF (X(J,1)> 43,43,60 

GO ■'^■■i^^^^W^&^^^^^k^^U 
Y { J , I )=X ( J, I )*A( I J /BBai- -^k=^<ZWJ=t 
CONTINUE 

r JOf 



11 
11 



) »F 
),X 



(11) 
(12, 



,Q( 
11 ) 



11) ,R{ 
,Y(12, 



11) 
11) 



•Z ( 11) 



,Z,P 



SALARIES) 




49 CONTINUE 

C AT THIS POINT THE X ARRAY CONTAINS NONSALARY+SAL ARY WITH NO GROWTH FACTU 

C AT THIS POINT THE Y ARRAY CONTAINS NONSALARY+SAL ARY TIMES GROWTH FACTOR 

C HERE GROWTH FACTOR MEANS ENROLLMENT 

DO 46 1=1, L 

46 R (I )=R(I )*SF 
DU 47 1=1, L 

S ( I ) =X ( 1 , 1 ) +X ( 2 , 1 ) +X ( 3 » I ) + X ( 4 , 1 ) +X ( 5 , I ) 

47 T( I )=Y(6, I )+Y(7,I )+Y(8,I )+Y(9, I ) + Y( 10, I )+Y( 11, I )+Y( 12,1 ) 
DO 48 1=1, L 

U i L }=B( I )+F ( I )+Q(I)+R (I )+S( I )+T( I )+E( I )+C0 

48 W (I J=V (I )-U( I ) : 
CALL LINK (F0R3) 

END , 





o 







C WHITMAN COLLEGE BUDGET FORECAST PROGRAM 

C VERSION 2 

C 

C THIRD LINK OF THREE LINK PROGRAM . 

DIMENSION A(11),B(11),C(11),D(11),E(11),F{11),G(11),R( 11) 
DIMENSION SU1),T<11),U(11),V( 1 1 ) , W ( 1 1 ) , X ( 12 , 1 1 ) , Y ( 12 , 1 1 ) , Z { 11) 

DIMENSION P( 11 ) 

DIMENSION FFB(ll), AFB(ll) 

DIMENSION ZZ( 11) 

COMMON A,B,C,D,E,F,Q,R,S,T,U,V,W,X,Y,Z,P 

COMMON L t K ...=~ . .-• . . . . -^~r:; ■ • • ■ ' 

COMMON H, HH, SF , CO 
106 FORMAT ( 1H0 , 10X , 1 3HRESEARCH EXP., 3X,22HTOTAL FACULTY SALARIES,' 3X, 

1 18HFACULTY RETIREMENT, 3X, 17HSCH0LARSHIP FUNDS, -5X, 

2 13HSPECIAL COSTS) ~" " ' 

107, FORMAT ( 5H0YEAR, 13, 4X , F10.0, 7X, F10.0,12X, F10.0,12X, F10.0, 
112X, F10.0) - ..'---I . -- - 

108 FORMAT { 1 2H0C0NT I NGENCY , 5X, F10.0) 

109 FORMAT ( 1H1, 18X, 51HBUDGET CATEGORIES NOT SUBJECT TO ENROLLMENT 
1GR0WTH/ 15X, 4HN0.1 f 10X, 4HN0.2,.10X, 4HN0.3, 10X, 4HN0.4, 

2 10X, 4HN0.5) " - . - 

110 FORMAT ( 5H0YEAR, 13, 3X , F10.0, 4X, F10.0, 4X, F10.0, 4X, F10.0, 
1 4X, F10.0) ~>. ' - 

111 FORMAT ( 1H0, 25X, 46HBUDGET CATEGORIES SUBJECT TO ENROLLMENT GROWT 
1H/ 15X, 4HN0.6, 10X, 4HN0.7, 10X, 4HN0.8, 10X, 4HN0.9, 9X , 5HN0.10 
1, 9X, 5HN0.11, 9X,.. 5HN0.12) :._ 

112 FORMAT (' 5H OY E AR ,~ T3 , 3X, F10.0, 4X, F10.0, 4X, F10.0,-4X, F10.0, 
1 4X, F10.0, 4X, F10.0, 4X, F10.0) " f is 

113 FORMAT < 1H1,1LX, -12HT0TAL INCOME, 5X, 18HT0TAL EXPENDITURES, 5X, 
1 7HSURPLUS, 5X, 15HSURPLUS/STUDENT, 5X, 15HEXPENSE/STUDENT ) 

114 FORMAT ( 5H0YEAR , 13, 5(7X, F10.0)) 

= — PRINT 106 . -^=—3; - • _. . _ _.. 

-" -" DO 52 1 = 1, L v-r_jr_^f^=:_" 

52 PR INT 107, T, B (I ) , : P„U ) , - Q ( I ), R(I), EU) .^11^^^:^ . „ 
^^r^^PR I NT^^ 

PRINT 109 

DO 53 1=1, L 

, 53 PRINT- 110,-1, X(1,I ),X(2,1-),X(3,I),X{4,I),X{5,_IJ- 

- • .... print ii i 1: -imm^m^^ 

DO 54 I=l~,t-: J™ IZ . _ 

54 PRINT U2,^-¥5G^^2^^ill^ 
PRINT 113 
DO 60 1*1, L 

Z ( I ) =W ( L) / AJJJ _____ -= _ ^ 

60 PR INT 1 l*- f >Ti iVdlr-UU) W ( \ )WH I ) , "IZ ("1 ) "~ 
CALL EXIT - _ -- - - - 7=- ------ - --" - 

END 




RESIDENT STUDENTS WHO 







STUDENTS 


STUDENTS 


BOARD ONLY 


RESEARCH 


TUITION+FEES 


ROOM 


BOARD 


GIFT INCOME 


ENDOWMENT 


YEAR 


1 


970. 


610. 


610. 


15000. 


1180. 


320. 


490. 


85500. 


4860000. 


YEAR 


2 


980. 


620. 


630. 


15000. 


1580. 


360. 


520. 


85500. 


5103000. 


YEAR 


3 


980. 


620. 


630. 


15000. 


1580. 


360. 


520. 


85500. 


5358000. 


YEAR 


4 


990. 


630. 


650. 


15000. 


1580. 


360. 


520. 


85500. 


5626000. 


YEAR 


5 


1000. 


640. 


670. 


15000. 


1580. 


380. 


550. 


85500. 


5907000. 


YEAR 


6 


1000. 


640. 


670. 


15000. 


1780. 


380. 


550. 


90000. 


6203000. 


Y FAR 


7 


x vj j \J . 






15000 . 




^ AO . 




~ UU \J\J . 


O D X J5UUU . 


YEAR 


8 


1200. 


750. 


790. 


15000. 


1780. 


400. 


580. 


90000. 


6838000. 


YEAR 


9 


1200. 


750. 


790. 


15000. 


1780. 


400. 


580. 


90000. 


7180000. 


YEAR 


10 


1210. 


760. 


800. 


15000. 


1780. 


400. 


580. 


90000. 


7539000. 


YEAR 


11 


1220. 


770. 


810. 


15000. 


1780. 


400. 


580. 


90000. 


7915000. 



RESEARCH OVERHEAD RATIO 1.32 
ENDOWMENT YIELD .0600 







RESEARCH INCOME 


TUITION+FEES INCOME 


ROOM INCOME 


BOARD INCOME 


ENDOWMENT INCOME 


TOTAL INCOME 


YEAR 


1 


19800. 


1144600. 


195200. 


298900. 


291600. 


2035600. 


YEAR 


2 


19800. 


1548400. 


223200. 


327600. 


306180. 


2510680. 


YEAR 


3 


19800. 


1548400. 


223200. 


327600. 


321480. 


2525980. 


YEAR 


4 


19800. 


1564200. 


226800. 


338000. 


337560. 


2571860. 


YEAR 


5 


19800. ; 


1580000. 


243200. 


' 368500. 


354420. 


2651420. 


YEAR 


6 


19800. 


1780000. 


243200. 


368500. 


372180. 


2873680. 


YEAR 


7 


19800. 


1869000. 


250800. 


379500. 


390780. 


2999880. 


YEAR 


8 


19800. 


2136000. 


300000. 


458200. 


410280. 


3414280. 


YEAR 


9 


19800. 


2136000. 


300000. 


458200. 


430800. 


3434800 . 


YEAR 


10 


19800. 


2153800. 


304000. 


464000. 


452340. 


3483940. 


YEAR 


© 


19800. 


2171600. 


308 0^^ 


469800. 


474900. 


35341^1 



o a 







FACULTY 


SALARY INCREASE 




NO. OF FACULTY 


STUDENT FACULTY RATIO 


YEAR 


1 




1.050 




67.0 


14.47 


YEAR 


2 




1.050 




68.0 


14.41 


YEAR 


3 




1.050 




70.0 


14.00 


YEAR 


4 




1.050 




75.0 


13.20 


YEAR 


5 




1.050 




76.0 


13.15 


YEAR 


* 




1.100 




76.0 


13.15 


YfAft 


T 




i.050 




78.0 


13.46 


YE Aft 


a 




1.050 




84.0 


14.28 


YEAR 






1.050 




86.0 


13.95 


YEAR 


10 




1.050 




90.0 


13.44 


YEAR 


11 




1.050 




98.0 


12.44 


COST OF LIVING FACTOR 
ADMINISTRATIVE SALARY INC* 
SCHOLARSHIP RATIO 


1.015 
1.044 
.250 












FACULTY 

FRINGE BENEFITS 
RATIO 




ADMINISTRATION 

CD IM^C HP MPC ITC 

RATIO 




YEAR 


1 




.070 




0.000 




YEAR 


2 




.070 




0.000 




YEAR 


3 




.070 




0.000 




YEAR 


4 




.070 




0.000 




YEAR 


5 




.070 




0.000 




YEAR 


6 




.070 




0.000 




YEAR 


7 




.070 




0.000 




YEAR 


8 




. .070 




0.000 




YEAR 


9 




.070 




0.000 




YEAR 


10 




.070 




0.000 




YEAR 


11 




.070 




0.000 





YEAR 
YEAR 
YEAR 
YEAR 

YEAR 
YEAR 
YEAR 
YEAR 8 
YEAR 9 
YEAR 10 
YEAR 11 



YEAR 
YEAR 
YEAR 
YEAR 
YEAR 
YEAR 
YEAR 
YEAR 8 
YEAR 9 
YEAR 10 
YEAR 11 
CONTINGENCY 



1 
2 
3 
4 
5 
6 
7 



AVERAGE FACULTY SALARIES 
10073* 
10576. 
11105. 
11660. 
12243. 
13468. 
14141. 
14848. 
15591. 
16370. 
17189. 

RESEARCH EXP. TOTAL FACULTY SALARIES 

674891. 



15000. 
15000. 
15000. 
15000. 
15000. 
15000. 
150pO. 
15000. 
15000. 
15000. 
15000. 



719212. 

777383. 

874556. 

930528. 
1023581. 
1103043. 
1247287. 
1340833. 
1473358. 
1684539. 



© 



FACULTY RETIREMENT 

47242. 

50344. 

54416. 

61218. 

65136. 

71650. 

77213. 

87310. 

93858. 
103135. 
117917. 



© 



SCHOLARSHIP FUND 
286150. 
387100. 
387100. 
391050. 
395000. 
445000. 
467250. 
534000. 
534000. 
538450. 
542900. 



SPECIAL COSTS 
0. 

100000. 
0. 

100000. 
0. 
0. 

100000. 
100000. 
100000. 
100000. 
100000. 



^^^^ 

















© 






BUDGET 

NU.l 


CATEGUR I ES NOT 
NO. 2 


SUBJECT TO 
NO. 3 


ENROLLMENT GROWTH 

NO. 4 NO. 5 






YEAR 


1 


589710. 


1066470. 


100000. 


0. 


0. 






YEAR 


2 


615367. 


1182467. 


101500. 


0. 


0. 






YEAR 


3 


642149. 


1200204. 


103022. 


110000. 


0. 






YEAR 


4 


670104. 


1218207. 


104567. 


111940. 


0. 






YEAR 


5 


699286. 


1236480. 


106136. 


113921. 


0. 






YEAR 


6 


729746. 


1255027. 


107728. 


115946. 


0. 






YEAR 


7 


761543. 


1273852. 


109344. 


118015. 


0. 






YEAR 


8 


7947 34. 


1292960. 


110984. 


120130. 


0. 






YEAR 


9 


829380. 


1312354. 


112649. 


122292. 


0. 






YEAR 


10 


865546. 


1332040. 


114338. 


124502. 


0. 






YEAR 


11 


903299. 


1352020. 


116054. 


12676 1. 


0. 










NO. 6 


BUDGET CATEGORIES SUBJECT 
NO. 7 NO. 8 


TO ENROLLMENT 
NO. 9 


GROWTH 

NO. 10 


NO. 11 


NO. 12 


YEAR 


1 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


YEAR 


2 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


YEAR 


3 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


YEAR 


.V 


. 


0. 


0. 




0. 


0. 


0. 


YE Aft 


5 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


YEAR 


6 




0. 


0. 


0. 


0. 


0. 


0. 


YEAR 


7 


150000. 


0. 


0. 


. 


0. 


0. 


0. 


YEAR 


8 


175657. 


0. 


0. 


0. 


0. 


0. 


0. 


YEAR 


9 


180022. 


0. 


. 


0. 


0. 


0. 


0. 


YEAR 


10 


186066. 


0. 


0. 


0. 


0. 


0. 


0. 


YEAR 


11 


192335. 


0. 


0. 


0. 


0. 


0. 


0. 







TOTAL INCOME 


TOTAL EXPENDITURES 


YEAR 


1 


2035600. 


2779463. 


YEAR 


2 


2510680. 


3170991. 


YEAR 


3 


2525980. 


3289276. 


YEAR 


4 


2571860. 


3546645. 


YEAR 


5 


2651420. 


3561489. 


YEAR 


6 


2873680. 


3763680. 


YEAR 


7 


2999880. 


4175262. 


YEAR 


8 


3414280. 


44 78064. 


YEAR 


9 


3434800. 


4640391. 


YEAR 


10 


3483940. 


4852437. 


YEAR 


11 


3534100. 


5150827. 



SURPLUS 


SURPLUS /STUDENT 


EXPENSE/ STUDENT 


-743863. 


-766. 


2865. 


-660311. 


-673. 


3235. 


-763296. 


-778. 


3356. 


-974785. 


-984. 


3582. 


-910069. 


-910. 


3561. 


-890000. 


-890. 


3763. 


•1175382. 


-1119. 


3976. 


•1063784. 


-886. 


3731. 


•1205591 . 


-1004. 


3866. 


1368497. 


-1130. 


4010. 


1616727. 


-1325. 


4221. 
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Program Abstract 

Austin College Symbolic Programming System 

David Ro Musser 
Austin College 
Sherman, Texas 75091 



Users Group Code : 5252 

Direct Inquiries to: Mr Robert F Randall, Director 

Hoblitzelle Computer Center 
Austin College 
Sherman, Texas 75091 
Phone: 214-^92-910 1 Ext . 40 

Subject Classification : 1.1 

Description/Purpose : 

AC-SPS is an SPS system designed to provide users of a 1620 
without a disk-file some of the versatility, speed and convenience 
of disk operation. This is achieved by modifying an existing SPS 
system, AFIT SPS (Id. 027 , written by R. L. Pratt), which is itself 
an improvement of the standard IBM system, especially in terms of 
speed of operation. Essentially all of the features and capabilities 
of AFIT SPS have been retained, but additional features are provided: 
(1) a supervisor urogram which controls processing by interpreting 
control cards; (2) load-and-go operation for small programs, which 
reduces the time required for debugging programs; (3) macro instruc- 
tion subroutines for input-output operations (READ, ACPT, PNCH, TYPE) , 
which provide an 1/0 system similar to that of Basic Fortran, greatly 
simplifying 1/0 programming . 



© 



Specifications : 
Storage : 
Equipment 



40K or larger (self-adjusting) 

Card system, TNS , TNF, MF, Auto-divide, Indirect 
addressing. Requirements for TNS, TNF, MF and 
Auto -divide could be removed fairly easily by 
reprogramming . 



Language : AFIT SPS. Writeup is in English. 



o 
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PROGRAM DESCRIPTION 



Introduction 



Austin College SPS is a modification of AFIT SPS (1,0.027— 
"AFIT Symbolic Programming System, Modified for 80- $0- Listing"). 
Originally AC-SPS was developed as a modification to IBM SPS, 
but the resulting system left something to be desired with regard to 
speed of operation . It was then found that the AFIT SPS system, 
written by Mr. R. L. Pratt, was much superior to IBM SPS, 
especially in terms of speed. In fact, AFIT SPS assembles three 
to four times as fast as IBM SPS, for the following reasons: 

A binary search of the opcode and label tables is used. 
Indirect addressing and W , TNF, TNS are used. 



$0-80 list board is produced (one listing card per 
source program card, usually); thus the deck requires 
only half as long to punch, half as long to list. 
During the second pass, the source program is retrieved 
from memory rather than reread (if it was possible 
to store the entire source program during pass one.) 

Aside from its greater speed, AFIT SPS also provides addi- 
tional language features, which in a number of instances make 
programming more convenient . 

AC-SPS retains the speed of operation and essentially all 
the capabilities of AFIT SPS, but makes it more versatile and 
more convenient to operate by the addition of the following 
features: 

(1) A supervisor program, which controls the processor by 
interpreting control cards included in the source pro- 
gram deck by the programmer. 

(2) Load -and- go operation for short programs. Only a few 
minor restrictions are placed on the source language 
permitted when the processor is operated in the load- 
and-go mode. All essential floating point macro in- 
structions are allowed, assuming that the special load- 
and-go subroutine deck has been loaded. 

(3) Input -output macro instruction subroutines, which pro- 
vide an I/O system similar to the system of Basic Fortran, 
thus greatly simplifying I/O programming. These macros 
are available for either load-and-go operation or inclu- 
sion in punched object decks. 




with a standard 



/// 



The first two of the?se features combine to. give the system 
some of the power and convenience of a disk SPS processor under 
the control of a monitor system . The third feature makes SPS a 
more useful language, as it becomes somewhat more competitive 
with Fortran with regard to ease of programming (especially for 
short programs which formerly often required more programming 
effort for I/O than for actual computation), while retaining its 
greater versatility. 

The present description of AC-SPS will deal only with the 
more significant aspects of the features which have been added 
to AFIT SPS . Familiarity with the AFIT SPS system would be 
helpful, but is not essential to the understanding and operation 
of AC-SPS, if the reader is familar with the Reference Manual , 
IBM 1620/1710 Symbolic Programming System , C 26-5600. When the 
program is submitted to the Users Group Library, a more thorough 
description of all features will be included in the writeup. 




Supervisor program 



The use of a supervisor program was motivated mainly by the 
desire for more convenient operation (especially since more options 
were provided, making too many to reasonably handle with switch 
settings). It is also felt that this system can be used to in- 
troduce programming students to the supervisor program concept, 
which, of course, they will need to become familiar with if -they 
later work with a larger machine or a 1620 with a disk-file. 

Most, but not all, of the control of the operation of the 
processor is determined by control cards prepared by the pro- 
grammer and included in his source deck. In general, all options 
over which it is felt the programmer should have advance control 
have been defined as control card options, and a few remaining 
options have b^en left under manual control of the console switches. 
The latter category consists of options which are not often used, 
and/or which require action by the console operator; e.g. input 
or error correction from the typewriter, and listing on the type- 
writer. These options are controlled by Switches 1, 2, and 3 
(see the operating instructions). If any of these options is 
desired, then the processor should be operated with switch 4 on , 
which causes the processor to halt after each stage of processing 
so that Switches 1-3 can be reset if necessary. Generally, how- 
ever, none of these options is desired, and the processor may be 
operated with all four switches off . In this case, there is usu- 
ally no need for operator intervention, and the machine does not 
halt unless there is. This generally means faster, more efficient 
operation when a number of short programs are being processed. 

The supervisor-processor system processes .jobs . A job is 
defined here to be deck of cards composed of either: 



Job card* 

Source Source program deck*, including source program 
Program control cards 

Job Data cards (if the program is to be executed and 

contains read statements) 
End-of-file card 

or 



Job card* 

Object Object program control cards 
Program Object program* 
Job Data cards 

End-of-file card 



*Must be present in all jobs of this type. 
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The control cards are described in the operating instruc- 
tions . It may be noted here that the simplest possible source 
program job consists of 

Job card 

Source program deck, without any source program 
control cards 

For such a job, the program will be assembled into memory, 
and if no errors are detected, it will be executed. This mode 
of operation, called the load-and-go mode, is described in the 
next section. An object deck, a listing deck, etc., may be 
obtained, and the operation of the system may be otherwise modi 
fied, by the inclusion of one or more of the source program 
control cards a 

Similarly, the simplest possible object program job is 

Job card 
Object program 

In this case the object program is loaded under control of 
the supervisor program loader, which protects the supervisor- 
processor system, and the subroutines (if they are in memory). 
Control cards are provided to remove some or all of this protec 
tion when necessary, or to allow the loading of a non-AC-SPS 
assembled object deck. 
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Load-and-go mode 



In this mode of operation no object deck is punched; rather, 
the object program is loaded into a work area of memory as it is 
assembled. This saves time and card handling, and permits complete 
processing of a job, including execution of the object program 
with data after it has been assembled. This is especially useful 
when debugging a program. 

A limited attempt has been made to protect the supervisor- 
processor from being damaged by errors in the program during 
execution, and to make operation as continuous as possible. 

(1) The area of memory into which the object program is 
loaded is initialized to flags and record marks, so 
that undefined variable errors will not clear memory. 

(2) No attempt is made to check for any of the countless 
possible errors which can cause check-stops, but an 
easy restart procedure is available when they occur. 

(3) The halt statement is modified so that control is re- 
turned to the supervisor when it is executed, rather 
than an actual halt occuring. Also, the input macros 
provided with the system can detect an end-of -file 
card and return control to the supervisor. Thus pro- 
grams terminated with either a halt statement or 
running out of data cards will not halt the automatic 
operation of the system. 



These and certain other aspects of the load-and^go mode are 
discussed more fully in the operating instructions. 

The remainder of this discussion is concerned with restric- 
tions on the size of a program run in this mode, and minor res- 
trictions on the language permitted in the source program. 

Size of the program 

On a 40K machine, a program which requires the subroutines 
is restricted to 7000 digits length; one which does not may be 
16000 digits long. 



Language 

TRA and.TCD statements and the DNB statement are not allowed 
in the load-and-go mode. ji • • * 

Generally, the DORG statement should not be used to define/ 
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the program itself, as the processor automatically assigns it 
properly, and changing this assign/ment may result in loading of 
the object program being inhibited. 

It is even more important that absolute addresses in the 
source program be avoided, since if such appear in instructions 
then the supervisor-processor may be clobbered when these 
instructions are executed. (The processor does not check for 
this . ) 

The DEND statement must have as its operand the label of the 
instruction where execution is to begin; if this operand is 
omitted then execution of the program will be bypassed. 

The load-and-go subroutine deck contains most, but not all 
of the standard floating point subroutines. DIV, AT AN, FSLS 
and FSRS have not been included. Subroutines which were added 
to the AFIT system, INC, OUTC, FC and CKPP, have also been 
deleted. 



Input - output macros 



The INC and OUTC macros were added to the AFIT SPS system 
to make I/O programming more convenient than it is using the 
hardware I/O commands. In AC-SPS these are replaced by an even 
more powerful set of macros, READ, ACPT, PNCH and TYPE, which 
provide an I/O system which is similar to the system used in 
Basic Fortran. No attempt will be made here to describe in 
detail the usage of these macros. The discussion will be in- 
stead limited to a simple example. 

Consider the following Fortran statements: 



Corresponding to these Fortran statements are the following 
AC-SPS statements: A 



Fl DSA 61004,61002,51407,90600,91000, 

Assuming that the quantities A,B, etc., were correctly 
stored, the output produced by the SPS statements would be 
exactly the same as that produced by the Fortran statements. 
This example illustrates several aspects of the AC-SPS I/O macros: 

1. The operands field contains first the label of a "Format 
declaration," then any number of labels (or absolute 
addresses or any other valid form of operand for a macro 
instruction), separated by commas. Note that the last 
operand also has a trailing comma , which cannot be omitted. 
The I/O subroutine itself uses the five digit zero field 
produced by this comma to determine the end of the list 

of variables. 

2. The "Format declaration" is a labeled DSA statement 
with five digit codes (these are not addresses) as its 
operands, representing format specifications. Permitted 
numeric conversion specifications include: 



PUNCH 1, A, B, C, I, J, P, Q 



1 FORMAT (F10.4,F10.2,E14.7,IS, 



110) 



PNCH 




Code 



Corresponding Fortran 
Specification 

E 



9 /xy 



There is also a TT D TT specification, code 4? which stores 
numbers in fixed point fields on input but inserts a 
decimal point in numbers output from a fixed point field. 
Code 7 can be used to obtain the same effect as the non- 
numeric X specif icatbn in Fortran. 

This code digit is the first digit of the five digit code. 
The next two digits specify the width of the field in 
which the number is contained in external form. The 
last two digits specify the position of the decimal 
point for the "D", Tr E TT , and n F TT specifications. Thus the 
form of the entire specification is 

cwwdd 

3. Note the presence of the trailing comma in the format 
declaration also. This has somewhat the same effect as 
the right parenthesis in the Fortran Format statement 
when the number of specifications is less than the number 
of variables in the list; i.e. a new record (card or type- 
writer line) is skipped to, and the beginning of the format 
declaration is returned to for specifications to use with 
the remaining variables. 

4. Format DSA T s, like most other constants in an SPS program, 
cannot be interspersed with instructions in the way that 
Format statements can in a Fortran program. 

The input routine reads free form data; that is, field widths 
are not observed and numbers to be input may be positioned in any 
manner, except that they must be separated by one or more blanks. 
(In fact, the input routine is closely patterned after the input 
routine in the PDQ Fortran Free Form Subroutines . ) The E and F 
specifications may be used interchangeably to convert the standard 
E and F forms to floating point fields, and the I and D specifica- 
tions interchangeably to convert numbers with or without a decimal 
point (if present, it is ignored) to fixed point fields. The I 
and D specifications use the dd code digits to determine the length 
of the field in which to store the numbers. 

Use of the output routine is very much the same as use of 
Basic Fortran output routines, except that when the width of an F 
specification is exceeded, the number is automatically output in 
E form. Care must be taken, however, when using the I and D forms 
that ww is sufficiently large for the length of the field being 

output. 
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Student Data Processing System at Christian 
Brothers College Employing an IBM 1620 

Brother Jerome David Wegener, Registrar 
Christian Brothers College, Memphis, Tennessee 



Christian Brothers College is a small all men's college. Last year we began 
with 950 students, and this year we expect to have about 1,050. We have an unusual 
situation in that nearly 60% of our students are majoring in engineering or science. 
It was for the instruction of these students that a 1620 was installed in January of 
1962. The computer center is run on a completely open shop basis both for students 
and faculty. As. a member of the faculty and Registrar, I learned Fortran and SPS, and 
we began using the computer center for records in January of 1963. 

The computer center consists of an IBM 1620, Model II, a 407, an 082 sorter, 
and an 026 key punch with a special Interspersed Gang Punch feature. Previous to 
this time no machinery at all was used in the Registrar's office. With the exception 
of a typewriter, all the work was done by hand. We did not graduate from a Unit 
Record operation, though our present system probably resembles that greatly since 
it is a completely card-oriented system. The computer is used as a substitute for 
many of the ordinary hand and machine operations. 

During the second semester of 1962-63, we paralleled the hand system with the 
computerized system, and in the fall of 1963 we went over to the computerized system 
completely. Besides taking care of the grades and other records, the system is also . 
used to make out bills for each student. 

I. INITIATING A STUDENT: 

The first step in getting a new student into the system is to prepare an admission 
card. This is done by information taken from the student's application blank and 
transcript. When several such cards are available, a special program is used to 
assign student numbers according to alphabetical order. It places the new student 
halfway between the numbers of two students who are already in the file. The output 
of this program is a special card which carries the same information as the admission 
card, with the addition of the student number. These are then interpreted and kept 
in a special file containing all students. 

By duplicating certain parts of this card, a #1 card is prepared. This is the 
student's name card, and later carries much additional coded information. For students 
who have been in school the previous semester, these cards can simply be updated. 
A special program has been written which updates the year and duplicates all the cards 
except those of graduating seniors (whose cards are not duplicated). Certain other 
changes must be made by hand, such as changing the probationary or boarding status 
of a student, or his major. 

With the #1 card as a starting point, the other basic cards are produced by key 
punching the appropriate information. These include the #2 card, used for billing 
address; the #8 card, used for parent's or guardian's address; and the #9 card, used 
for local address and telephone number. A special program produces from the #1 
card, the #7 card, which contains the name of the student's major and the high school 
he attended. Besides being used for addressing labels, all these cards are used in 
printing the student information form. 
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At the present time, the Master Schedule is made out completely by hand. No 
attempt has been made to computerize this or to do sectioning, mainly because no 
one has had the time to devote to it. 

Once the Master Schedule is made out, the list of teachers is alphabetized, 
and assigned numbers in order beginning with number one. The program at present 
allows for 99 instructors. Next the courses are numbered, alphabetically by depart- 
ment, beginning with number one. At present, the program allows for 240 course 
sections. The above mentioned numbers change from semester to semester and are 
used mainly by the computer in keeping track of totals. Department numbers and 
major field numbers remain the same from semester to semester. 

With this information, the Master Course cards are made out, one for each 
section listed on the schedule. These are identical to the Student Grade Card (#3 
card) with the exception of putting the teacher's name in place of the student's name. 
These are punched on right corner cut cards in order to make use of the interspersed 
gang punch. 

Blank Grade Cards are counted out by the computer for each section according 
to the estimated number of cards for each section. Each group is placed after the 
proper Master Course Card and these are then gang punched, leaving the place for 
the student's name and number blank. 

III. REGISTRATION: 

At the time of registration, the students fill out a single registration card 
listing their schedule and certain information. With this in hand they proceed to 
the gym or other designated area, pick up their name card and move around the gym 
picking up a Student Grade Card for each course that they have listed on their schedule 
At the end of the line, the cards are checked to see if that student has made no mis- 
takes. This entire process takes about ten minutes. 

The name card is a right-corner cut card, and is kept in front of the student's 
Grade Cards. These are then gang punched making use of the Interspersed Gang 
Punch features on the key punch. Once the student's name and number have been 
punched on all the cards, the name cards are separated from the Grade Cards, and 
these latter are sorted into alphabetical order by course along with the Master Course 
Cards, and the temporary class lists are printed. Once copy is given to each instruc- 
tor for each of his classes, and a complete set to the Dean's office. 

The cards are then re-sorted into alphabetical order by student. In this sort, 
all the cards numbered 1,2,3,7,8, and 9 are used, and the student's schedule and 
information sheet is printed. One copy is given to the Dean's office and one to the 
Business Office. A copy of the information sheet is given to each student so that he 
can report any errors in addresses, or schedule to the Dean's office. The students 
have ten days in which to finalize their schedules. During this time, any changes 
are recorded on Change Slips and the appropriate change is made among the student's 
Grade Cards. At the end of this time, the Final Class Lists are printed, as well as 
new student information sheets. 
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After the registration has been finalized , distributions are calculated according 
to the various categories, such as by Major, by state or country, by year in school, 
by boarding status, and by high school attended. A special program is used to do 
this. Listings are also made according to these various categories. 

Another program is used to count the number of Grade Cards that each student 
has, and at the same time determine those who are special students. A new name 
card is punched which contains this additional information. This is done so that, 
with all the handling of the cards, we can check to see that none is lost. 

IV. AT THE MARKING PERIODS: 

A week or so before marks are due, the Grade Cards are sorted into order by 
courses and distributed to the teachers along with the Grade Sheets. The teachers 
record the marks on the cards and on the sheets, and return them to the Registrar's 
office. The cards are hand sorted into groups by grade, and the computer is used to 
gang-punch the grades into the cards. The F's are done first so that the failure list 
can be printed. 

The Grade Cards are then sorted into order by student, along with his name and 
address cards. These are run through the computer which, under the control of a 
special program, calculates the quality points and quality point ratio. This program 
has an elaborate system of checks built into it to see that all the cards are present, 
that each has a proper grade on it, etc. The computer re- punches the Grade Cards 
with the proper quality points on it. Credit hours are corrected in the case of a 
failure or withdrawal. The computer then punches the #4 card which gives the total 
hours and quality points for the semester, and the #5 card which gives the cumulative 
hours and quality points, adjusting for repeated courses when necessary. Finally, 
the #6 card is punched giving the student's current qtiality point ratio along with any 
appropriate action such as "Dean's List," "Probation Removal/ 1 etc. These are then 
printed on the student's report card. 

Four copies of the report card are made: one for the student, one for his parents, 
one for his high school, and one is kept in the Registrar's office. Special copies of 
the reports for students in each major field are printed and given to the head of each 
department. Copies of Boarders' marks are given to the dorm prefects, and labels 
are printed for attaching to the student's permanent record. 

The computer has kept track of the total number of grades for each department, 
course and teacher. These are printed and dittoed, and distributed to appropriate 
offices and teachers. From the #6 cards, Probation Removal list, Probation List, 
and Probation Failure List are prepared, as well as the Dean's List. 

V. BILLING PROCEDURE 

The same basic set of cards is used in preparing each student's bill. A special 
program has been written, which, from the coded information on the student's name 
card, causes the computer to punch out special cards for the various charges such 
as tuition, room and board, matriculation and graduation fees , etc. These are then 
sorted with the students Grade Cards which contain any appropriate lab fees. Special 
students are charged by the semester hours, and the computer program calculates the 
tuition for those special cases. The computer keeps totals for each of the various 
charges, and for the lab fees by department, so that the business office has a check 
on the amount and distribution of the incoming funds. 
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Any credit that the student has with the business office by means of previous 
payments is entered on the bill by hand. A copy of the bill is sent to the student's 
parents, and a copy kept in the business office. 

VI. CONCLUSION 

We are very pleased with the system just described. The 1620 is not exactly 
a data processing machine nor do we have all the equipment that we might have. We 
do not feel that our small enrollment justifies a larger data processing machine , but 
we do feel that we are making very good use of the equipment that we do have on 
hando Although it still requires a large amount of hand work, it is not nearly as 
much as we did previously, and the amount of time spent on the various records has 
been cut quite drastically. 




MATRIX STRUCTURAL ANALYSIS 



Everett L. Cook 

Wichita State University 
Wichita , Kansas 

The program described herein may be used to analyze struc- 
tures by either the force method or the displacement method. 
The method of analysis is specified by a control card at execution 
time. The input data is, of course, a function of the method to 
be used. 

The program has two unique features: 

1. A recursive technique is used to solve for 

the redundants or kinematic deficiencies 
rather than a direct matrix inversion. 

2. The flexibility or stiffness matrix of the 

unassembled structure is stored in a 
special compressed format. 

The program presented here was compiled in PDQ Fortran and 
run on a 40k 1620 with card input-output. However, the original 
version was a four-part program for a 2 0k 16 2 (Ref . 2) . 

Matrix Force/Displacement Structural Analysis 

The analogy between the matrix force and matrix displacement 
methods of structural analyses was first tabulated by Argyris 
and Kelsey (1). Pestel and Leckie (3) have pursued this analogy 
and have developed an excellent notation which will be used here. 

The analogy is based on the inverse relationship of ~ the 
following basic parameters: 

External Forces = {f} {d} = External Displacements 

Internal Fore js = {p} (v} = Internal Displacements 

Element Flexibility = [F] [K] = Element Stiffness 



Sponsored by the U. S. Army Research Office - Durham under 
Grant No. DA-ARO (D) — 31-124-C482 . 
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Since the force method is easier to explain, the basic 
equations for this method will be stated and explained, and then 
the analogy with the displacement method will be shown. The 
two equations of final interest in either analysis are 

{d} = [P d ]{f> (1) 

and 

{p} = [B]{f}, .(2) 

since Equation (1) expresses the external displacements in terms 
of the known external forces and Equation (2) gives the internal 
forces in terms of these same external forces. 

The matrix equations which must be solved in order to obtain 
the two rectangular matrices, [F^] and [B] , are as follows: 

[D 1Q ] = [B 1 ] T [P v ] [B Q ] (3) 
[D 11 ] = [B 1 ]T [F v ] [B 1 ] (4) 

[x] = -[ D nl" 1 [ D io 1 (5) 

[B] = [B ] + [B x ] [X] (6) 

[F d ] = [B] T [F v ] [B] (7) 

Three of the matrices in the above equations must be formu- 
lated from the physical and geometric properties of the structure: 

[B Q ] = a set of internal forces in the structure due to 
unit values of the external forces. Each column 
represents a set of internal forces due to one 
of the external loads. 

[B 1 ] = a set of internal forces in the structure due to 

unit values of the redundants. Each column represents 
a set of internal forces due to one of the redundants. 
For a statically determinate structure [B 1 ] = [0] 



and [B] = [b q ] . 
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[F ] = the flexibility matrix of the unassembled structure, 
This is a diagonally partitioned matrix where the 
submatrices are the element flexibility matrices. 

These matrices are the input required for a matrix force 
analysis utilizing the program being discussed. 



Example: 



EI = 50 x 10 6 lb. in. 2 



30 



F = 0.0002 in. /lb. 



15 



»P< 



<3 



Assume that the force in the spring, p 5 , is the redundant. 



[B Q ] 



^1 
30.0 

0.0 

0.0 

0.0 

0.0 



2 

-1.0 
1.0 
0.0 
0.0 
0.0 



X. 



= 10 



-6 



0.20 
-0.10 



P3 
P4 
P5 
-0.10 
0.20 



0.10 
-0.05 



-45. 





Pi 


15. 





?2 


-15. 





P3 


0. 





P4 


1. 





P5 



-0.05 
0.10 



200, 



Pi 
P2 
P3 
P4 
J p 5 



The matrix displacement equations equivalent to Equations 
(3) - (6) are 

[C 10 ] = [Aj^Kp] [A Q ] (8) 

[C X1 ] = [A 1 ] T [K p ] [A x ] (9) 

[Y] - -tC 11 ]" 1 tC 10 ] (10) 

[A] = [A Q ] + [A x ] [Y] (11) 

[K f ] = [A] T fF v ] [A] (12) 

Two additional equations are required in order to obtain the 
external displacements and the internal forces with Equations (1) 
and (2) . 

[P d l = [Kj]" 1 (13) 

[B] = [K ] [A] [K J" 1 (14) 

The three input matrices required for a displacement 
analysis are 

[Aq] = the internal displacements due to unit displacements 
at the external loads. 

[A^] = the internal displacements due to unit displacements 
at the kinematic deficiencies. A kinematic defi- 
ciency exists for each degree of freedom for -which 
there is no external load. If there are loads 
corresponding to all unconstrained external displace- 
ments, [A.J = [0] and [A] = [A Q ] . 

[K ] = the stiffness matrix for the unassembled structure. 
P 
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Example: 




[A Q ] = 



[K ] 
P 



d l 


d 2 






Y3 


*4 




0.033333334 


0.0 


v l 




0.0 


0.0 


v i 


0.033333334 


1.0 


V 2 




0.0 


d.o 


V 2 


-0.066666667 


1.0 


V 3 


[A x ] « 


0.066666667 


0.0 


V 3 


-0.066666667 


0.0 


V 4 




0.066666667 


1.0 


V 4 


0.0 


0.0 


V 5 




1.0 


0.0 


V 5 



= 10' 



6.6666666 
3.3333333 



3.3333333 
6.6666666 



13.333333 6.6666666 
6.6666666 13.333333 



0.005 



v. 



v . 
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Recursion Analysis 

The direct approach to the solution of the matrix force or 
displacement equations requires the inversion of either tD^l*' 
or [C-q] • Since the order of [D^] is equal to the number of 
redundants and the order of [C-q] is equal to the number of 
kinematic deficiencies , the required inversion may be relatively 
large. Pestel and Leckie (3) have suggested a recursive tech- 
nique whereby the redundants f or deficiencies , are divided into 
a number of subgroups, and an equal number of lower-ordered 
inversions are performed. In attempting to program this recursive 
technique, it was discovered that, not only is the programming 
simplified, the data storage requirements are minimized if the 
redundants, or deficiencies, are considered one at a time. 
Therefore, if there are n redundants, or deficiencies, there 
are n one-by-one inversions or, in fact, no inversions at all. 
It should be noted that a displacement analysis does require an 
inversion; however, the maximum order is five for the present 
program. 

The basic recursive equations for the force method are as 
follows: 

{D oi } - {B i }T[F v ][B oi ] (15) 

{X 1 } = -{Dq }/D^ (16) 

fB^ 1 ] = [Bj] + {BjHX 1 } 1 (17) 

Equations (3) and (4) have been combined to obtain Equation 
(15). These equations could be written 

[D 10 l D lll = [B 1 ] T [F V 3[B |B 1 1 (18) 

For the first recursion, only the last column of the 
matrix [B-jB*] is taken as the effective {B*}. As a result, 



both { D 01 } and ^ x } are column vectors rather than rectangular 
matrices. The divisor in Equation (16) is the last element 
in {Dq^} and the redundants are actually computed and stored in 
{D^}. The [Bq^ 1 ] calculated with Equation (17) has one 
less column than [B^] and f after n recursions, [Bq^] is the 
unit internal force matrix [B] . 

The recursion equations for the matrix displacement method 
can be written directly from the analogy as 

(19) 
(20) 
(21) 



{C 01 } = {A^^KpltAjj 

{y 1 } = -{cj}/cj; 

{A^ 1 } = [Aj] + {aJHY 1 } 1 



Storage of [F ] or [K ] 
- — X E — 

One of the major problems encountered in matrix structural 
analyses of the type being described here is the large size of 
the stiffness or flexibility matrix for the unassembled structure. 
However, since these matrices are diagonally partitioned, it is 
possible to store them in a compressed format. The present 
program uses an array with the number of rows equal to the maxi- 
mum number of internal forces (50) and five columns. The number 
of columns restricts the size of the element flexibility, or 
stiffness, matrix which can be accommodated. Five columns were 
chosen in order to be able to analyze structures which contain 
rectangular or trapezoidal plate elements. 

The element stiffness, or flexibility, matrices are right- 
justified in the array and the routine for calculating the 
product {A^} T [K ], or {B^} T [F v ], determines the element size by 
searching a row until it finds a nonzero coefficient. Since the 
premultiplication matrix is a row matrix, the multiplication is 
reasonably simple. 
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Results for Example 



FORCE METHOD 



DISPLACEMENT METHOD 



INTERNAL FORCES 



INTERNAL FORCES 



+ • 1 2445820E + 0-2 
+.58513932E+01 
-•58513932E+01 
+*00000000E~50 
+ #39009288E4-00 



+#30960000E-02 
+.66563470E+00 
+/33436533E+00 
+#000D0000E-50 
-•22291022E-01 



+• 12445820E+02 
+#58513930E+01 
-•58513927E+01 
+ # 57l76470E-06 
♦•39009287E+00 



♦•30959600E-C2 
+ •66563470E+00 
+ • 33436528E+00 
-•875294 1 OE-07 
-•22291022E-01 



DISPLACEMENTS 



DISPLACEMENTS 



+.57120742E-04 



1 9783282E-05 



+.57120743E-04 



19783282E-05 



1 9783281E-05 +• 198761 62E-06 1 978328 1 E-05 +• 1 9876 1 6 1E-06 
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MATRIX FORCE/DISPCACEMENT PROGRAM 
E.L. COOK AND G.E. LAMBERT 
REVISED MARCH 1966 



THIS PROGRAM MAY BE USED TO ANALYZE STRUCTUTES BY EITHER THE 
MATRIX FORCE METHOD OR THE MATRIX DISPLACEMENT METHOD • THE STRUCTURE 
MAY HAVE ANY COMBINATION OF DIFFERENT TYPES OF ELEMENTS AS LONG As NO 
ELEMENT HAS A FLEXIBILITY OR STIFFNESS MATRIX LARGER THAN FI VE-BY-Fl VE. 

OTHER RESTRICTIONS ARE- 



1 • MAX I MUM NUMBER OF ELEMENTS * SO 

2. MAXIMUM NUMBER OF INTERNAL FORCES « 50 

3. MAXIMUM NUMBER OF EXTERNAL FORCES « 5 

4. MAXIMUM NUMBER OF EXTERNAL FORCES 

PLUS REDUNDANTS OR KINEMATIC DEFICIENCIES « 28 



THE INPUTS TO THE PROGRAM ARE AS FOLLOWS— 



************ PARAMETER CARD ~ NUMBER OF ELEMENTS ( K ) « NUMBER OF 

* * INTERNAL FORCES <L>. NUMBER OF EXTERNAL FORCES CM>« 

* CARD 1 * NUMBER OF REDUNDANTS OR DEFICIENCIES (N)t NUMBER OF NON— 

* * ZERO ELEMENTS IN BO OR AO ( NZEBO ) « NUMBER OF NONZERO 
************ ELEMENTS I N Bi OR Al <NZEB1)« AND A DESIGNATION FOR THE 

METHOD BEING USED ( KK ) • FOR THE FORCE METHOD* KK» 1 • AND 
FOR THE DISPLACEMENT METHOD* KK*2. THESE PARAMETERS MUST BE IN FIXED 
POINT NOTATION (NO DECIMAL POINTS) AND MUST BE SEPARATED BY AT LEAST 
ONE BLANK • I.E.* 



12 24 2 11 T 54 1 



************ ELEMENT FLEXIBILITIES OR STIFFNESSES - THESE CARDS 

* * CONTAIN THE ELEMENT FLEXIBILITIES OR STIFFNESSES IN A 

* CARDS 2 * SPECIAL FORMAT* EACH CARD MUST CONTAIN FIVE NUMBERS IN 

* THRU L+l * FLOATING POINT NOTATION (DECIMAL POINTS REQUIRED) • A 

* * TWO-FORCE MEMBER WILL HAVE A SINGLE CARD WITH FOUR ZE*OS 
************ FOLLOWED BY THE FLEXIBILITY OR STIFFNESS OF THE ELEMENT— 



0.0 0.0 0.0 0.0 .66666667E-06 



A BEAM ELEMENT ( SAY PINNED AT THE ENDS WITH MOMENTS AT BOTH ENDS) WILL 
HAVE TWO CARDS WITH THREE ZEROS FOLLOWED BY THE ELEMENT FLEXIBILITY 
OR STIFFNESS* I.E.* 

0.0 0.0 0.0 .111111 11E-06 -.05555555E-06 
0.0 0.0 0.0 -.05555SB5E-06 .1111111 1E-06 

AND SO ON UP TO A PLATE ELEMENT WHICH WILL HAVE FIVE CARDS CONTAINING 
AN ELEMENT FLEXIBILITY OR STIFFNESS MATRIX OF THE FORM - 

.40OOOE-05 .29999E-05 .99999E-06 -.33333E-06 -.33333E-06 
•29999E-05 .B6666E-05 .56666E-05 .29999E-05 -.29999E-05 
•99999E-06 .56666E-05 .66666E-05 .26666E-05 -.33333E-05 
-.33333E-06 .29999E-05 .26666E-05 .39999E-05 -. 19999E-05 
-.33333E-06 -.29999E-05 -.33333E-05 -.19999E~05 .39999E-05 



************ NONZERO ELEMENTS IN BO OR AO - EACH OF THESE CARDS CONTAINS 

* NEXT * ONE NONZERO ELEMENT OF BO OR AO PRECEEDED BY NUMBERS 

* NzEBO * INDICATING THE ROW AND COLUMN IN WHICH IT APPEARS. 

* CARDS * THE ROW AND COLUMN ARE IN FIXED POINT NOTATION AND THE 
************ ELEMENT IS IN FLOATING POINT NOTATION, I *E* « 

10 3 -.70710680E+01 



************ NONZERO ELEMENTS IN Bl. OR A I - THESE CARDS ARE REQUIRED 

» NEXT * ONLY FOR STATICALLY INDETERMINATE OR K I NEMATI CALLY 

* NZEBl ♦ DEFICIENT STRUCTURES. THEY CONTAIN THE NONZERO ELEMENTS 

* CARDS * OF Bl OR Al IN THE SAME FORMAT AS THE CARDS IN THE 
************ PRECEDING GROUP. 



GENERAL OPERATING INSTRUCTIONS 

TO LOAO A NEW DECK INTO THE COMPUTER • PLACE THE DECK IN THE 
1622 READ HOPPER* DEPRESS INSTANT STOP AND RESET ON THE 1620* AND 
DEPRESS LOAD (NOT READER START) ON THE !622» 

THE 1622 ALWAYS STOPS PRIOR TO READING THE LAST TWO CARDS OF A 
DECK. DEPRESS READER START ON THE 1622 TO REAO THESE LAST TWO 
CARDS • 

TO PUNCH CARDS* BLANK CARDS MUST BE PLACED IN THE PUNCH HO**PER 
OF THE 1622 AND PUNCH START MUST BE DEPRESSED* WHEN THE 1622 STOPS 
PUNCHING* THERE ARE STILL THREE CARDS IN THE MACHINE (ONE PUNCHED AnO 
TWO BLANK). TO GET THESE CARDS * EMPTY THE PUNCH HOPPER AND DEPRESS 
NONPROCESS RUNOUT ON THE PUNCH END OF THE 1622* 

TO STOP EXECUTION OF A PROGRAM* DEPRESS INSTANT STOP ON THE 1620. 

TO RESTART A PROGRAM. DEPRESS RESET* INSERT* RELEASE* AND START 
ON THE 1620. 

to empty the read un it*. remove all cards from the read hopper* 
depress reader stop on the 1622* and hold the 1622 reader nonprocess 
Runout switch in the on position until two cards fall into the reject 

STACKER. 
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STEP-BY-STEP PROCEDURE 

ON THE 1620 CONSOLE* SET THE PARITY AND I/O SWITCHES TO STOP AND 
THE OFLOW SWITCH TO PROGRAM. SET PROGRAM SWITCHES 1* 2. 3* AND 4 
TO OFF FOR NORMAL OPERATION. IF SWITCH 4 IS SET TO ON* A PROGRAM 
TRACE WILL BE PUT INTO OPERATION. 



2. PLACE THE PROGRAM DECK IN THE READ HOPPER OF THE 1622. 

3. PLACE THE INPUT DATA ON TOP OF THE PROGRAM DECK AND LOAD THESE TWO 
DECKS ( SEE THE GENERAL OPERATING INSTRUCTIONS). AFTER THE PROGRAM 
DECK HAS BEEN READ* THE DATA WILL BE READ AND PUNCHED. THE PROGRAM 
MAY TAKE SEVERAL MINUTES FOR LARGE STRUCTURES. THEREFORE* THE 
PROGRAM SHOULD NOT BE STOPPED UNLESS IT IS OBVIOUS THAT IT IS 
TAKING TOO MUCH TIME. 

4. AFTER THE INTERNAL FORCES AND DISPLACEMENTS HAVE BEEN CALCULATED 
AND PUNCHED. THE MESSAGE - EMPTY PUNCH AND LOAD DATA - WILL BE 
TYPED. REMOVE THE PUNCHED OUTPUT FROM THE PUNCH STACKER AFTER 
THE PUNCH UNIT HAS BEEN EMPTIED. 

5. LIST THE OUTPUT ON THE 407 ACCOUNTING MACHINE. 

6. IT IS NOT NECESSARY TO RELOAD THE PROGRAM DECK IN ORDED TO RUN 
ANO THE SET OF DATA. PLACE 5 THE NEW DATA IN THE 1622 READ HOPPER 
AND DEPRESS READER STARTCNOT LOAD) AND PUNCH START ON THE 1622. 
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C C MATRIX FORCE /DISPLACEMENT PROGRAM 

C E.L. COOK AND G.E. LAMBERT - REVISED MARCH 1966 

C 

DIMENSION FV(50*5) *B01 (50*28) *B1TFV(50) *D01 (28 ) *FDC5t5> 

C 

C PART I - INPUT 

C 

BEGIN TRACE 

C 

C READ AND PUNCH PARAMETERS 

1 REAO lOO*K*L*M«N«NZEB0*NZEBl *KK 
GO T0C2*3)*KK 

2 PUNCH 301 
GO TO 4 

3 PUNCH 302 

4 PUNCH 30O*K*L.M.N*NZEB0*NZEBl 

C 

C REAO AND PUNCH FV OR KP BY ROWS* FIVE VALUES PER CARD 

GO TO(5*6>*KK 

5 PUNCH 311 
GO TO 7 

6 PUNCH 312 

7 DO 8 I»1*L 

REAO 101 *FV« 1*1) *FV< 1*2) *FV< 1*3) *FV( 1*4) *FV( 1*5) 

8 PUNCH 101 *FVf I* 1 > .FV( I *2 > *FV( I *3) *FVC I*4)*FV<I.5> 

C 

C CLEAR BO 1 OR AO 1 

MN«M+N 
DO 11 I"1*L 
DO 11 J«1*MN 

11 BOKI* J) "0.0 

C 

C READ ANO PUNCH NONZERO ELEMENTS IN BO OR AO 

GO TO < 12*13) »KK 

12 PUNCH 321 
GO TO 14 

13 PUNCH 322 

14 DO 15 IA«1«NZEB0 
READ 102*I*J.B01( I*J) 

15 PUNCH 1 02. I ♦ J. BO I ( I .J' 
IF(N> 16*31 .17 

16 PRINT 333' 
STOP 

C 

C READ AND PUNCH NONZERO ELEMENTS IN Bi OR A 1 

17 GO TO( I8t 19) *KK 

18 PUNCH 331 
GO TO 20 

19 PUNCH 332 

20 DO 21 IA«1 *NZEBI 
READ 102*I*J*B1IJ 
PUNCH 102* I ♦ J«8l IJ 
JM*J+M 

21 BOl ( I ♦ JM)=Bl U 
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C PART II - RECURSION ANALYSIS 

31 DO 49 IR»1 *N 

Bl TRANSPOSE TIMES FV OR Al TRANSPOSE TIMES KP 
DO 32 IM*L 

32 B1TFVC I )«0.0 
J*l 

DO 45 IS«1*K 

IF(FV< J.t> > 50* 35* 33 

33 J4*J+4 

DO 34 IA»1.5 
IAJ*IA+J~1 
DO 34 JA«J»J4 

34 BITFV< IAJ>«B1TFV< I A J>+B01 ( JA«MN)*FV( JA. I A) 
J»J+5 

GO TO 45 
33 IF(FV(J*2) >50*38*36 

36 J3-J+3 

DO 37 IA«2*5 
IAJ=!A+J-2 
DO 37 JA=J*J3 

37 B1TFV< IAJ)«B1TFVC IAJ)+60l < JA*MN)*FV< JA. I A> 
J=J+4 

GO TO 45 
36 IF(FV( J.3) >50«41«39 

39 J2*J+2 

DO 40 IA«3*5 
IAJ*IA+J~3 
DO 40 JA»J«J2 

40 BlTFVC IAJ>»B1TFV( !AJ>+BOl ( JA« MN) *FV( JA* I A) 
J»J+3 

GO TO 45 

41 IF(FV< J*4> )50*44*42 

42 Jl«J+l 

DO 43 IA«4.5 
IAjalA+J-4 
DO 43 JA«J«J1 

43 B1TFVC IAJ)*B1TFVC IAJ1+B01 < JA*MN>*FV( JA* I A) 
J»J+2 

GO TO 45 

44 B1TFVC J)«80l (J**N>*FVt J«5> 
J*J+1 

45 CONTINUE 

CALCULATION OF D0MB1TFV TIMES B01 OR A1TKP TIMES AO 1 ) 
DO 46 1 = 1 *MN 
D0l<D*0»0 
DO 46 J»l *L 

46 DO 1(1 >«001 < I )+BtTFV< J)*801 ( J, I ) 

CALCULATION OF REDUNDANTS OR DEF I C I ENC I ES C STORED IN D01 ) 
DO 47 I«1*MN 

47 DO I < I >*-D01 ( U/DOI (MN> 
CALCULATION OF NEW B01 OR A,0l 
DO 48 1*1 .L 

DO 48 J*1*MN \ 

48 BOKl • J>=B0l ( I * J *+B0 1 ( I •MN>*DOl < J> 

49 MN«MN-1 

GO TO 71 

50 PRINT 334 
STOP 
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C C PART III - CALCULATION OF FLEXIBILITY OR STIFFNESS MATRIX 

C CALCULATION OR FV TIMES B OR KP TIMES A t STORED IN B> 

71 M1«M+1 
MMsM+M 

DO 72 1=1 *L 
DO 72 J'MUMM 

72 BOH I «J>=0*0 
DO 65 IRM»M 
IRM*!R+M 
J»l 

00 85 IS»1*K 
IF<FV<J*1> J50«75*73 

73 J4»J+4 

DO 74 !A*l«5 

IAJ-IA+J-1 

DO 74 JA«J«J4 

74 B01 ( t A J * I RM y *B01 i IAJ* IRfO+BOl ( JA, IR>*FV< JA, I A) 
J»J+5 

GO TO 85 

75 IF(FV( J«2> ) 50*78* 76 

76 J4-J+3 

DO 77 IA«2«5 
IAjaIA+J-2 
DO 77 JA*J«J4 

77 BOH IAJ*IRM>*B01C!AJ«IRM)+B01CJA«1R)*FV<JA*IA> 
J»J+4 

GO TO 85 

78 IFCFV(J*3) 150*81 *79 

79 J4=J+2 

DO 80 IA«3»5 
IAJ-IA+J-3 
DO 80 JA«J#J4 

80 B01< IAJ«IRM)-B0M IAJ* J»M>+BOl(JA*!R)*FV<JA»!A> 
J«J+3 

GO TO 85 
61 IF<FV< J*4) )50«84*82 

82 J4*J+1 

DO 83 IA«4*5 
IAj*IA+J-4 
DO 83 JA-J.J4 

83 B01 ( IAJ.IRM)«B01 ( I A J* IRMI+BOl ( JA* IR!*FV< JA* I A) 
J*J+2 

GO TO 85 . 

84 B01 ( J* IRM)«B01 C J* IR)*FV< J*5> 
J=J+1 

85 CONTINUE 

C 

C CALCULATION OF FD OR KP 

DO 88 I«1«M 
IM*I+M 
DO 88 J»1*M 
FD( I . J>*0.0 
DO 88 K»l *L 

FD( I « J>-FD< I ♦J)*B01(K, IM>*B01 <K»J) 
IF ( ABS(FD< I.J) >-• 10**20) 87*87* 88 

87 FD(I*J)-0*0 

88 CONTINUE 
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PART IV - 



INTERNAL FORCES AND DISPLACEMENTS 



- DISPLACEMENT METHOD 



GO TO ( 98 * 6O0 ) » KK 

C 

C INVERSION OF KF TO OBTAIN FD 

600 DO 60« 1=1 «M 
STORE=FD( I • I > 
FD( I • I > = I «0 
DO 60 1 J=1«M 

601 FD( I * J > =FD( I »J)/STORE 
DO 604 K>1 ,M 

IF (K- 1)602. 604 ,602 

602 STORE=FD(Kt I > 
FD(K» I ) =0.0 
DO 603 J'lfM 

603 FD(K« J)=FD(Kt J>-STORE*FD( I ♦ J> 

604 CONTINUE 

C 

C CALCULATION OF B 

DO 99 I = 1 ♦ L 
DO 99 J=1*M 
BO 1 ( I iJ)=0.0 
DO 99 K-l «M 
IM=K+M 

99 SOI < ! • J> S B01 ( !« J>+B01 ( I* IM)*FD(K« J> 

C 

C FORMAT STATEMENTS 

C 

100 FORMAT (715) 

101 F0RMAT(5E16»8> 

102 FORMAT(2I5tE16.8> 

300 FORMAT (6 15/) 

301 FORMAT ( 26H FORCE METHOD - PARAMETERS/) 

302 FORM AT ( 33H DISPLACEMENT METHOD PARAMETERS/) 

311 FORMAT ( 22H ELEMENT FLEXIBILITIES/) 

312 FORMAT (20H ELEMENT STIFFNESSES/) 

321 FORMAT (/23H NONZERO ELEMENTS IN BO/) 

322 FORMAT (/23H NONZERO ELEMENTS IN AO/) 

331 FORMAT ( /23H NONZERO ELEMENTS IN Bl/) 

332 FORMA T(/23H NONZERO ELEMENTS IN Al/) 

333 FORMAT ( 1 4H N IS NEGATIVE) 

334 FORMAT ( 29H NEGATIVE FIRST ELEMENT I N FV ) 

335 FORMAT (26H EMPTY PUNCH AND LOAD DATA) 

400 FORMAT (5E16» 8) 

401 FORMAT ( / 1 6H INTERNAL FORCES/) 

402 FORMAT (/14H DISPLACEMENTS/) 
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C C PART V - INTERNAL FORCES AND DEFLECTIONS - OUTPUT 



C INTERNAL FORCES 

98 PUNCH 401 

GO T0(61 1 »62I ,631 ,641 ,651 ) »M 

611 DO 612 1=1, L 

612 PUNCH 400*SOKI«I> 
GO TO 700 

621 DO 622 J-l.L 

622 PUNCH 400*B0l ( I * l > *Q01 ( I <2) 
GO TO 700 

631 DO 632 1*1 ,L 

632 PUNCH 400 • BO 1(1,1), BO 1(1*2), BO 1(1,3) 
GO TO 700 

641 DO 642 I*I,L 

642 PUNCH 400,B01 ( I , 1 ) ,B01 ( I ,2) »B01 ( I ,3) ,B01 ( 1,4) 
GO TO 700 

631 DO 652 I*1,L 

652 PUNCH 400 *80 1(1*1) *BQ 1 ( I ,2 ) • BO 1 ( I • 3 ) ,80 1 ( I * 4 ) • BO 1 ( I • 5) 

C 

C DISPLACEMENTS 
700 PUNCH 402 

GO TO ( 7 1 1 t 72 I • 73 1 , 74 1 , 75 1 ) • M 

711 DO 712 1*1, M 

712 PUNCH 400*FOCI*1) 
GO TO 800 

721 DO 722 1*1, M 

722 PUNCH 400* FD (1,1), FD (1,2) 
GO TO 800 

731 DO 732 1*1 ,M 

732 PUNCH 400«FD< I • X > *FD< I *2> *FD< I «3> 
GO TO 800 

741 DO 742 I*1,M 

742 PUNCH 400,FD( I • 1 ) ,FD( I ,2),FD( I ,3) »FD( I ,4) 
GO TO 8O0 

75 1 DO 752 I * I , M 

752 PUNCH 400, FD( I * 1 ) ,FD( I *2>,FD( I ,3) ,FD( I ,4) ,FD( I ,5) 

C 

800 PRINT 335 
GO TO 1 

C 

END TRACE 

C 

END 
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DISPLACEMENT METHOD - PARAMETERS 

+ 3 +5 +2 + 2 +6 +A 



FORCE METHOD - PARAMETERS 

+3 +5 +2 +1 +3 +4 



ELEMENT STIFFNESSES 



ELEMENT FLEXIBILITIES 



♦•0O0OOOO0E-50 
+ • 0O00OOO0E-50 
♦•OOOOOOOOE-50 
♦•0C00C000E-50 
+•000000006-50 



+.0OO0OO00E-50 

+.0000000CE-50 
♦•0000000CE-50 
♦•OOOOOOOCE-50 
+♦000000005-50 



♦•OOOOO000E-50 
♦•OOOOOOOOE-50 

+«0O000OO0E-50 
+.OOO0OOOOE-50 
+.OOO0OOOOE-50 



♦ • 66666666 E+ 07 

♦•33333333E+07 
+ • 13333333E+03 
+ • 66666666E+07 
+.O0O0OOO0E-50 



♦•33333333E+07 
+ « 66666666E+0 7 
i .&6666666E+07 
+• I3333333E+08 
♦•50000000E+Q4 



+♦ 000OOO00E-50 
+.00000000E-50 

+*0GOQ0OO0E-50 
+.000OOO00E-50 
+.OO0OOOO0E-50 



♦•O000OO00E-50 
+.0OO0OOO0E-50 
+.000OOO00E-50 
+»0OO0OOO0E~50 
+»OOOO0O00E-50 



♦•O0O00OO0E-50 
+.0O0OO0O0E-50 
+.OOOOOOO0E-50 
+.00000O00E-50 
+#OOO00Q0OE-50 



♦•20000000E-06 
1OOO00OOE-06 
+.100OOQOOE-06 
-•50OOOOO0E-07 
+.OOOO00O0E-50 



-•iOOOOOOOE-06 
+♦ 2000O0O0E-06 
-•5000OO00E-07 
+• 10000000E-06 
+.20OOOOO0E-03 



NONZERO ELEMENTS IN AO 



NONZERO ELEMENTS IN BO 



+ 1 


+ 1 


+ .33333334E-0 i 


+ 1 


+ 1 


♦♦30000000E+02 


+ 2 


+ 1 


+.33333334E-0 1 


+ 1 


+2 


10000000E+01 


+ 2 


+2 


+• 10000000E+0 1 


+ 2 


+2 


+• lOOOOOOOE+Ol 


+ 3 


+2 


+• 1OOO0OCCE+0 1 








+ 3 


+ 1 


-.66666667E-01 


NONZERO ELEMENTS IN Bl 


+ 4 


+ 1 


- . 66666667E-0 1 


+ 1 


♦ 1 


-•45000000E+02 


DNZERO ELEMENTS IN A I 


+ 2 


+ 1 


+ • I5000O00E+02 








+ 3 


+ 1 


150OOOO0E+02 


+ 3 


+ 1 


♦•66666667E-01 


+ 5 


+ 1 


+.10000000E+01 


+ 4 


+ 1 


+ • 66666667E-0 1 








+4 


+2 


+. lOOOOOOOE+Ol 


INTERNAL 


FORCES 


+5 


+ 1 


+ • IO0OOOO0E+0 1 









INTERNAL FORCES 

♦♦12445820E+02 
+.58513930E+01 
-•58S13927E+0I 
+.57176470E-06 
+.39009287E+00 

DISPLACEMENTS 

♦♦57120743E-04 
-.1978328 IE -05 



♦•3095960CE-02 
+.66563470E+00 
+.33436528E+00 
-.875294 10E-07 
-•22291022E-01 



-.19783282E-05 
19876 16IE-06 



+ . 1 2445820E+02 ♦ • 30960000E-02 

+ .58513932E+01 ♦♦66563470E+00 

-•58513932E+01 +»33436533E+00 

. f •OC'JOOOOOE-50 +.00000O00E-50 



+.39009288E+00 
DISPLACEMENTS 



-• 22291 022E-01 



+.57120742E-04 -• 19783282E-05 
19783281E-05 +• 19876162E-06 



o 



o 



c 



o 



COMPUTER DESIGN 



OF 



HEAVY MULTISTORY RIGID FRAME 



INDUSTRIAL STRUCTURES 



By 

G.L. Kirby 
N.M. Dudnikoff 
V.L. Arndt 



Compliments 



of 

Stearns* 



per 



Presented at the 
Western Region Summer COMMON Meeting 
Denver, Colorado 
July 7, 1966 



10 



TABLE OF CONTENTS 

Page 

Introduction . I 

Ore Storage Bin Support Design . i 

Design Considerations . 2 

Rigid Frame Supporting Structure .......... 2 

Preliminary Design 2 

Frame Layout • 2 

Loads ....... . 3 

Center Foundation Size 3 

Scheme I - Top Girder Pin-Connected to Column . . 3 

Scheme I - Data Input Sheets . 4 

Scheme I - Analysis Results . 10 

Scheme II - Column Pin-Connected to Top Girder 15 

Comparison of Results 16 

Final Design 16 

Final Member Sizes 1 5 

Final Computer Run Ii3 

Miscellaneous Details . 16 



- i i - 



LIST OF FIGURES 



Figure Page 

1. Elevation of 10,000 Ton Ore Storage Bin 1 

2. Plan of Ore Storage Bins 2 

3. Plan Showing Rigid Frame Layout and Bin Floor Framing .... 3 
^. Frame and Loads - Scheme I 

5. Moments/Reactions for Combined Vertical and Horizontal Loads 

Scheme I ........... 15 

6. Moments/Reactions for Combined Vertical and Horizontal Loads 

Scheme II . . 15 

7. Elevation of Frame Showing Final Design Member Sizes .... 16 

8. (a) Top Girder 17 

(b) Column Splice 17 

(c) Girder to Column Connection 17 

(d) Base Plate Detail 17 



/ L /6~ 



LIST OF TABLES 



Table Page 

1. Data Sheets for Support Structure - Scheme I ....... . 5 

(a) Frame Parameters ......... . 5 

(b) Geometrical Properties 6 

(c) Beam Loads 7 

(d) Column Loads . 8 

2. Computer Results for Support Structure - Scheme I .... . II 

(a) Member End Rotations/Deflections II 

(b) Member End Moments/Shears - 13 



-iy- 



Introduction. « During 1964, a Southwestern Copper Plant, undertook a major improvement and 
plant expansion program. Stearns-Koger Corporation, tngi neers-Con tractors, of Denver, Colorado 
designed and constructed the necessary facilities for this expansion. 

The facilities consisted of: 

1. Crushing plant. 

2. Expansion of concentrator plant. 

3. One new reverbatory furnace. 

H. Associated conveyors and transfer towers. 

5. Two new I0 T 000 ton live storage copper ore bins. 

6. Modification of two existing 5,000 ton live storage bins. 

Ore storage bin support design. « This paper will cover the application of Stearns-Roger • s 
Multistory Rigid Frame Program for analyzing and opt imum i zing the design of the two ore storage 
bin support frames. Figure I shows a typical bin with head house and support structure. 
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Figure !. Elevation of 10,000 Ton Ore Storage Bin. 
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Design considerations. — Reference should be made back to Figure I to fully comprehend the 
following design criteria: 

1. Bin capacity - 10,000 ton established by customer. 

2. Bottom elevation of bin established by both railroad clearance and arrangement of 
chutes used for loading railroad cars. 

3. Bins were to be located between two existing bins as shown in Figure 2. 




Figure 2. Plan of Ore Storage Bins. 



4. Top elevation of bin established by 10,000 ton capacity, since all other dimensions 
were set. 

5- Most column depths were limited by railroad clearance. 

6. Railroad traffic had to be maintained to an ore concentrator located 7 miles away. 
It was mandatory that foundations, supports, and bins be constructed during normal 
rail traffic hauling 30,000 tons of ore per day (20 trains of 30 cars each). 

Rigid frame supporting structure. — A rigid frame design was selected because: 

1. Flexibility provides better resistance to seismic activity (UBC-Zone I ). 

2. Flexibility provides better resistance to shock loading from blasting done inside the 
bins to break stoppages in plugged ore chutes. 

3. Transverse bracing would limit accessibility to operating platforms located within 
the structure. 

A. Steel construction was selected to speed erection. Aiim steel was utilized to 
reduce member (especially column) size and weight, thereby minimizing handling 
problems during erection. 

B. Due to the limited railroad clearance, pier size was kept to a minimum by using 
5000 psi concrete. 

4. The bin floor was constructed of cast in place concrete. The floor was connected to 
the top steel girder, thereby obtaining economy and rigidity from the composite con- 
struction. Floor beams were made continuous over supports where possible, for eco- 
nomical reasons, and also to decrease deflections. 

I . Prel imi nary Des i gn 

Frame layout. « Figure 3 shows the layout of the bin floor framing and the rigid frame sup- 
ports. Frames A and C were analyzed on the IBM 1620 Computer. The results from these runs were 
ratioed to provide design moments and forces for other frames. The analysis and design of Frame C 
is covered in the remainder of this paper. Frame A was analyzed and designed similarly to Frame C. 
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Figure 3- Plan Showing Rigid Frame Layout and Bin Floor Framing. 



Loads. — Vertical loads were distributed to the rigid frames from associated bin floor framing, 
and by assuming contributary widths over which distributed loads act. Dead and live loads amounted 
to almost 10 ksf. Zone I seismic and wind produced about the same lateral loading. A total lateral 
load of 230 kips was used for design. Foundations were located on rock with a 25 ksf allowable 
bearing pressure. 

Center foundation size. — Since railroad traffic could not be interrupted during construction, 
it was necessary to limit the size of the base pads for center foundations. It was also desirable 
to make the top girder continuous over the center column in order to minimize stresses and deflec- 
tions in the bin floor. However, it was assumed that this configuration would produce greater loads 
on the center column than if the top girders were interrupted at this point. These two schemes were 
investigated on the computer to determine which design was more suitable, considering girder stress/ 
rigidity vs foundation loads. 

Scheme I * top girder pin-connected to column. — Figure u. shows the ideal i zed frame which was 
analyzed on the computer. Loads, dimensions, and assumed moments of inertia are shown. At this 
stage of design, column bases were assumed pinned. Vertical and horizontal loads were assumed to be 
load conditions I and 2 respectively. 
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Figure 4. Frame and Loads - Scheme I* 
Scheme I - data input sheets. — The completed data sheets are given in Table I. 
Table I (a) describes the structure. 

Table I (b) gives structural properties - geometry, stiffness, connection rigidities, 
etc. 

Table I (c) describes beam loads. 
Table I (d) describes column loads. 
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Table I (a). Data Sheet for Support Structure - Frame Parameters. 
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Table I (b). Data Sheet for Support Structure - Geometrical Properties. 
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Table I (c). Data Sheet for Support Structure - Beam Loads. 

IBM DATA FORMAT SHEET - BEAM DATA SHEET NO.i 0F_£ 
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Table I (d). Data Sheet for Support Structure - Column Loads 
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Table i (d). (Continued) 
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COLUMN 
LENGTH 
» L(FT) 
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R 


UNiF. LOAD 
W (K/FT) 
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LOAD 
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POSITION 
0,(FT) 
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LOAD 
H 2 (K) 


POSITION 
o 2 (FT) 
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LOAD 
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POSITION 
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CONCEN 
LOAD 
P 2 (K ) 
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b 2 (FT) 


position 
c 2 (fti 




COL. 
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NO 

L_l 






SUM 
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(K) 
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FMB (K-FT) 


BOTT FIXED 
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SEE 
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Scheme I - Analysis Results. — The data sheets were keypunched and the problem analyzed on 
the computer. Table 2 (a) gives member end rotations/deflections for each load condition. Table 
2 (b) gives member end moments/shears. The final moment/ reaction diagram is given on Figure 5. 
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JOB NO, B-27500 CUSTOMER- PROJECT- 

STE ARNS-ROGER CORPORATION MULTISTORY RIGID FRAME PROGRaM - G7 

SUBJECT- BIN SUPPT. FRAME C- SCHEME I 8Y NMD DATE 5-6-64 



MEMBER ROTATIONS AND DEFLECT I QNS 

MEA 1 ?M?h.M. -^ftk 

MEA 2 -.00069 ff^^fv -.00006 " 

MEA 3 0.00000 D*$L 0.00000 



MEA 


4 


-.00069 


-.00006 


MEA 


5 


0.00000 


^-.00083 


MEA 


6 


0.00000 


0.00000 


MEA 


7 


0.00000 


-.00083 


MEA 


8 


.00069 


-.00006 


MEA 


9 


0.00000 


0.00000 


MEA 


10 


.00069 


-.00006 


MEA 


11 


.00034 


-.00002 


MEA 


12 


0.00000 


0.00000 


MEA 


13 


-.00002 


-.00039 


MEA 


14 


.00007 


-.00069 


MEA 


15 


0.00000 


0.00000 


MEA 


16 


.00007 


-.00069 


MEA 


17 


0.00000 


-.00083 


MEA 


18 


0.00000 


0.00000 
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l 2. 



MEA 


19 


0.00000 


-.00083 


MEA 


20 


-.00007 


-.00069 


MEA 


21 


0.00000 


0.00000 


MEA 


22 


-.00007 


-.00069 


MEA 


23 


.00002 


-.00039 


MEA 


24 


0.00000 


0.00000 


MEA 


25 


-.00034 *p frt. 


-.00002 


MEA 


26 


-.00002 MK Rrt* 


-.00039 


MEA 


27 


.03172 Pe$L 


.06168 


MEA 


28 


-.00069 


-•00006 


MEA 


29 


.00007 


-.00069 


MEA 


30 


.03172 


.06168 



2 2 



Z4 



k VS 



I 

* 



MEA 


31 


0.00000 


-.00083 


MEA 


32 


0.00000 


-.00083 


MEA 


33 


.00001 


.10957 


MEA 


34 


.00069 


-.00006 


MEA 


35 


-.00007 


-.00069 


MEA 


36 


-.03171 


.06164 


MEA 


37 


.00034 


-.00002 



TabU 2 (a). Member End AoUtle*t/D«f tactions. 
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MEA 


38 


.00002 


-.00039 


MEA 


39 


-.03171 


.06164 


M EA 


40 


-.00002 


-.00039 


MEA 


41 


.00035 


-.00157 


MEA 


42 


-.03168 


.16245 


MCA 


43 


•00007 


-.00069 


MEA - 


44 


.00030 


-.00141 


MEA 


45 


-.03168 


.16245 


MEA 


46 


0.00000 


-.00083 


MEA 


47 


0.00000 


-.00083 


MEA 


48 


.00001 


.11457 


MEA 


49 


-.00007 


-.00069 


MEA 


50 


-.00030 


-.00141 


MEA 


51 


.03175 


.16249 


MEA 


52 


.00002 


-.00039 


MEA 


53 


-.00035 


-.00157 


MEA 


54 


.03175 


.16249 



If 



11 
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It 
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J* ft. 



. 6* I&4 



T.bl. 2 (a). (Continued) 
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JOB NO. 8-2.7300 CUSTOMER- PROJECT- 

STEARNS-ROGER CORPORATION MULTISTORY RIGID FRAME PROGRAM - <\'{ 

SUBJECT- BIN SUPPT. FRAME C- SCHEME I BY NMD DATE 5-6-64 



FINAL MOMENTS AND SHEARS u „j c*nA.Z 

EMS 1 v -,21 1****1, -70.59 

EMS 2 -1801.10 -91.76 |( 

EMS 3 -47 1* 67 faS^f -11.59 

EMS 4 1945.50 -47.43 

EMS 5 0.00 0.00 \Z 

EMS 6 -392.22 -2.37 

EMS 7 0.00 0.00 

EMS 8 -1945.50 -47.37 |3 

EMS 9 -587.78 -2.36 



EMS 10 1801.10 -91.67 

EMS 11 .14 -70.54 

EMS 12 -214.34 -11.58 



14 



EMS 13 2.55 -127.40 

EMS 14 11.03 -153.78 K 

EMS 15 .97 -20.08 

EMS 16 0.00 0.00 

EMS 17 0.00 0.00 

EMS 18 0.00 0.00 

EMS 19 0.00 0.00 

EMS 20 0.00 0.00 23 

EMS 21 0.00 0.00 

EMS 22 -11.07 -153.76 * 

EMS 23 -2.59 -127.38 Z4 c& 

EMS 24 -.97 -20.08 

EMS 25 .20 *r-w*>*. 70.59 

EMS 26 24.23 bitt***"' 43.86 

EMS 27 2.22 **«*C 10.40 £ 

EMS 28 -144.12 139.20 f -2 

EMS 29 42.55 -13.38 11 ^ 

EMS 30 -9.23 11.43 

EMS 31 0.00 .03 

EMS 32 0.00 .01 \J 

EMS 33 0.00 0.00 

EMS 34 144.11 139.07 

EMS 35 -42.61 -13.52 |4 

EMS 36 9.22 11.41 

EMS 37 -.20 70.54 i«r 



Table 2 (b). Member End Momenti/ Shears. 



A57 



EMS 


38 


-24.2 3 


43.01 


EMS 


39 


-2. 22 


1 0.39 


EMS 


40 


-26.79 


83.53 


EMS 


41 


0.00 


-.01 


EMS 


42 


-2. 32 


7.26 


EMS 


43 


-53.57 


167. 16 


EMS 


44 


0.00 


-.04 


EMS 


45 


-4.65 


1 4. 54 


EMS 


46 


0.00 


.03 


EMS 


47 


.01 


.05 


EMS 


48 


0.00 


0.00 


EMS 


49 


53.66 


167.34 


EMS 


50 


-.02 


• 11 


EMS 


51 


4.66 


14.56 


EMS 


52 


26.83 


83. 58 


EMS 


53 


0.00 


.02 


EMS 


54 


2.33 


7.26 



lb 
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Table 2 (b). (Continued) 



o 



15 



Scheme 11 - column p in-connected to top girder. — All cards required to reflect continuity 
over the center column were keypunched and inserted in the original data deck, and the problem 
reanalyzed on the computer. Figure 6 shows the idealized frame and also the final moment/ react ion 
diagram for the computer rerun. 




Figure 5. Moments/ React ions for Combined Vertical and Horizontal Loads - Scheme I. 




Figure 6. Moments/ React ions for Combined Vertical and Horizontal Loads - Scheme II. 
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Comparison of results. — For combined vertical and lateral loading, it can be observed that 
Scheme II produced lower moments in the top girder (I7U7 ft-kip vs 1993 ft-kip for Scheme l). 
However, Scheme I produced a considerably lower center column reaction (785 kip vs 1010 kip for 
Scheme ll). 

2- Final Design 

Final member sizes, — Members were sized, based on the computer results for Scheme I, be- 
cause of foundation size limitation. With the exception of the center column, base plates and foun 
dations were designed for an assumed 50% fixity in order to reduce wind caused sway. The center 
column was designed as pinned at top and bottom, so that only axial load could be transmitted. 
Any column fixity would have produced overturning moments, necessitating larger center column 
foundations. All members were fabricated from plate. The steel was fabricated in Corpus Christi, 
Texas. Columns were stress relieved according to ASME Specifications. They had to be shipped to 
the nearest stress relieving furnace at Houston, Texas. There are approximately 250 tons of steel 
in each bin support (including floor framing), most of this is A44I steel. 

No recognized standard procedures covering the design of the composite top girder could be 
found in the building codes, or available literature. This was due mainly to the massive size of 
the girders, and also to the heavy loading pattern. Much guidance and help was received from the 
American Institute of Steel Construction in designing shear connections for these girders. 

Figure 7 shows the final frame with member sizes noted. 




Figure 7. Elevation of Frame Showing Final Design Member Sizes. 

Final computer run. — Since Scheme I computer results were predicated upon assumed member 
stiffnesses and pinned column bases, the decision was made to rerun the frame using actual member 
stiffnesses. The purpose of this run was to give a final design check on the frame and to produce 
computer results compatible with the final design. All member sizes were found to be acceptable 
and no design changes were made. 

Miscellaneous details. — Figure 8 shows various details of the frame. 
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Figure 8 (d). Base Plate Detail, 
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LIBRARY, MAGAZINES, AND A COMPUTER 



By Gale Ahlborn 
Assistant Research Engineer 
Institute of Transportation and Traffic Engineering 
University of California, Berkeley 

There has been a considerable amount of information published concerning libraries 
and automation, the term automation referring to anything from manual manipulation of 
punched cards to the use of computers with disk and tape files and remote communication 
equipment. We were interested in the use of computers in library operations but were 
unable to find out what has actually been accomplished or completed and whether these 
operations have been beneficial or not. There is little doubt that large libraries can make 
good use of computer equipment in their operations, but there has been some question 
about the use of computers with smaller libraries. A few examples of the use of com- 
puters with libraries are: Key-Word-In-Context indexes, catalog card production, pro- 
ducing book catalogs, circulation control, lists of recent acquisitions , and operations with 
library serials. 

From the standpoint of making use of a computer, library activities can be generally 
considered to be the manipulation of many files of information — a searching and sorting 
of various items of information for the purpose of performing some control on a library 
operation, producing an information list, or a combination of the two. These various li- 
brary information files are of varying size and, most important, require updating at various 
intervals of time. The frequency of updating is a function of the particular file and the opera- 
tion being performed with this file. As an example, if the operation is control of circulation, 
that is control and information concerning publications being checked in and out of the library, 
this file requires updating at least daily if not immediately as each operation occurs. On the 
other hand, a file concerning annual publications can be updated at a much less frequent in- 
terval, possibly only once per year. 

Each computer installation differs from another in equipment, size,, personnel, and 
operation. The same applies to libraries. Each library is different in size, equipment, 
personnel, and operation. 

Our particular computer equipment is a 1620 Model I system, twenty thousand posi- 
tions of core storage, one disk drive with punch card input and output. Supporting off-line 
equipment for this system are key punches, sorter, and a 407 accounting machine for 
printing. Our partner in this operation is the Transportation Engineering Library, a 
special purpose library. Both computer and library are located at the Richmond Field 
Station about five miles from the Berkeley campus of the University of California. The 
primary use for this computer is educational. It is used by scheduled engineering classes 



and graduate students in their advanced degree research. The development of the use of 
this equipment in the library's operation has had a low priority and therefore has progressed 
rather slowly. 

Our first venture into this kind of library automation has been to work with the 
library's serial file. Since I am a library user rather than a librarian, I am substituting 
the word "magazine" in place of "serial, " although I have been informed that this is not 
entirely correct. There are various items of interest concerning magazines, such as, 
the frequency of issue, the general subject matter of each magazine, the cost of sub- 
scription, the length of subscription,, to whom the magazine is circulated, and in what 
order. At the end of each year some magazines are discarded, others are held a certain 
number of years, while others are retained indefinitely. For those magazines that are 
retained, there are other items of interest, such as, the number of issues to be bound 
together, who does the binding and what type of binding is used. The magazine file of 
the Transportation Engineering Library consists of about five hundred and fifty maga- 
zines, and there are about thirty possible items concerning each magazine that are of 
interest to the library. 

When we first started, we considered developing a rather complete system of 
computer operation with the magazine file. This system would have the magazine file 
on the disk storage at all times, available for on-line updating and producing various 
reports to provide information and control of the library's magazine operation. Although 
we felt this system had merit, we were concerned as to whether it would be justified for 
a small computer-small library operation and also concerned about the lack of experience 
on both sides in this area of library automation. We finally decided that it would be better 
to develop a manual system of operation first, use this system for a period of time, re- 
view our experience, then determine whether to continue with the same or a different 
manual system or to switch over to a complete computer operation. 

The manual system consisted of having the magazine file punched into a deck of 
cards, then using the mechanical sorter and accounting machine to produce the desired 
reports. With this system, the computer was to be used only as necessary to supple- 
ment the operation. Three cards were used for the information concerning each maga- 
zine. The first two cards contained information that was of interest to both the library 
and library users. The information on the third card was only of interest to the library. 
Numeric codes were duplicated into a common area on all three cards for the purpose of 
sorting. The information punched into the cards outside of the sorting area was either 
in plain language or a form of alphabetical coding so that the information is relatively 
easy to understand when printed directly on the accounting machine. 
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Five different reports are produced using this manual system. The largest 
report is the "Users Report" that contains nine different listing from the magazine file. 
These lists include an alphabetical list, newsletters, statistical sources, foreign 
periodicals arranged by country of origin and language, subject index, indexing and 
abstracting services. In all, the report is sixty-five pages in length. Since there are about 

250 copies of this report distributed, offset masters for the press are printed directly 

j 

from the cards on the accounting machine. In addition to the five reports presently pro- 
duced, there are at least five more that are desired which we are hoping to avoid producing 
until we revise our present magazine file system. 

We have used this manual system for approximately a year with no major changes 
being made to the system. Several items of information about the magazines were de- 
leted from the file, while other items were added. It became necessary to make use of 
the computer almost from the very beginning of this operation for several reasons. 
First, the computer is very efficient at checking the validity of codes punched in the 
cards. We use both the computer and a manual checking of the list to detect errors. 
Second, we found that it was desirable to use various items of information in various 
combinations which would involve a fairly complex mechanical sorting. The computer 
was used quite easily for doing this type of operation. 

We found that updating of the magazine file consisted mostly of making many 
small changes throughout the file rather than having to make large changes. The frequency 
of updating with this file is dependent on the frequency of producing reports. Considering 
all the various types of reports that our library would desire from this file* it would not 
be necessary to produce reports any more frequently than once per month. 

Our conclusion, after working with this manual system, is that the magazine file 
should be changed over to a computer system. The system that we have devised and 
is being developed at the present time is as follows: The magazine file will be retained 
on cards in the library where they can be easily referred to and updated as needed. 
When reports are required, the card file will be loaded onto the disk working area and 
then the reports generated. While this system sounds very similar to what we had 
originally considered, there actually are large differences between the two systems. We 
know now that if we had developed the original system, that there would have been considerable 
difficulty with it, some portions being practically unworkable. We fee now that we are 
developing a very practical and workable system that will have sufficient flexibility. 
This system will set the pattern for the additional work we expect to do with the library. 
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Accurate Solution of 
Systems of Linear Equations 



This paper is concerned with the solution of the matrix 



equations 



Axsb, 



(1) 



AD* I, 



(2) 



where A is a general n x n matrix, assumed to be non-singular. 

Eqn. (l) is the matrix formulation of the system of n equations 
in n unknowns 



Eqn* (2) is the problem of finding the inverse DsA" of A, 

Although the solution of (1) is given by xrA^b, methods 
of computing x exist which involve less calculation than the cal- 
culation of A"**- and multiplication by b* On the other hand, the 
solution of Eqn. (2) itself may be treated as an extension of 
Eqn. (1). Again, however, there are methods of computing A~* 
which require less calculation than this extension. Thus optimum 
efficiency requires separate methods of dealing with these two 
problems, although there are many operations which are basic to 
both, (Considerations of efficiency will, in this paper, be almost 
entirely from the standpoint of the number of calculations --and 
thus the time—required to obtain the solution to a given degree 
of accuracy, rather than in terms of storage requirements • ) 

Even if a method is chosen which is algebraically best suited 
to the given problem (solution of systems of equations or the inverse 
problem), it may not jgive the accuracy desired in a particular 
application, because of the nature of the matrix A — i.e. A may be 
almost singular, or "ill -conditioned." Worse still, the results 
calculated may offer no indication of the condition of A, although 
in physical applications it is often important to know how close 
to singular A is, so that the effect of imprecision in the elements 
of A can be estimated. If A is ill-conditioned then even slight 
imprecision in the elements of A may mean gross imprecision in the 
solution. 

A standard method of obtaining high accuracy in computer 
solutions of matrix equations is to carry out all calculations in 
double precision. However, this requires about four times as much 
calculation as single precision and twice as much storage. Also 
no knowledge of the condition of A is obtained, and the amount of 
calculation is the same whether A is well- or ill-conditioned. 





n. 



4- 1 




A better procedure than blanket use of double precision appears 
to be to calculate by a direct method an initial approximation to 
the solution in single precision, and then improve this solution by 
a method of successive refinement • The latter method does require 
some double precision, but it turns out that only double precision 
additions are required, and none of the elements of the matrix A 
or the solution need be stored in double precision. 

With this method the number of iterations and thus the amount 
of calculation required depends on the condition of A and will be 
less for a well-conditioned matrix* In addition, the number of 
iterations does give a valuable indication of the condition of A. 

A minor modification to this procedure is to do some of the 
calculation in the direct method in double precision so that a 
better initial approximation is obtained • In this case, for a given 
matrix fewer iterations will be required. Again, accuracy can be 
significantly improved using only double precision additions, if 
a proper ordering of the calculations is made. 

A direct method which accomplishes this ordering is the Crout 
method, one of several "compact" schemes, so-called because they 
require no intermediate modification and restorage of the elements 
of A. 

This paper will not go into the mathematical details of either 
the Crout method or the iterative improvement method; rather, its 
remainder will be devoted to a general description of the Fortran 
II subroutines which have been written to implement this method on 
the IBM 1620 computer • Special attention will be given to the sub- 
program which is the heart of the method: the DOT function sub- 
program which is coded in SPS to do the necessary double precision 
arithmetic. Using these routines solutions can usually be obtained 
which are correctly rounded to single precision (eight digits), even 
when A is very ill-conditioned. 

The mathematical details of the methods used are contained in 
a thesis which the author wrote at Austin College. The essential parts 
of this thesis will be included in the program writeupwhen the set of 
routines is submitted to 1620 program library. The ideas behind these 
routines and the thesis were, of course, not original. The literature 
on the subject is extensive. Primary stimulus for the project came 
from a report by Cleve Moler of the Jet Propulsion Laboratory in 
California, entitled "Numerical Matrix Inversion" (JPL Technical 
Report No. 32-39U, March 15, 1963). 




General description of the Fortran II programs 
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The structure of the programs was designed to make them flexible 
and convenient to use without making it necessary to waste time cr 
storage by providing a more general program than is necessary for 
a particular application . 

Thus the implementation has a hiearchial structure, being com- 
posed of nine subprograms, some of which are used primarily to 
organize the others for the special task to be performed. The 
basic subprograms and their uses are as follows t 



FACTOR: factor matrix A into product LU using 
Grout method with partial pivoting* 

SOLSYS: solve Lya b and Ux=y. 

SIMFRV* improve the accuracy of the solution 
vector x (uses SOLSYS ) • 

INVA; form inverse D of A from L and U. 

IMPRUV: improve the accuracy of the inverse D. 

INOUT: input and/or output a matrix. 



These basic routines may be used directly by the programmer in 
his own Fortran II program, or more conveniently by use of one 
the following subprograms J 



SOLVE; solve Ax*b for x (uses FACTOR, SOLSYS, and 
SIMPRV). 

INVERT*, find inverse D of A (uses FACTOR, INVA, and 
IMPRUV). 



Finally, a main program is provided which can input a matrix, 
perform either the inversion or solution of equations, or both, 
and output the results. The basic subprograms are used directly 
instead of SOLVE and INVERT however, to avoid repetition of the 
factorization of A. 



DOT : 



accumulate inner products in double 
precision and round the result to single 
precision (used by all of the basic 
routines below except INOUT), 




DOT function subprogram 

The methods used in the Fortran routines were chosen so that 
all of the critical arithmetic would be concentrated in a single 
type of operations the computation of inner or M dot fl products of 
vectors t 



In practice there are two ways of computing such stuns of products* 

(1) Round (or truncate) to single precision after each multipli- 
cation and each addition. 

(2) Retain the sum of products in a double length accumulator 
until the entire dot product has been formed, then round to 
single precision. 

It should be obvious that the amount of round-off error incurred 
would be much less in the latter method than the first. Yet only 
double precision additions are required, since the multiplications 
are of single precision operands (giving a double length product 
as with any single precision multiplication). 

The special function subprogram called DOT was written for the 
1620 to compute inner products using the latter method of computation* 

The structure of DOT was patterned after that of a routine 
written for the IBM 7090, described by Moler. It is structured so 
that it is possible that one or both of the vectors of which the 
dot product is to be formed may be a row or column of a matrix. 
(See Figure 1 for an example of the use of DOT to perform a 
matrix multiplication.) Provision is made for adding a single 
number to the accumulated product (actually the accumulator is 
set equal to this number before the product accumulation begins). 
If the number of terms in the summation is specif ied to be zero, 
then this single number becomes the value of the summation. These 
features make the DOT suprogram a valuable programming aid, since 
a single use of it can replace several statements (as in Figure l), 
and the special case where the number of terms in the summation is 
zero is automatically taken care of. 

All of these structural features can be implemented by a simple 
subprogram in the Fortran II language. In most Fortran II systems, 
however, double precision arithmetic is not allowed. This was the 
case with the 1620 Fortran II system, even though the 1620 is a 
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variable word length and the Fortran II system does to an extent 
take advantage of this feature. In 1620 Fortran II programs it is 
possible to vary the length of the floating point word mantissas 
from 2 to 28 digits, but unfortunately this variation is permitted 
only from program to program and not with a single program or set 
of subprograms. That is, if 8 digit mantissas are desired for 
most numbers stored by the the program, there is no way to use any 
numbers with longer mantissas. Thus a special machine language 
routine had to be written to implement DOT on the 1620. 



Conclusion 

The DOT subprogram and other subprograms have been extensively 
tested at the Austin College Computing Center. These tests bear out 
what is mathematically evident; that the use of the DOT' method of 
accumulating inner products contributes significantly to the 
accuracy obtainable from the Crout method, and is essential to the 
proper convergence of the iterative improvement scheme. 

Thus the DOT function would seem to be a very useful tool to 
have available for matrix work, where solutions of known high 
accuracy are desired. In fact the DOT function potentially has 
many other uses, in any calculations in which the critical arith- 
metic can be concentrated in the calculation of inner products. 
As Moler points out, for example, the complete eigenvalue problem 
is thus formulated by Householder's method. 



Algebraic notation; 




Fortran coding: 



10 SUM* SUM*A(I,K)*B(K,J) 
C(I,J)* SUM 



SUMaO. 
DO 10 K*1,N 



Fortran with DOT function t 



C(I,J)« D0T(N,A,If MA,MA,B,14 J*MA,1,0.) 



* Figure 1 * 




